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ABSTRACT 

Quasar absorption lines provide a precise test of whether the fine-structure constant, a, is the 
same in different places and through cosmological time. We present a new analysis of a large 
sample of quasar absorption-line spectra obtained using UVES (the Ultraviolet and Visual 
Echelle Spectrograph) on the VLT (Very Large Telescope) in Chile. We apply the many- 
multiplet method to derive values of Aa/a = (a^ — ao) /ao from 154 absorbers, and combine 
these values with 141 values from previous observations at the Keck Observatory in Hawaii. In 
the VLT sample, we find evidence that a increases with increasing cosmological distance from 
Earth. However, as previously shown, the Keck sample provided evidence for a smaller a in 
the distant absorption clouds. Upon combining the samples an apparent variation of a across 
the sky emerges which is well represented by an angular dipole model pointing in the direction 
RA = (17.3 ± 1.0) hr, dec. = (-61 ± 10) deg, with amplitude O.QTto^^l x 10"^ The dipole 
model is required at the 4. Icr statistical significance level over a simple monopole model 
where a is the same across the sky (but possibly different to the current laboratory value). 
The data sets reveal remarkable consistencies: /) the directions of dipoles fitted to the VLT 
and Keck samples separately agree; //) the directions of dipoles fitted to z < 1.6 and z > 1.6 
cuts of the combined VLT-nKeck samples agree; ///) in the equatorial region of the dipole, 
where both the Keck and VLT samples contribute a significant number of absorbers, there is 
no evidence for inconsistency between Keck and VLT. The amplitude of the dipole is clearly 
larger at higher redshift. Assuming a dipole-only (i.e. no-monopole) model whose amplitude 
grows proportionally with 'lookback-time distance' (r = ct, where t is the lookback time), the 
amplitude is (1.1 ± 0.2) x 10^^ GLyr^^ and the model is significant at the 4.2(7 confidence 
level over the null model (Aa/a = 0). We apply robustness checks and demonstrate that the 
dipole effect does not originate from a small subset of the absorbers or spectra. We present an 
analysis of systematic effects, and are unable to identify any single systematic effect which 
can emulate the observed variation in a. To the best of our knowledge, this result is not in 
conflict with any other observational or experimental result. 

Key words: quasars: absorption lines, cosmology: observations, methods: data analysis 



1 INTRODUCTION 



shifts (z), certain UV metal transitions are redshifted into the opti- 
cal and therefore are observable with ground-based telescopes. The 
wavelengths of these transitions depend to various degrees on the 
fine-structure constant, a = e^/(47reoftc). By comparing observa- 
tions of these transitions with precision laboratory measurements, 
one can determine if the fine-structure constant is the same in dif- 
ferent places in the universe. 



Quasar absorption lines, generated by the absorption of light along 
the line of sight to quasars by intervening gas clouds, provide a 
method through which variation in certain fundamental dimension- 
less constants can be constrained. For high enough absorption red- 
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In this section, we describe the many-multiplet method and 
discuss its previous application. In section |2] we discuss the origin 
and treatment of the spectral data used. We describe our analysis 
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methodology in section[3] We present our new VLT results in sec- 
tion |4] We combine the VLT results with previous Keck results in 
section [5] and explore different models for the observed values of 
Aa/a. We discuss potential systematic errors in section[6] 

1.1 The many-multiplet method (MM method) 

Fine structure splitting of atomic transitions is observed as a con- 
sequence of Q being non-zero. For alkali-doublet (AD) type transi- 
tions, the separation in the wavelength of transitions is proportional 
to o? . More generally, one can derive the observed change in the 
fine-structure constant, 

Aq/q = (of^ - qo)/qo (1) 

where az is the value measured at absorption redshift z and ao 
is the present value of a. For a given transition, the observed 
wavenumber, Wz , depends on a as 

Uz = + qx (2) 

where 

provided that | Aa/a <C 1|. g is the "sensitivity coefficient" which 
determines how sensitive a particular transition is to a change in a 
( |Dzubaetal.[r 999b, Murp hy et al.|2003[ l. 

Comparing the wavelengths of AD absorption transitions to 
laboratory wavelength values for those transitions yields the AD 
method. The AD method does not, however, make full use of the 
available information in the spectrum, as it only compares transi- 
tions which arise from the same ground state of the same atomic 
species. By examining transitions from the ground state of differ- 
ent atomic species, one gains a statistical improvement (from the 
use of more transitions) and a sensitivity improvement (as transi- 
tions from different atomic species may display widely differing 
sensitivities to a change in a). Comparing many transitions from 
different multiplets, possibly from different atoms/ions, yields the 
MM method ( |Dzuba et al.|1999b|a[[Webb et al.|1999) , which typi- 
cally yields an order of magnitude improvement in sensitivity over 
the AD method, as well as significant resistance to certain system- 
atics. 

The transitions used by the MM method have q coefficients 
that are approximately zero ("anchor transitions"), positive ("posi- 
tive shifters") or negative ("negative shifters"). It is not the absolute 
values of the q coefficients which are important, but instead the dif- 
ferences in q between the transitions used; transitions with greater 
differences in q provide increased sensitivity to detect a change in 
a. The values of q must be calculated using quantum many-body 
methods ( |Dzuba et al..l999b.ai ). 

1.2 Previous results 

In the first application of the MM method, Webb et aT] \\999) 
examined 30 absorbers using spectra from HIRES (High Resolu- 
tion Echelle Spectrometer) on the Keck Telescope, in Hawaii. This 
yielded evidence that a may have been smaller in the absorption 
systems (compared to the current laboratory value) at the 3a level 
using Fe II and Mg I/II transitions. The addition of more absorbers 
by 2001 jWebb et al.|[200T| [Murphy et aT]|2001a| l increased the 
significance of the detection to 4c7. In addition to the use of Fe II 
and Mg I/II transitions, these works made use of the Ni Il/Cr Il/Zn II 



transitions. The two different sets of transitions yielded consistent 
evidence for a change in a. The sample size was increased to 128 
absorbers by |Murphy et al. 1 12003), who found > 5a evidence for a 
change in a. By 2004, this was increased to 143 absorption systems 
with further data from the Keck telescope ( [Murphy et al.|[2004l l. 
[Murphy et aT] ( |2004) give a weighted mean over all the systems of 
Aq/o? = (—0.57 ± 0.11) X 10~^ - evidence that the a may have 
been smaller at high redshift at the 5a level. However, as the Keck 
telescope is situated in the northern hemisphere (at ~ 20°N), it is 
unable to observe much of the southern sky, and therefore Aa/a is 
preferentially measured in one celestial hemisphere. Additionally, 
the use of only a single telescope and instrument is undesirable; 
one can postulate that the measured results are due to some instru- 
mental systematic. 

Although it is possible to conjecture that the Keck results are 
due to some instrumental systematic, we note that there have been 
thorough analyses of many potential systematic effects (Murphy] 
[et al. 2001b, 2003). At present, there is no evidence or analysis in 
the literature which identifies any specific systematic of sufficient 
magnitude and character to explain the Keck results. 

The only other statistical sample to search for a change in a 
using the MM method is the work of [Chand et al.[(2004 ). They 
analysed 23 absorption systems using spectra from UVES (the Ul- 
traviolet and Visual Echelle Spectrograph) on the VLT (Very Large 
Telescope), in Chile. Although they reported a statistical error on 
Aa/a of just 0.06 x 10"'', it has been shown that the errors were 
under-estimated by a factor of 6 and that the optimisation algorithm 
was demonstrably unreliable, failing to reach a solution in most of 
the absorbers (Murph y"iraLl[2007l [2008l l. [Srianand et al.[(2007l l 
reply to Murphy et al.l( [2007^ ^ and state that an updated analysis of 
the |(Chand et al.| sample shows adequate performance of the optimi- 
sation algorithm, though no updated values of Aa/a for individual 
absorbers are provided and so we are unable to comment further on 
the robustness or otherwise of that work. 

1.3 Objective & overview 

In section [2] we describe the data that form the basis of our anal- 
ysis. In section [3] we describe our analysis methodology in detail, 
with a particular focus on improvements made to our application of 
the MM method since the results of [Murphy et al.[ ( (2004l >. In that 
section, we also discuss certain aspects of the MM method which 
we believe should be re-emphasised. We further introduce a modi- 
fication of the Least Trimmed Squares method - a robust statistical 
method which assists with the detection of outliers in a dataset. In 
section |4] we present the results of our MM analysis, giving esti- 
mates of Aa/a for our sample. We also give weighted mean and 
dipole fits to our VLT data. In section |5] we combine our results 
with those of Murphy et al.[ ( [2004^ to derive the main results of this 
work. In section|6] we discuss potential systematic errors. 

This paper analyses the same data set presented in [Webb et al.[ 
( [20 10^ , however we provide substantially more detail here. 



2 QUASAR SPECTRA AND LABORATORY 
WAVELENGTH TRANSITIONS 

We have used publically available spectra from VLTAJVES, in 
Chile (at 25° S), to generate a MM sample similar in size to the 
Keck sample. To do this, we have searched the ESO Science 
Archive for VLTAJVES observations which might be used to con- 
strain Aa/a. Priority was given to spectra with high expected 
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signal-to-noise, although obtaining a sufficient number of high red- 
shift objects and maximising sky coverage was also a consider- 
ation. The ESO MIDAS pipeline was used for the first stage of 
the data reduction, including wavelength calibration, although en- 
hancements were made to derive a more robust and accurate wave- 
length solution from an improved selection of thorium-argon cal- 
ibration lamp transitions (Mu rphy et al.|2007] >. In the next stage, 
echelle spectral orders from several exposures for each quasar were 
combined using the program UVES_POPLEfQ which was specifi- 
cally designed for this purpose. A separate paper will be presented 
which describes the data reduction and processing in considerable 
detail (Murphy et. al., in prep). 

The signal-to-noise ratio (SNR) of the combined data is gen- 
erally significantly better than the Keck data used by |Murphy et al.| 
1 2004 1. The SNR per 2.5 km s^^ pixel in the VLT sample ranges 
from ~ 8 up to ~ 190 (see figure [9](. The higher SNR for the 
VLT sample allows better determination of the velocity structure 
required for a particular absorber, and reduces the statistical error 
on derived values of Aa/a. 

UVES has three CCD chips (1 in the blue arm and 2 in 
the red arm), which allows almost complete wavelength coverage 
(~ 3, 000 A to ~ 10, 000 A) to be obtained with a single setting for 
the echelle grating. A gap occurs at either ~ 4, 500A or ~ 5, 500A 
depending on which dichroic is used. The Keck spectra in the sam- 
ple of [Murphy et al] p004l > were taken when HIRES had only 
1 chip, thereby necessitating multiple exposures to achieve good 
wavelength coverage. In many of the Keck spectra, wavelength 
coverage gaps are visible, whereas in the VLT spectra there are al- 
most no spectral gaps. In principle, this allows the use of a greater 
number of absorbers per sightline. Similarly, this in principle al- 
lows the use of more transitions per absorber 



2.1 Extraction problems 

The MIDAS pipeline appears to incorrectly estimate the errors as- 
sociated with the flux data points in the base of saturated lines. In 
particular, the dispersion of the flux data points is too large to be 
accounted for by the statistical error. Fitting a straight line through 
the base of saturated lines typically produces ^ 2, where is 
measured with respect to the weighted mean flux and divided 
by the number of degrees of freedom, v. The problem is somewhat 
more noticeable in the blue end of the spectra. Although it is dif- 
ficult to determine precisely what happens in regions of low, but 
non-zero flux, we believe that the errors there are also underesti- 
mated. The effect of this is to give falsely high precision on any 
quantity derived from these data points (including Aa/a). Addi- 
tionally, one cannot fit plausible models to data involving regions 
of low or negligible flux. We demonstrate this effect in figure[T| and 
give more details there. 

When we combine individual exposures into a final spectrum 
using UVES_POPLER, UVES_POPLER provides a check on the con- 
sistency of the exposures contributing to each pixel. As the com- 
bined value of each flux pixel is given as a weighted mean of 
pixels from the contributing spectra, this consistency check is a 
value of for each pixel about the weighted mean. We attempt 
to correct for the problem described earlier, and any other prob- 
lems which cause inconsistency between contributing spectra, by 
applying the following algorithm. For each pixel in the combined 



Murphy, M. T., |http : / /astronomy . swin . edu ■ au/ -mmurphy/| 
|UVES_popler ■ html] 



O 




3500 



4000 



4500 



5000 



Wavelength (A) 



Figure 1. Example of the under-estimation of flux uncertainties in the base 
of saturated lines for the spectrum of J0407 18— 441013. The quantity R 
plotted on the vertical axis is the ratio of the RMS of the normalised flux 
array to the average value of the RMS array in selected regions of appar- 
ently zero flux. The RMS an'ay is a modified version of the error array 
produced by tJVES_POPLER which attempts to account for inter-pixel cor- 
relations. The quantity R should be approximately 1 if the RMS array is 
a good representation of the statistical uncertainty, but will be significantly 
greater than 1 if the eiTors are under-estimated. One can easily see from this 
plot that all the values are greater than 1. Here, the mean value is R = 2.08, 
with a standard error of 0.10, demonstrating inconsistency with i? = 1 at 
the > lOo- level. 



spectrum, we take a region of five pixels centred on that point. We 
then take the median of the x^ values just described that are asso- 
ciated with those five points. We then multiply the error estimate 
for that spectral pixel by the square root of that median value (that 
is, CTi — )• CTi X Y/iiiedian[x?]). 

However, we only increase the error array if median(x^) > 
1. This is because the spectral errors are generated from photon- 
counting. Poisson statistics requires that the error from photon- 
counting is a lower limit on the true error. Thus, it is difficult to 
justify decreasing the error array without particular evidence that 
the spectral uncertainty is systematically under-estimated. 



2.2 Other data and fitting problems 

Inspection of the spectra which contribute to our exposures reveal 
problems which could prevent acceptable fits being achieved, and 
which might interfere with the determination of Aq / a. We discuss 
these here. 

Cosmic rays. When the raw spectra are combined, 
UVES_POPLER attempts to automatically remove cosmic rays 
through cr-clipping. That is, pixels from contributing spectra 
which are more than some number of standard deviations from 
the weighted mean (calculated excluding the pixel that is being 
tested) are discarded from the combination. The default rejection 
threshold is 3(t, although this is user-adjustable. A second round of 
spectral cleaning is done manually. More details will be provided 
in Murphy et al. (in preparation). 

Sky absorption lines. Sky absorption lines are seen with vary- 
ing intensities depending on the orientation of the telescope. Tran- 
sitions from quasar absorbers sometimes fall in regions affected 
by sky absorption. In principle, atmospheric absorption could be 
fitted simultaneously with the absorption caused by the quasar ab- 
sorber, allowing the use of these transitions. However, our spec- 
tra have been converted to heliocentric wavelengths. As a result, 
the atmospheric transitions are shifted by velocities of up to ± ~ 
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30 kms~^. The co-addition of exposures taken at significantly dif- 
ferent times of the year therefore creates a complex absorption pat- 
tern. Rather than attempt to model this, we simply do not use tran- 
sitions which fall in regions of spectra which appear to be affected 
by sky absorption. 

Sky emission lines. The MIDAS program should automatically 
remove sky emission as part of the pipeline. However, we have no- 
ticed that strong sky emission lines are not necessarily properly 
subtracted. This can be seen visually in the spectrum, with sharp 
dips and spikes at the location of strong sky emission lines. Tran- 
sitions from the quasar absorbers are sometimes affected by this 
problem. Where we suspect that the spectral data are affected by 
sky emission lines, we either clip out the affected pixels (if they 
cover a small region of spectrum we want to use) or do not use the 
affected transition. 

Gravitational lenses. We discovered a small number of ab- 
sorbers for which the velocity structure appeared to be the same for 
transtions arising from the same ground state, but where the line 
intensities differed substantially. A particular example of this is the 
z = 0.82 absorber along the line of sight to J081331+254503. Fur- 
ther investigation demonstrated that this quasar is a known grav- 
itational lens. The complex line-of-sight geometry causes the ef- 
fect described. We discarded any system for which this problem 
appeared, and where we could identify the quasar as a known grav- 
itationally lensed system. 
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Figure 2. Our MM fit to the absorber at 2abs = 1.018 toward 
J220852— 1934359, which is apparently well fitted by a single component. 
The horizontal scale indicates the velocity difference from the arbitrary red- 
shift stated at the bottom for the given data points. The black line indicates 
the observed normalised flux, with the green line indicating our best fit solu- 
tion. At the top of each box, the black line indicates the standardised resid- 
uals (that is, [data - modelj/en'or), with the red lines indicating itlcr. The 
position of the blue tick marks indicates the fitted position of the single 
component. 



2.3 Atomic data 

In table |Bl| we present the atomic data and q coefficients used for 
our analysis. Our g-coefficients are from |Dzuba et al.^ ( ,1999b|aj> and 
[Berengut et al.|p004a|b^ and we provide all references for the other 
atomic data in the table caption. 



3 ANALYSIS METHODOLOGY 
3.1 Voigt profile fitting 

For each absorption system, we attempt to model all available MM 
transitions simultaneously (see appendix [B] for a list of these tran- 
sitions). We have used the non-linear least squares Voigt profile fit- 
ting program VPFIlj^ which was used to produce the Keck Aa/a 
results. Aq/q is included explicitly as a free parameter in the fit. 
Statistical uncertainties are determined from the diagonal terms of 
the covariance matrix at the best-fittin g solution, but are multiplied 



by y Xv to yield confidence intervals ( Press et al 



19921. We have 



verified through Markov Chain Monte Carlo simulations that the 
statistical uncertainties produced by VPFIT for Aa/a, crAa/a^ are 
robust ( [King et al.|2009[ l. We emphasise that the statistical infer- 
ence techniques used to measure Aa/a from each spectrum are 
standard and well-established. 



3.2 Modelling the velocity structure 

Some absorption complexes are well-modelled by a single Voigt 
profile, such as the one shown in figure|2] However, it is clear that 
almost all absorption profiles are not well modelled by a single 
Voigt profile. However, it is possible to obtain statistically accept- 
able fits by using multiple Voigt profiles ("velocity components"). 
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Unfortunately, there is no way of knowing a priori how many com- 
ponents must be used to obtain a statistically acceptable fit. Thus, 
the process of modelling the absorption profile amounts to adding 
components until a physically realistic, statistically acceptable fit is 
reached. 

Although we tie corresponding velocity components for 
species with similar ionisation potentials in our modelling proce- 
dure, this does not equate to the assumption that our results rely 
inherently on an assumption of identical velocity structure amongst 
all species fitted. In fact the assumptions we make in deriving the 
results presented in this paper are: (a) that the velocity structures 
for species with similar ionisation potentials are sufficiently similar 
such that it is statistically valid to model multiple species simul- 
taneously with one velocity structure (that is, any kinematic seg- 
regation which may be present is small relative to the observed 
linewidths), and (b) that any kinematic segregation between dif- 
ferent species is random in nature between different absorbers i.e. 
the sign of any velocity shift which may be present as a result of 
kinematic segregation is random between any pair of species in dif- 
ferent absorbers). Assumption (b) is justified because we are view- 
ing absoiption systems, where the quasar light passes through the 
absorption cloud along the line of sight, rather than optically thick 
emission systems, where we might receive light preferentially from 
specific components of the emitting region. If assumption (b) were 
found to be violated this would mean that kinematic segregation in 
different absorbers would be spatially correlated over cosmological 
scales, implying that the atoms in different absorbers must "know" 
about the presence and location of the atoms in the other absorbers, 
and that there is coherent organisation of the different atoms over 
cosmological distances. We do not consider such an arrangement 
to be plausible. 

We have three criteria for a statistically acceptable fit. 

i) Normalised of order unity. As a basic statistical test, we 
require that ~ 1- This follows directly from the fact that the 
X^ distribution with v degrees of freedom has mean v. However, 
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this criterion is not tlie only criterion that must be used. Adding 
components until 1 only suggests that the dispersion of the 

spectral data about the model is what one would expect for a good 
model. 

[Murphy et al.| ( (2008) l demonstrate through simulations that, at 
least for one synthetic spectrum considered, "underfitting" of spec- 
tra may lead to significant differences between the input and output 
Aa/a values, whereas "overfitting" seems not to induce a signifi- 
cant difference. We are therefore particularly cautious about under- 
fitting spectra. 

A model with x2 J; 1 ni''y be a reasonable model, but this 
test cannot determine whether it is the optimal model. Under the 
maximum likelihood method, if one can find a model which better 
explains the data, then this model should be preferred. The addi- 
tion of components will always improve x^, and therefore one con- 
siders whether the improvement is significant. A rigorous way to 
proceed is to to perform a statistical signifiance test on every com- 
ponent added, and only accept components which are statistically 
significant. However, this is laborious. Instead, other heuristics are 
available which tend to lead to reasonable choices of models. One 
of these is to try to find the model which minimises x^- If one adds 
a component and xi increases, this suggests that the extra compo- 
nent is not supported by the data. 

Other methods are available, which penalise free parameters 
more or less strongly. We have chos en to use the Akaike Informa- 
tion Criterion (AIC) ( |Akaike|l974| , defined as AIC = + 2p, 
where p here is the number of free parameters. In fact, the AIC is 
only correct in the limit of large n/p (where n is the number of 
spectral points included in the fit). This is not generally true for our 
fits. Thus, we use the AIC corrected for finite sample sizes, defined 
as 



AICC = X + 2p + 



2p(p + 1) 



(4) 



l |Sugiura|1978^ . A significant advantage of the AICC is that it al- 
lows multiple models to be compared simultaneously. If several 
competing models are being considered, one chooses the model 
which has the lowest AICC. The actual value of the AICC is not 
important; only relative differences matter. The suggested inter- 
pretation scale for the AICC is the Jeffreys' scale (Jeffreys|196I[ l, 
where AAICC = 5 and AAICC = 10 are considered strong 
and very strong evidence against the weaker model respectively. 
AAICC < 5 means that the difference between the models is 
"barely worth mentioning". These considerations lead to our sec- 
ond and third criteria: 

ii) Choose the model which has the lowest AICC. In this way, 
we attempt to find a model which best explains the data. 

Hi) No substantial correlations in the residuals. For a well- 
fitting model, the residuals of the fit (defined as Vi = [data — 
model] /error) should be random. This means that no significant 
correlations should occur in the residuals. Such correlations sug- 
gest a deficient model. Although statistically improbable runs of 
residuals can be detected through formal tests (Wald & Wolfowitz| 
|I940[ >, in practice this is not necessary. The addition of components 
in regions where such correlations are visible generally removes the 
effect, and decreases the AICC as well. We accept the addition of 
components which decreases the AICC. 

Thus, in fitting the observed absorption profiles, we attempt to 
obtain a fit which has x^ ~ 1 ''"^ the minimum AICC possible 
and has no substantial correlations in the residuals. We regard this 
fit as our final fit. We show in figure[3]an example of a complicated 
absorption system, which requires many components to achieve a 



statistically acceptable fit. We do not include Aa/a as a parameter 
while establishing the final fit. That is, we avoid biasing Aa/a by 
only adding it once the final model is found. 

Determining the best-fitting number of velocity components, 
i.e. minimising the AICC, can be difficult. Considerable effort was 
expended to find the lowest possible AICC; that is, for each ab- 
sorber we could not find a model with a lower AICC. Thus, we can 
say that our models should be close to the best model for the data. 
Any inaccuracy in the measurement of Aa/a which is induced by 
our model being sub-optimal will have a random magnitude and 
sign, and therefore average out when considering a large number of 
Aa/a values. 

3.3 Random and systematic errors 

Unfortunately, the Voigt profile decomposition of any given absorp- 
tion system, with any particular choice of fitted transitions, is not 
unique. Errors in modelling the velocity structure may impact the 
fitted value of Aa/a. We thus distinguish between three different 
types of errors which may affect Aa/a: 

Statistical errors. These errors are simply the errors on Aa/a 
which derive from the propagation of uncertainty from the flux er- 
ror array via the Voigt profile model. These errors are simply the 
errors produced by VPFIT from the covariance matrix at the best- 
fitting solution. 

Random errors. Random errors are any effects which might 
cause Aa/a to be measured inaccurately when considering a sin- 
gle absorber. Significant errors made in determining the correct ve- 
locity structure could cause an error of this type. [Murphy et al.| 
(f2008 1 demonstrated that an under-fitted spectrum (i.e one with a 
deficient model for the velocity structure) gives erroneous values 
of Aa/a. Other potential causes of random errors include: i) kine- 
matic segregation of different elements (which cannot be preferen- 
tially biased along the radial sightline over a large number of ab- 
sorbers); ii) random blends with other transitions arising in differ- 
ent absorbers (at different redshifts); Hi) random departures of the 
wavelength calibration solution from the true wavelength scale; iv) 
cosmic rays and other uncleaned data glitches; v) incorrect determi- 
nation of the broadening mechanism (turbulent or thermal) for any 
component. Importantly, the effect of random errors will average 
to zero when considering an ensemble of velocity components in a 
given absorber, and also an ensemble of absorbers. This is because 
these errors will displace Aa/a to be more positive as often as they 
will displace it to be more negative. When considering only a single 
absorber, this type of effect must be considered a systematic error. 
However, when considering an ensemble of absorbers, this effect is 
a random effect, which adds extra scatter into the data. Any small 
kinematic shifts between different species in corresponding veloc- 
ity components are of random sign and magnitude, and therefore 
any effect this has on Aa/a will tend to average out even when 
considering many components within one absorber. When consid- 
ering many absorbers, each with multiple velocity components, any 
bias introduced into Aa/a must average to zero. This point was 
also discussed by [Murphy et al.|pOOTbl l and [Murphy et al.[ ( |2003l l. 
We explicitly constrain potential kinematic segregation of different 
species for a subsample of absorbers used in our analysis in Section 

Systematic errors. This is any error which systematically af- 
fects the values of Aa/a measured over an ensemble of absorbers. 
Such effects would include: i) inaccuracies in the laboratory wave- 
lengths for particular MM transitions; ii) a different heavy isotope 
abundance for Mg in the clouds relative to terrestrial values; Hi) 
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Figure 3. Part of our MM fit to the absorber at Zabs = 2.152 toward J233446— 090812. This is a complex absorption system, requiring many components 
in order to achieve a statistically acceptable (xj ^ 1) fit. See figure|2|for a description of the various lines, marks and labels. We draw the reader's attention 
to the presence of a wide range of transitions, some with relatively small magnitude q coefficients (Nin A 1709 and the two Mgll transitions), some with 
large magnitude, positive q coefficients (Fe II A2382, 2344, 2374, 2260, the Zn II transitions and the Mn II A2576 transition) and some with large magnitude, 
negative q coefficients (Fe II A1608, Ni II A1741, 1751 and the Cm transitions). The Al III transitions fitted are not shown, but these have a minimal impact 
on the value of and precision for Ao/a. Note that for the stronger species, such as the Fe II A2382 and Mg II transitions, the centre regions of the profile are 
saturated, and thus a constraint on Ao/a only comes from the optically thin wings. Conversely, for the weaker species (for example, Fe II A1608, 2260 and 
the Ni II transitions) most or all of the profile is optically thin, and thus a constraint on Aa /a is derived across the whole profile. Importantly, a single velocity 
structure model provides a good model to all the observed MM transitions. This serves to validate an assumption underlying the MM method, namely that 
kinematic segregation of the different species - if present - must be relatively small. We note the presence of two weak interlopers in Fe II A2260 and one in 
Fe ll A1608, yielding an additional three tick marks. The regions of the Mg II transitions near t) fa have one less component than the corresponding regions 
in Fe II . This is because this region is saturated, and VPFIT would not support so many strong components in the saturated region. This does not significantly 
affect Aa/a because changes in parameters for lines in the middle of saturated regions have a marginal impact on and therefore on Aa/a. We have 
clipped pixels out of the spectrum for the region containing Cr II A2062 and Zn II A2062 because of a problem affecting the spectrum in that region. 
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systematic blends with nearby but unmodelled transitions in the 
same absorber; iv) time-invariant differential light paths through 
the telescope for different wavelengths; v) atmospheric dispersion 
for spectra taken without an image rotator; vi) differential isotopic 
saturation (see [Murphy et al.|2003| l, and; vii) wavelength miscali- 
bration due to thorium-argon line list inaccuracies. [Murphy et al!] 
( [2003^ consider many possible systematic effects in detail. Note 
that these effects will not necessarily produce the same spurious 
shift in Aa/a in every absorber. For example, if some transitions 
have inaccurate laboratory wavelengths, the effect on Aa/a will 
depend on which other transitions are fitted in the same absorber. 
Nevertheless, the above effects are considered systematic errors be- 
cause they will cause similar or correlated shifts in certain subsets 
of absorbers. 

There is a crucial distinction between these effects: some ef- 
fects may be considered systematics in single absorbers, but are not 
systematics in an ensemble of absorbers. For this reason, we are 
cautious against placing too much emphasis on the interpretation 
of the value of Aa/a from any individual absorber. We demon- 
strate later that the impact of random effects is non-negligible. We 
explain our treatment of random effects in section [J.5.3| 

3.3.1 Physical constraints 

As in the previous Keck analyses, we require that the velocity width 
(6-parameter, with b — \/2c") of each modelled component for a 
particular species in the fit is related to the corresponding compo- 
nents for other species. The two extreme cases are wholly thermal 
broadening (6^ = 2kT/M for a species with atomic mass A/) and 
wholly turbulent broadening (b^ = b^urh)- In general, there will be 
contributions from the two mechanisms (6^ = 2kT/M + btu^b)' 
however we have found that most systems are generally well-fitted 
with turbulent broadening. As noted by |Murphy et al.| ( |2003| l, it is 
possible to explicitly determine the degree of thermal and turbulent 
broadening, however in this circumstance the fe-parameters are gen- 
erally poorly determined, which makes the optimisation difficult. 

It turns out that the turbulent fit is preferred on the basis of the 
AICC in 7 1 percent of the absorbers, and the thermal fit in 29 per- 
cent of the absorbers. However, it should be noted that the fits were 
initially constructed with turbulent broadening and then converted 
to thermal broadening. It may be that if the fits were constructed 
thermally and then converted to turbulent that these figures might 
change significantly. We emphasise that mistakes made in choos- 
ing turbulent or thermal fitting may bias Aa/a for a single ab- 
sorber, but these effects must average to zero over a large number 
of absorbers due to the random nature of the bias from absorber to 
absorber. 

The previous analyses of the Keck results required that Aa/a 
calculated using both thermal and turbulent fits differed by no more 
than la for that absorber to be included in their final analysis, 
where the difference is considered only in terms of the statistical 
error. However, the generally higher SNR of the VLT spectral data 
(compared to the Keck spectral data) often leads to very precise 
statistical bounds on Aa/a. This makes the Icr-difference crite- 
rion difficult to fulfil in a significant number of cases. We describe 
below how we resolve any potential inconsistency between Aa/a 
values from the thermal and turbulent fits. 

In determining how to resolve any potential inconsistency, 
there are three cases to consider, i) Where the difference between 
the fits is substantial, as measured by the AICC, one wants to take 
the statistically preferred fit. ii) Where the quality of the fits is sim- 
ilar (AlCCturbuicnt ~ AlCCthcimai), and the values of Aa/a are 



the same, then it does not matter which fit is used. Hi) If the val- 
ues of the AICC for the thermal and turbulent fits are similar, but 
the values of Aa/a produced by those fits differ significantly, then 
the statistical precision accorded to Aa/a should be reduced to 
account for the conflicting evidence, and value of Aa/a should be 
somewhere between the two cases. 

To resolve this problem, and to be more accurate for our 
meaning in "similar" and "significant" above, we use a method- 
of-moments estimator which takes into account the relative differ- 
ences in the AICC and the agreement, or otherwise, of the values 
of Aa/a. We estimate the underlying probability distribution of 
Aa/a for the absorber in question as the weighted sum of two 
Gaussian distributions (one for the thermal result, one for the tur- 
bulent), with centroids given by the best fit value of Aa/a for each 
fit, and a equal to cr/^^/^ for each fit. We weight the sum by the pe- 
nalised likelihood of the fits, via the AICC (see |Liddle|2007l l. That 
is, if 

k = exp(- AlCCturbuicnt /2) + exp(-AICCthermal/2), (5) 



jl = exp(-AICCturbulont/2)/fc, (6) 

j2 = exp(-AICCthcrmal/2)/fc, (7) 

ai — Aa/aturbulent, (8) 

a2 = Aa/athermal, (9) 

51 = cr(Aa/aturbuiont), and (10) 

52 = Cr(Aa/athcrmal) (H) 

then matching the first two moments of our weighted sum of distri- 
butions with a Gaussian yields 

m — Aa/a = jiai + j2a2, and (12) 
o-Aq/q = V jiSi + J2si +jiaj + j2al - m?. (13) 



This covers all the cases described above. In particular, where the 
AICC is similar but Aa/a differs significantly between the tur- 
bulent and thermal fits, the estimated error increases with the dif- 
ference between them, providing resistance to incorrectly deter- 
mining the line broadening mechanism. To see this, note that with 
ji + j2 ~ 1, equation[T3]reduces to 

<yc^c/c = + (1 - ji)si + ji(l - - a2Y- (14) 

Thus, errors only ever increase from our smallest error estimate, 
and therefore this method could be considered conservative. In the 
event where one broadening mechanism is significantly preferred, 
then our result will be effectively the same as if only that broaden- 
ing mechanism was considered. For the case where the fits are sta- 
tistically indistinguishable (ji = J2), Aa/a is given by the simple 
mean of the two values of Aa/a, and the variance is the simple 
mean of the individual variances plus 0.25(ai — a2Y . 

3.3.2 Aim 

In principle the Al III transitions can be included in a MM fit, how- 
ever its ionisation potential is somewhat different to the other MM 
transitions described. Therefore, the Al III transitions may not arise 
from the same location, and therefore velocity, as the other MM 
transitions. If the Aim transitions arise at significantly different 
velocities to the other MM transitions then an error would be intro- 
duced into Aa/a for a system with Aim included (although this 
effect must average to zero over a large number of absorbers, as 
there is no reason for a systematic bias in the centroid of the Al III 
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transitions with respect to the other MM transitions along a line of 
sight to Earth). 

Generally, the profiles for different transitions for the other 
MM transitions used correlate well with each other. By this, we 
mean that the relative column densities between corresponding ve- 
locity components are similar for different MM transitions. How- 
ever, we have noticed that the absorption profiles for some Al III 
transitions in some absorbers differ significantly in the relative line 
strengths between components, when compared to other MM tran- 
sitions. Importantly, we found some absorbers where it was diffi- 
cult to apply the same velocity structure model to Al III transitions 
and the other MM transitions simultaneously. For this reason, we 
are therefore cautious in fitting Al III together with the other MM 
transitions. 

Therefore, we include and model Al III if both the transitions 
are available, and allow the spectral data to contribute to Aq/q: de- 
rived from the other MM transitions for that absorber, but do not 
constrain the modelled structure with the velocity structure from 
other MM transitions. Given the small difference in the q coeffi- 
cients between the two Al III transitions (Ag ~ 250) the statistical 
contribution of Al III to Aa/a is low, however given that the ex- 
posures have already been obtained it is prudent to try to maximise 
our use of the existing data. 

As an example: if Al III is not utilised in the z — 1.857 ab- 
sorber towards J013105— 213446, the turbulent fit value of Aa/a 
changes from (0.30 ± 1.43) x 10"^ to (0.43 ± 1.44) x lO"'-'. 
Similarly, if Al III is not utilised in the z = 1.71 absorber towards 
J014333-391700, the turbulent fit value of Aa/a changes from 
(-2.20 ± 2.67) X 10"^ to (-2.32 ± 2.68) x 10"^ 

We note that previous works have included the Al III tran- 
sitions as part of the MM analysis, and it was demonstrated by 
[Murphy et al.| ( |2003) that the inclusion of these transitions did not 
significantly alter the Keck results. Nevertheless, the approach we 
have adopted is conservative. 

We have also observed less substantial relative line strength 
differences between the Mg I transitions and other MM transitions, 
but in no case did we find a system where we could not apply the 
same velocity structure model to the Mg I and MM transitions, and 
so we include the Mg I transitions in the full MM analysis. 
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Figure 4. The Si II transitions from our MM fit to the absorber at z^bs = 
2.638 toward J212912— 153848, shown without an included interloper. The 
horizontal scale indicates the velocity difference from the arbitrary redshift 
stated at the bottom for the given data points. The black line indicates the 
observed normalised flux, with the green line indicating our best fit solu- 
tion. At the top of each box, the black line indicates the standardised resid- 
uals (that is, [data - modelj/en'or), with the red lines indicating Her. The 
position of the blue tick marks indicates the fitted position of the single 
component. The strength of the Sill A1808 transition can be predicted 
from the model for Sill A1526. Sill A1808 shows excess absorption at 
V ft! — 40kms~^ which cannot be explained by Sill A1526. To account 
for this, we include a single, unconstrained interloper. After including the 
interloper, the fit is statistically acceptable. 



interloping transitions. This means that the interloping transition 
can be modelled if it is an MM transition from an absorber at a dif- 
ferent redshift, or if it is Si IV AA1393 or 1402. If the interloper is 
from the C IV doublet, the contamination can also be modelled de- 
spite the fact that the rest wavelengths for this doublet are relatively 
poorly known. This is done by allowing the C IV transitions to have 
a separate value of Aa/a, which is then discarded. This extra pa- 
rameter effectively absorbs any error introduced through inaccurate 
knowledge of the rest wavelengths. 



3.4 Interlopers 

Some transitions display excess absorption beyond what is pre- 
dicted from other transitions from the same atomic species. This ex- 
cess absorption is caused by another absorber located along the line 
of sight to the quasar, which is usually extragalactic. Even where 
a prediction cannot be made from another transition of the same 
atomic species, interlopers can still be detected when many transi- 
tions are fitted together, as the redshifts and fe-parameters of each 
component are constrained. Although the contaminated sections of 
spectrum can be discarded, it is often possible to adequately model 
the contamination, thereby maximising use of the spectral data. 

We distinguish two types of interlopers: identifiable and 
tmidentifiable. 



3.4.1 Identifiable interlopers 

In some cases, the interlopers can be identified and modelled si- 
multaneously with the MM transitions. By "identified" we mean 
that the redshift, atomic species and wavelength of the transition 
which causes the excess absorption can be determined. In princi- 
ple, one needs accurate rest wavelengths and g-coefficients for the 



3.4.2 Unidentifiable interlopers 

In many cases, however, the interloping transition can not be identi- 
fied. In this case, our decision as to how to proceed depends on the 
degree of contamination. If the degree of contamination is small, 
and confined to a small area of the observed profile, we can in- 
clude unknown interloping transitions where the residuals of the 
fit ([data - model]/error) are bad until a statistically acceptable fit is 
achieved. Doing this provides a statistically acceptable model of the 
contamination. Note that the contribution to Aa/a of the affected 
MM transition will be reduced as a result of this, as the interlopers 
included are unconstrained by other spectral regions. We show an 
example of this in figure|4] 

As the degree of contamination begins to increase, the poten- 
tial error introduced into Aa/a may grow. Our treatment of tran- 
sitions affected by significant contamination depends on whether 
there are other transitions from the same species available which 
can be used to constrain the velocity structure for that transition. In 
the case of Fe II, a wide variety of transitions are often available. 
Effectively, a fit to several other transitions of the same species may 
allow the structure of the contamination to be determined, partic- 



ulaiiy if the SNR is high. However, in some cases, there may be 
no other transitions which can be used to obtain the velocity struc- 
ture. This may occur as a result of all of the transitions suffering 
from contamination from different absorbers, or because the spec- 
tra only includes one of the transitions of that species (due to gaps 
in the spectral coverage), or because the species has no other tran- 
sitions which could be used for that purpose. The last case is partic- 
ularly problematic for Al II, for which we only use the Al II A1670 
transition - there are no other Al II transitions which can be used 
to directly constrain the Al II structure. This issue also occurs for 
Si II - although in theory both the Si II A1526 and 1808 transitions 
can be used to constrain the velocity structure of Si II, the oscillator 
strength for the A1526 transition (/ ~ 0. 13) is much larger than for 
the A1808 transition (/ ~ 0.002). For many systems observed to 
have Si II A1526 absorption, the column density of Si II is not large 
enough to detect the A 1808 transition. Even if the A 1808 transition 
is detected, it may be too weak to provide a meaningful constraint 
on the Si II structure. 

In cases where we are unable to obtain a good constraint on 
the velocity structure of a particular species from other unaffected 
transitions of the same species, and the degree of contamination is 
not small, we are very cautious about inserting interlopers, due to 
the potential bias this could introduce into Aa/a for this absorber. 
Note that any bias introduced here is a random effect, and there- 
fore will average out when considering an ensemble of absorbers. 
Nevertheless, we wish to avoid introducing extra scatter into the 
Ao/q values where possible. In these cases, we clip out the pix- 
els which appear to be affected by contamination, leaving a wide 
buffer on either side. Note that we can only do this where another 
transition from the same species exists. Otherwise, components sit- 
uated in the middle of the clipped pixels might have very little, if 
any spectral data to constrain them, and thus their column densi- 
ties could take on values which would not be consistent with the 
general model used for the absorber. If there is no other transition 
for the species in this case, we simply do not use the transition. In 
figure |5] we show an example of where we have clipped out pixels 
because of contamination by an interloper. 

Fitting contamination has the potential to introduce a ran- 
dom bias for individual absorbers, however this may nevertheless 
reduce systematic effects. To see this consider a hypothetical ab- 
sorber where the Mgll AA2796, 2803, AllI A1670 and the Fell 
AAA2383, 2600, 2344 transitions are available, but that the AllI 
A 1670 transition suffers from some minor contamination in part of 
the observed profile, and the absorber of the contaminating transi- 
tion cannot be identified. One could simply ignore the Al II A1670 
transition and fit the Mgll and Fell transitions. Deriving Aa/a 
from just the Mg Il/Fe II combination is not very robust to a simple 
stretching or compression of the wavelength scale (Mur phy et al.| 
|2003( l. The Mg Il/Fe Il/Al II combination just described is much 
more robust against this effect, as the Mg II and Al II anchor tran- 
sitions are positioned on other side of the high-g Fe II transitions. 
Modelling the contamination of the Al II profile here may introduce 
a random error, but this can be averaged out in the context of many 
absorbers. The potential introduction of this random error is eas- 
ily justified by the increased resistance to systematic effects (such 
as wavelength scale distortions, which in principle could be com- 
mon to many absorbers or spectra). In the situation just described, 
the only choices are to discard the Al II transition, or to model the 
contamination. We choose the latter option where the degree of 
contamination is minor, and the former where the contamination 
is severe. 
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Figure 5. Part of our MM fit to the absorber at Zabs = 1-268 toward 
1005758—264314. The horizontal scale indicates the velocity difference 
from the arbitrary redshift stated at the bottom for the given data points. The 
black line indicates the observed normalised flux, with the green line indi- 
cating our best fit solution. At the top of each box, the black line indicates 
the standardised residuals (that is, [data - model]/error), with the red lines 
indicating itlcr. The position of the blue tick marks indicates the fitted po- 
sition of the single component. We show here the model from Mg II A2803 
plotted over the spectral region of Mg II A2796. One can see that there is 
significant excess absorption in Mg II A2796 for < —40 km s~^. There 
also appears to be excess absorption in the range — 40kms~'^ J; ^ J; 
—20 km ^ . Due to the wide range in velocity over which the absorption 
occurs, instead of attempting to model the absorption we clip away all pixels 
fori> < -20kms-i in the Mgll A2796 region for our actual fit. 

3.4.3 Transitions in the Lyman-a forest 

It also happens that some transitions fall in the Lyman-a forest 
("the forest"), a dense series of absorption lines blueward of the 
quasar Lyman-Q emission line. These transitions are caused by Hi 
absorption along the line of sight to the quasar. We are cautious 
about using MM transitions which fall in the forest, due to the un- 
certainties in determining the structure of the forest. Neverthless, 
the use of MM transitions in the forest may afford significantly 
better constraints on Aa/a. Although this often occurs with low-z 
Mg Il/Fe II absorbers, with some Fe II transitions falling in the for- 
est, it also occurs in high-z systems, where the Sill A1526/A1II 
A1670/FeII A1608 combination is common. Where the SNR is 
high for the transitions which fall in the forest, we model the forest 
structure with H I absorption. If the SNR ratio is low, determina- 
tion of the forest structure can be difficult, and therefore we do not 
utilise the contaminated transitions. Again, we emphasise that al- 
though this has the potential to introduce bias for a single absorber, 
because the contamination is random from absorber to absorber, it 
must average out over a large number of systems. We have used 
transitions which fall in the forest in 27 of the 1 142 spectral fitting 
regions in the VLT sample (2.4 percent). 

3.4.4 Cr II A2062, Zn II and Mg I A2026 

The Cr II A2062 (2062. 24A) and Zn II A2062 (2062. 66A) fines 
are relatively closely spaced, being separated by ~ 62kms~^. 
For narrow absorption systems, one can distinguish between these 
transitions as they do not overlap. However, these transitions are 
most commonly associatd with damped Lyman-a absorbers, where 
the velocity structure is generally complicated, and the system dis- 
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plays absorption over tens to hundreds of km/s. In this case, the Cr 
II A2062 and Zn II A2062 transitions often overlap. The velocity 
structure for these transitions can be determined by simultaneously 
modelling these transitions with the Cr II A2052, 2056, 2066 and 
Zn II A2062 transitions. 

There is one point of caution here. For high column density 
systems, a potential blend exists with Mg I A2026. Mg I A2026 is 
weak, with oscillator strength / = 0.113, and is rarely seen. One 
can in principle use the Mgl A2852 (/ = 1.83) information to 
constrain the Mg I structure. In this case, a joint fit of Zn II A2026 
and Mg I AA2026, 2852 will ensure that the Zn II A2026 results are 
not biased by any absorption due to Mg I A2026. However, the ab- 
sorbers for which Mg I A2026 might be detected are often at high 
redshift, in which case Mg I A2852 is often unusable, either due to 
heavy contamination by sky emission or absorption, or because it 
is located out of the red end of the spectral coverage. In this par- 
ticular circumstance, we are generally cautious about fitting Zn II 
A2026. Where we consider that the Zn II A2026 transition might be 
affected by Mg I A2026, and we are unable to utilise Mg I A2852, 
we do not include the Zn II A2026 transition. 



3.5 Aggregation of Aa/a values from many absorbers 

3.5. 1 Weighted mean 

If one assumes that all the Aq/q values are described by a constant 
offset of Aq/q from the laboratory values, one can combine the 
Aq/q values together using a weighted mean. This process is valid 
provided that the Aq/q values support a constant value of Aq/q. 
If Aq/q 7^ 0, then this implies that there must be a transition at 
some point from laboratory conditions (Aq/q = 0), and therefore 
the Aq/q values should be inspected to see if a transition point can 
be identified. Additionally, one must examine the residuals about 
the fit for a weighted mean, plotted against various parameters of 
interest (e.g. redshift and sky position) to determine if unmodelled 
trends exist. 



3.5.2 Dipolefit 

A dipole-l-monopole model constitutes the first two terms of the 
spherical harmonic expansion. The simplest dipole model is of the 
form 

Aq/q = Acos(e) + m, (15) 

where is the angle between the pole of the dipole and the sky 
position under consideration, A is an angular amplitude and m (the 
monopole) represents a possible offset of Aa/a from the labora- 
tory value. An equivalent (and more computationally convenient) 
form is 

AQ/Q = c.x + m, (16) 

where x is a unit vector pointing towards the direction under con- 
sideration and c contains the amplitude and direction information 
of the dipole. The components of c {cx, Cy and Cz) are easily related 
to the right ascension (RA) and declination (dec.) of the direction of 
the dipole, and |c| gives the magnitude of the dipole. In this form, 
Aq/q is linear in the Ci and so the Ci may be determined through 
weighted linear least squares. 

Although naively we might expect that m — 0, some theories 
contemplate otherwise. This could be possible if a depends on the 
local gravitational potential or density ( [Khoury & Weltman|2004| 



|Mota & Shaw|2007[ pTve & Pospelov|2008| l - laboratory condi- 
tions differ quite significantly in this regard to the conditions in the 
quasar absorbers. Note that by including the m term, one gets an 
explicit statistical test of this idea. Additionally, in the presence of 
temporal evolution of q, m amounts to its average value. In any 
particular redshift slice, m therefore represents the average angle- 
independent value of Aq/q. 

Equations [15] and [T6] take no account of any redshift depen- 
dence. Clearly, A = A{z). Nevertheless, a model which includes 
angular dependence is useful because it provides a method of de- 
tecting spatial variations in a which does not require the specifi- 
cation of a functional form for A{z). Use of this model is valid 
under several possible circumstances. One is where any variation 
in Q with redshift in the sample along a particular direction is small 
compared to variation in a in the opposite direction. This might be 
possible in our sample, depending on how a might vary, as we typ- 
ically probe lookback times of greater than 5 gigayears. Another 
is where a does vary significantly with redshift in our sample and 
the distribution of absorber redshifts does not vary greatly with sky 
position. With enough data, one could simply take redshift slices 
and apply this model to each redshift slice, thereby building up the 
functional form of A{z) in model-independent manner. However, 
given that we have only ~ 300 absorbers between the Keck and 
VLT samples, we simply cannot slice the data enough to do this 
for more than ~two redshift bins. Another consideration is the ef- 
fect of choosing a particular form of A{z). An incorrect choice of 
A{z) may reduce sensitivity to detect an effect, and could lead to 
the wrong conclusion if the choice is sufficiently bad. As a result, 
we explore an angle-dependence model initially, and later consider 
explicitly including distance dependence. 

We derive our uncertainties on the direction of dipole vectors 
by tranforming the covariance matrix from rectilinear coordinates 
{cx,Cy,Cz)lo spherical coordinates {r,(j},6), using the standard Ja- 
cobian matrix. That is, if J is the Jacobian matrix of the transforma- 
tion from rectilinear to spherical coordinates, and C is the covari- 
ance matrix calculated from the fit, then C = J • C ■ gives the 
approximate covariance matrix in spherical coordinates. The radial 
component, r, corresponds to the amplitude of the dipole. Our er- 
rors given on RA and dec. are thus linearised approximations based 
on the covariance matrix at the best-fitting solution, and should be 
regarded only as approximate. These error estimates will be inac- 
curate if they subtend a large fraction of the sky. 

Note that, by virtue of the fact that r ^ in spherical coordi- 
nates, the dipole amplitude. A, is not Gaussian. Thus, we perform 
a resampling bootstrap analysis jPress et al.|1992|f to derive an un- 
certainty for a dipole amplitude. Similarly, one cannot use a t-test 
to determine if A is significantly different from zero. Thus, we cal- 
culate the statistical signifiance, 1 — p, of the dipole model over 
the monopole model using a bootstrap method where we randomise 
Aq/q over sightlines, and from the observed distribution of de- 
termine the probability that a value of as good or better than that 
given by our observed dipole fit would occur by chance. One can 
also use analytic methods (C ooke & Lynden-Bell|2010j ), if desired. 
These methods should yield similar answers for large sample sizes. 
However, for small sample sizes, the results may differ somewhat 
(especially if the statistical uncertainties vary significantly in mag- 
nitude between the Aq/q values). As the dipole model will always 
improve the fit over a monopole, the statistical test is one-tailed, and 
so when we state the a-equivalence of a statistical significance, this 
is calculated as |probit(p/2)j in order to accord with the conven- 
tional usage, probit is the inverse normal cumulative distribution 
function. 
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Unless otherwise mentioned, we multiply uncertainty esti- 
mates on monopole values and sky coordinates by as a first- 
order correction for over- or under-dispersion about the fitted model 
( [Press et al.|T992l i. 

3.5.3 Robust estimate of random errors 

It may happen that the values of Aa/a demonstrate excess scat- 
ter about a particular model. Although this might be due to model 
misspecification, it is also due to the random effects described in 
section [33] Previous works (e.g. [Murphy et al.|2003| [2004^ have 
added a constant, (Trand, in quadrature with the uncertainty esti- 
mates of each Aa/a value, cTstat, until = 1 about the fitted 
model (i.e. a^^^ = errand + '^EtsLt)^ thereby making the data statis- 
tically consistent with each other under that model. This constant 
is an estimate of the aggregation of any errors which average to 
zero over a large number of systems (i.e. random errors), under the 
assumption that the magnitude of this effect is the same in each 
absorber. 

We are cautious, however, about not over-estimating the term 
f7"rand required on account of apparent outliers in the Aa/a values. 
We define an outlier as any point which has a standardised residual 
\ri\ — |(Aa/a — model prediction)/(Jtot| > 3, even after adding 
an appropriate (Jrand. The existence of outliers even after increas- 
ing the error bars implies that these systems may be affected by 
random or systematic errors which occur infrequently, or with sig- 
nificantly lower magnitude in most of the systems. Adding a con- 
stant in quadrature with the error estimates until xt ~ ^ will there- 
fore over-estimate the random error associated with most points, 
and under-estimate the random error associated with high scatter 
points, increasing the probability of a false negative (a failure to 
detect Aa/a ^ when it is in fact 7^ 0). 

To alleviate this problem, we use a modification of the Least 
Trimmed Squares (LTS) method ( Rousse euw| 1 984| l . Instead of fit- 
ting all n Aa/a values using weighted least squares, the LTS 
method traditionally fits only k — (n + p + l)/2 Aa/a values 
(where p is the number of parameters fitted) using weighted least 
squares, and searches for the combination of k data points and fit- 
ted model that yields the lowest weighted sum of squared residuals 
(X^). If one is more confident that the data contain only a small 
number of outliers, then one can choose k somewhere between 
(n+p+ 1)/2 and n. For n — >■ 00, the use ofk — {n+p+ 1) /2 will 
produce a very robust fit, where the fit is resistant to < 50 percent 
of the data being erroneous. However, for small n (e.g. n < 20) the 
resultant fit parameters can be significantly affected by small per- 
turbations to the fitted data points because of the strong sub-sample 
sensitivity of the LTS in small sample sizes. To help remedy this, 
we choose k = 0.85n when applying the LTS method. We utilise 
the Fast-LTS method jRousseeuw & Driessen|2006) l, an algorithm 
which determines the LTS solution fairly rapidly. 

The LTS method has traditionally been applied to fit models 
to data for which uncertainty estimates are not initially available, 
leading to the assignment of constant magnitude error estimates to 
all data points. However, our Aa/a values already have statisti- 
cal error estimates, which serve as lower bounds on the true errors. 
We propose a new variant of the LTS method to increase the error 
bars in quadrature with some (Jiand, but that addresses our con- 
cerns about false negatives. This method was inspired by [Pitsoulisl 
[& Zioutas]j2009^ , who use linear scaling of existing errors rather 
than our quadrature error addition. 

To obtain our robust estimate of (Trand, we increase the error 
bars until a robust scatter measure is what we would expect for 



a Gaussian distribution. We calculate our robust scatter measure, 
xi{k), as the sum of squared residuals per degree of freedom about 
a fit taken over only the k Aa/a values with the smallest squared 
residuals, rather than over all n Aa/a values. To estimate (Jrand! 
we slowly increase (Jrand from zero and add (Jrand in quadrature 
with the existing error bars. For each value of (Jrand, we calculate 
the LTS fit, then calculate xt{^)- We continue to increase (Jrand un- 
til xlW is equal to the expected value for a Gaussian distribution 
with large numbers of degrees of freedom, {xt (fc)) . If we define / 
as fc/n, then {xt{k)) is calculated as 



(xlik))^ J ^ dx (17) 



where a = probit[(l + ./)/2]. We take the value of (Jrand derived 
in this way as our estimate of the additional random error for the 
data given the model. For the standard least squares case (given by 
k — n), this yields (x^('^)) = 1> the well-known result that the 
expected value of xt is 1 for large numbers of degrees of freedom. 
Thus, if the data are contaminated by a few outliers, these will not 
impact the estimate of the random error which affects most points. 

After applying the LTS method to estimate the random error 
term, we then discard all points with \ri\ > 3 about the LTS fit, 
but only if we are applying the method to a full sample of Aa/a 
values (i.e. the whole VLT or Keck sample, or a combination of the 
two). This is because in small-n fits one does not have much data, 
and so it is not clear whether outliers would become inliers with 
more data. 

If we remove outliers, we then reapply the LTS method to 
check that no more outliers are unmasked, and to re-estimate (Jrand. 
The LTS fit is statistically inefficient because it ignores some good 
data (15 percent for k = 0.85n if all the remaining points are in- 
liers). Therefore, after we discard high residual points, we apply a 
normal weighted least squares fit to the remaining data to estimate 
the parameters and achieve the best possible confidence limits on 
our modelled parameters. 

The benefits of the LTS method can be summarised as follows. 

i) Robust estimate o/(Jrand. If we calculate (Jrand by increas- 
ing it until = 1' then even a single, arbitrarily large outlier 
can increase (Jrand without bound. This is much less likely with the 
LTS method. A more appropriate estimate of (Jrand means that false 
negatives are less likely. 

ii) Robust detection of outliers. In a standard x^ minimisation 
fit, residuals with larger magnitude \ri\ are weighted as r^, which 
distorts the fit towards them. This tends to mask outliers. By dis- 
torting the fit, one might incorrectly decide that some good points 
are in fact outliers. Similarly, the existence of one outlier tends to 
conceal the existence of additional outliers (a masking effect). 

Hi) Objectivity. Manual outlier rejection is often characterised 
as subjective. The rule provided here provides an objective method 
of classifying data points as outliers, thereby removing this objec- 
tion. 

iv) More robust parameter estimates. Even a few outliers can 
substantially distort the fit. This biases parameter estimates away 
from their underlying values. We are interested in the underlying 
values, not the values given by a blind least squares fit. The rate of 
false positives (detections of Aa/a 7^ when Aa/a — 0) should 
also be decreased, as false positives can be caused by outliers. 
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4 VLT RESULTS 

We present here the results of our analysis of the VLT MM ab- 
sorbers. The values of Aa/a for each absorber are given in ta- 
ble |Al| The estimated additional error term, drand, depends on the 
model used, and we give values of (Trand used below according to 
the model in question. The key for the transitions given is defined 
in table [T] 

The frequency with which certain transitions are fitted is also 
given in table [T] We show the distribution of observed wavelengths 
for certain representative MM transitions in figure[6] In figure]?] we 
show the relationship between the g-coefficients and the observed 
wavelength of the transitions used. This figure demonstrates that, 
although q is correlated with wavelength for the low-2 Fe Il/Mg II 
combination (and therefore this combination is susceptible to low- 
order wavelength scale distortions), when the full MM sample is 
considered there is a wide variety of combinations of q and ob- 
served wavelength used. This means that the MM method is gen- 
erally resistant to simple wavelength scale distortions. We give the 
distribution of the signal-to-noise ratio (SNR) for the spectral re- 
gions used in the VLT sample in figure|9] 

Our VLT results are summarised as samples 3, 5 and 6 in table 
|2] We show the values of Aa/a for the VLT sample against redshift 
in figure[8] 

Appendix [C] (in figures |C1| and \C2\ gives two examples of 
the Voigt profile fits for the VLT absorbers. The fits to all of the 
absorbers and an ASCII version of table |Al| can be found atjhttpT] 
[7 /astronomy . swin . edu . au/ ~mmurphy/pub . html| 



4.1 Weighted mean for the VLT Aa/a values 

The LTS method applied to a weighted mean model indicates that 
the z = 1.542 absorber toward J000448-415728 is an outlier, 
with a residual of 4.2cr about the LTS fit, and so we remove this 
Aa/a value. If we do not remove this value, the weighted mean 
after increasing errors is Aa/a — (0.154 ± 0.132) x 10~^, with 
Xl = L06 (cr.and = 0.951 X 10-5). 

After removing this value when calculating our weighted 
mean result for the VLT sample, a weighted mean fit with our raw 
statistical errors yields Aa/a = (0.229 ± 0.095) x 10"^, with 
xt = 1.78. Applying the LTS method to this data set yields a ran- 
dom error estimate of (Jrand = 0.905 x 10"^. 

After accounting for this extra random error, the weighted 
mean becomes Aa/a = (0.208±0. 124) x 10"^ with = 0.99. 
This result differs from that of [Murphy et 31.^2004} at the ~ 4. To- 
level. 



4.1.1 Distribution of Aa/a values with redshift and validity of a 
weighted mean model 

In the bottom panel of figure [s] we show binned values of Aa/a 
plotted against redshift for the VLT sample. For z < 1.5, 3 of the 
5 binned points fall in the region Aa/a < 0. For z > 1.5, 6 
of 7 points in the binned plot fall in the region Aa/a > 0. This 
trend with redshift is different to that seen in fig. 6 of [Murphy et al.| 
( [2004^ ; for z < 1.6, all 7 points fall in the region Aa/a < 0, 
whereas for z > 1.5 all 6 points also fall in the region Aa/a < 0. 
The apparent change in sign of Aa/a with z in the VLT sample 
suggests that a weighted mean model is not a good description of 
the VLT data. 
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Figure 6. Distribution of the observed wavelength of certain representative 
transitions utilised in the MM fits for the VLT sample. The vertical scales 
alternate. 
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Figure 7. Relationship between g-coefficients and observed wavelength for 
all utilised transitions in all absorbers in the VLT sample. Although the low- 
z Fe II/Mg II combination is sensitive to low-order wavelength distortions 
because the g-coefficients for this combination are correlated with wave- 
length (see fig. 1 of Murphy et al.|2003 1, one can see that for the full sample 
there is little correlation between observed wavelength and q, making the 
MM method resistant to systematics when many absorbers at different red- 
shifts are used. 
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Figure 8. Values of Aa/a for the VLT sample. The top panel shows all values of Aa/a with error bars increased in quadrature with o"iand = 0.905 X 10~^. 
The midel panel shows the same data as the top panel, with with the vertical range restricted for better viewing of the higher statistical weight points. Both 
of the panels have been shaded according to a greyscale as the logarithm of the uncertainty estimate, with lower uncertainty points being darker. The bottom 
panel shows binned values of Aa/a where approximately 12 points contribute to each bin. At z > 1.5, Aa/o > 0. This contrasts with |Murphy et al.]j2004| 
(see fig. 6 of that paper), where Aa/« < for z > 1.5. This suggests that a weighted mean model is not a good description of the Aa/a values. 



4.2 Angular dipole fit applied to tiie VLT Aa/a values 

In this section, we fit the angular dipole model of equation[T5]to the 
VLT sample. 

Inspection of the residuals about the LTS angular dipole fit, 
plotted as a function of redshift, reveals no obvious trend for higher 
scatter at higher redshifts and therefore we treat all absorbers 
the same in attempting to estimate (Trand- We again identify the 
z = L542 system toward J000448— 415728 as an outlier, with a 
residual of 4.6(t about the LTS fit, even after increasing the error 
bars. Thus, we remove this system from our sample when fitting 
the dipole model, and re-estimate (Trand = 0.905 x 10~*. We call 
this sample "VLT-dipole". 

Our dipole fit parameters after adding errand ~ 0.905 x 10^^ 
in quadrature to all error bars are: m = (—0.109 ± 0.180) x lO"'', 
A = 1.18 X 10"^ (la confidence limits [0.80, 1.66] x 10"^), 
RA = (18.3 ± 1.2) hr and dec. = (-62 ± 13)°. For this fit, 
= 141.8 and xS: = 0.95. 

To assess the dipole fit compared to a monopole-only 
(weighted mean) fit, we compare a weighted mean fit with errors 
adjusted according to the same (Trand as used for the dipole fit, 
in order to ensure consistency of the data points used when cal- 
culating the statistical significance. As the weighted mean fit has 
X^ = 149.8, the dipole fit yields a reduction in x^ of 7.9 for 
an extra 3 degrees of freedom, when a reduction of ~ 3 would 



be expected by chance. Our bootstrap method yields a significance 
for the dipole-l-monopole model over the monopole-only model at 
the 97.1 percent confidence level (2.19(t), indicating marginal evi- 
dence for the existence of a dipole when considering only the VLT 
data. We demonstrate this fit in figure [70] 

We also give the parameters for a dipole-only (no monopole) 
fit in table |2] 



4.2.1 Ejfect of the choice of the method-of-moments estimator 

In section [T.3.1| we suggested that a method-of-moments estimator 
was preferable in attempting to reconcile Aq/q values from tur- 
bulent and thermal fits. It is legitimate to ask whether our results 
differ if we simply choose that fit (turbulent or thermal) which has 
the lowest xi, instead of applying our method-of-moments estima- 
tor The results for a VLT dipole model if we do this are: (Trand = 
0.928 X 10~^m = (-0.112±0.184) x 10"^ A = 1.15 x lO"'^ 
(Ict confidence limits [0.76, 1.63] x IQ'''), RA = (18.2 ± 1.2) hr, 
dec. = (—62 ± 14)°. The dipole model is preferred over the 
monopole model at the 96 percent level (2.1(t). Thus our choice 
of the method-of-moments estimator does not change the results 
significantly, although arand is mildly larger if we simply choose 



those fits which have the lowest Xv 
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Table 1. The frequency of occurrence for each MM transition in our fits. 
Note that the Mg I A2026 transition is fairly weak compared Mg I A2852, 
and so it is included in few fits. Where we have included Mg I A2052 in 
our fit, and Zn II A2026, is included in our fit, Mg I A2026 will also be 
modelled, although the contribution may be extremely minor. Nevertheless, 
we count this as an occurrence of Mg I A2026, as that transition is included 
in our model. The transition key provides a convenient, short-hand way of 
referring to a particular transition. This key is used in table |Al| 



Transition 


, 1 . 

q (cm ) 


Key 


Frequency of occuiTence 


Mg I A2026 


87 


ai 


3 


Mgl A2852 


86 


a2 


53 


Mgll A2796 


211 


bi 


88 


Mg II A2803 


120 


b2 


86 


Al II A 1670 


270 




60 


Aim A 1 854 


464 


di 


25 


Al III A 1 862 


216 


d2 


25 


Sill A1526 


50 


ei 


57 


Sill A 1808 


520 


62 


31 


Cr II A2056 


-1110 


hi 


21 


Cr II A2062 


-1280 


h2 


15 


Cm A2066 


-1360 


ha 


17 




-1300 


ji 




Fell A1611 


1100 


h 


9 


Fe II A2260 


1435 


h 


12 


Fell A2344 


1210 


ji 


97 


Fe II A2374 


1590 


is 


51 


Fe II A2382 


1460 


36 


100 


Fell A2587 


1490 


jr 


74 


Fe II A2600 


1330 


js 


97 


Mnll A2576 


1420 


h 


13 


Mn II A2594 


1148 


12 


9 


Mn II A2606 


986 


h 


9 


Ni II A 1709 


-20 


ki 


22 


Nill A1741 


-1400 


k2 


24 


Nill A1751 


-700 


k-i 


21 


Ti II A3067 


791 


91 





Ti II A3073 


677 


92 





Till A3230 


673 


93 





Till A3342 


541 


94, 


1 


Till A3384 


396 


95 


1 


Zn II A2026 


2479 


h 


9 


Zn II A2062 


1584 


h 


13 



4.3 VLT results summary 

In this section, we have described the analysis of 154 new MM ab- 
sorbers. The Aa/a values from the VLT sample at z > 1.5 tend to 
show Aa/a > 0, which contrasts with Aa/a < from the Keck 
sample (Murp hy et al.|2004) . We showed that, in the VLT sample, 
an angular dipole model is preferred over a weighted mean model 
at the 2.2a level, which seems to suggest angular (and therefore 
spatial) variations in a. The direction of maximal increase in a is 
found to be RA= (18.3 ± L2) hr and dec. = (-62 ± 13)° under 
a dipole-l-monopole model. We therefore explore the consistency of 
our Aa/a values and model parameters with those derived from 
the same models applied to the Keck sample, and a Keck -I- VLT 
sample, in the next section. 

We have shown that the VLT Aa/a values display excess 
scatter (xt > 1) about the simple models described. This is 




Figure 9. Distribution of the signal-to-noise ratio (SNR) for different spec- 
tral regions in the VLT sample. The histogram gives the distribution of SNR 
values, and the red line gives the cumulative distribution function. The blue 
line indicates the median SNR of 43. The SNR ratio was calculated as 
the inverse of the median of the ratio of the flux array to the error amy for 
the (up to) 1 1 pixels on either side of the spectral region. By "up to", we 
note that for some spectral regions less than 1 1 pixels are available on ei- 
ther side due to gaps in the spectral coverage, or the fact that pixels had been 
clipped out. The en'or array has been modified according to the algorithm 
described in section [TT] and so the SNR estimates given here are more ap- 
propriate than what would result from considering the raw error array. Note 
that the median SNR is somewhat higher than for the sample in [Murphy] 
|etal]{2003^, which was Ri 31. 
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Figure 10. Binned values of Aa/a plotted against angle to the best-fitting 
dipole for the VLT sample. The red, solid line is the model Aa/a = 
A cos{0)+m, and the dashed, blue lines indicate the Icr uncertainty on the 
dipole fit. Statistical eiTors have been increased prior to binning as described 
in the text. The dipole-l-monopole model is preferred over monopole-only 
model at the 2.2a level (see section [4!2) . The parameters for this fit are: 
m = (-0.109±0.180) X 10'^, A = 1.18 x IQ-^ (Icr confidence limits 
[0.80, 1.66] X 10-5), RA = (18.3 ± 1.2) hr and dec. = (-62 ± 13)°. 



likely due to both model mis-specification (from the use of simple 
weighted mean and angular dipole models) as well as unmodelled 
uncertainties. We described in section [33] a number of potential 
random effects which could give rise to excess scatter in the data, 
even if our model for Aa/a were correct. It is difficult to deter- 
mine the contribution of each of these effects to the error budget. 
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Table 2. First table: Constraints from different models for Aa/a from Keck and VLT spectra. / gives a number to identify tlie sample + model. The 
samples are as described in the text. A^abs gives the number of absorption systems used for the fit. m gives the monopole component of the dipole fit 
for dipole models, and the weighted mean of the Aq/q values for a weighted mean model. The RA and dec. columns specify the right ascension and 
declination of the fitted pole. The A column gives the amplitude of the dipole, and the S{A) colunin give s the la confidence limits on the dipole amplitude. 
Raw statistical uncertainties have been multiplied by to obtain confidence limits (Press et al. [1992}, with the exception of the uncertainty on the dipole 
amplitude, which is derived from bootstrapping. A is unitless, except for the r-dipole where A has units of GLyr"^. The column labelled "significance" 
gives the significance of the dipole model over the monopole model, as both a probability and its a equivalent, assessed using a bootstrap method. However, 
for models with no monopole, the significance is given with respect to the null model (Ao/o = 0). Second table: The values of (Tjand used for each 
model. For the Keck Aa/a values, LC and HC refer to low contrast and high contrast respectively. 





Sample + model 




m (10-5) 


RA(hi-) 


dec. (°) 


A (10-5) 


SA (10-5) 


significance 


1 


Keck04-dipole 


140 


-0.465 ± 


0.145 


16.0 ±2.7 


-47 ± 29 


0.41 


[0.29,0.78] 


36 percent (0.5(t) 


2 


# 1 with no monopole 


140 




N/A 


16.4 ± 1.2 


-56 ± 12 


1.06 


[0.82, 1.34] 


72 percent (l.lcr) 


3 


VLT-weighted mean 


153 


0.208 ± 


0.124 


N/A 


N/A 


N/A 


N/A 


N/A 


4 


Combined weighted mean 


293 


-0.216 ± 


0.086 


N/A 


N/A 


N/A 


N/A 


N/A 


5 


VLT-dipole 


153 


-0.109 ± 


0.180 


18.3 ± 1.2 


-62 ± 13 


1.18 


[0.80, 1.66] 


97.1 percent (2. 19o-) 


6 


#5 with no monopole 


153 




N/A 


18.4 ± 1.3 


-58 ± 15 


0.99 


[0.70, 1.37] 


98.4 percent (2.39o-) 


7 


Combined dipole 


293 


-0.178 ± 


0.084 


17.3 ± 1.0 


-61 ± 10 


0.97 


[0.77, 1.19] 


99.995 percent (4.06(t) 


8 


#7 with no monopole 


293 




N/A 


17.4 ±0.9 


-58 ±9 


1.02 


[0.83, 1.24] 


99.996 percent (4.14o-) 


9 


Combined r-dipole 


293 


-0.187 ± 


0.084 


17.5 ± 1.0 


-62 ± 10 


0.11 


[0.09,0.13] 


99.997 percent (4.15(7) 


10 


#9 with no monopole 


293 




N/A 


17.5 ± 0.9 


-58 ±9 


0.11 


[0.09,0.14] 


99.998 percent (4.22(t) 


11 


zl^ dipole, /3 = 0.46 ± 0.49 


293 


-0.184 ± 


0.085 


17.5 ± 1.1 


-62 ± 10 


0.81 


[0.55, 1.09] 


99.99 percent (3.9cr) 



/ 


Sample + model 


(T.andCVLT) (10-5) 


(7,and(KeckLC)(10-5) 


<7,and(KeckHC)(10-5) 


1 


Keck04-dipole 


N/A 





1.630 


2 


#1 with no monopole 


N/A 





1.668 


3 


VLT-weighted mean 


0.905 


N/A 


N/A 


4 


Combined weighted mean 


0.905 





1.743 


5 


VLT-dipole 


0.905 


N/A 


N/A 


6 


#5 with no monopole 


0.882 


N/A 


N/A 


7 


Combined dipole 


0.905 





1.630 


8 


#7 with no monopole 


0.882 





1.668 


9 


Combined r-dipole 


0.858 





1.630 


10 


#9 with no monopole 


0.858 





1.630 


11 


dipole 


0.812 





1.592 



and so we have assumed that all absorbers are affected by the same 
processes, and therefore increased our error bars conservatively in 
quadrature with a (Jrand term. If the extra scatter in the Aq/q val- 
ues is due to inaccuracies in modelling the velocity structure of the 
absorbers, it may be that observations at higher signal-to-noise ra- 
tios and higher resolving powers might help reduce the scatter. On 
the other hand, if the inter-component spacing is comparable to the 
intrinsic line widths then this may not be the case. 

We consider the specific effect of wavelength scale distortions 
on the VLT sample in sections [6.5| and [6/7] and show there how 
such distortions can give rise to extra scatter in the Aa/a values. 



5 COMBINATION AND COMPARISON WITH 
PREVIOUS KECK RESULTS 

We give a summary of all of the results presented here in table|2] 



5.1 Previous Keck results 

We note briefly that we noticed significant wavelength calibra- 
tion problems in the spectrum of Q220 6- 1958 (J220852- 194359) 
from sample 3 of Murphy et al. 12003') for A > 5000 A of the order 
of^ 



10 kms-^ at the time of that analysis. The two absorbers con- 
tributed by this spectrum were erroneously included in that paper, 
and so we remove them from the sample. 

[Murphy et al.|(2004[ ) divide their high-z sample into two por- 
tions, a high-contrast sample and a low-contrast sample. The high- 



contrast sample was defined by 27 absorbers where there were 
significant differences between the optical depth in the transitions 
used. [ Murphy et ar| ( |2003^ give arguments as to why this might be 
expected to generate extra scatter in the Aa/ot values. Due to the 
fact that many of the high redshift {z > 1.8) absorbers considered 
in |Murphy et al.| ( (2004^ are associated with damped Lyman-a sys- 
tems, this effect manifests itself as extra scatter in the Aa/a values 
about a weighted mean at high redshifts. For the VLT sample, we 
note that there is no evidence for excess scatter at higher redshifts 
compared to lower redshifts. 

We can examine the differences between the Keck and VLT 
samples in terms of the prevalence of weak species as follows. 
Firstly, define the following transitions as weak: Mg I A2026, Si II 
A1808, the CrII transitions, Fell AAA1608, 1611, 2260, the Mnll 
transitions, the Nill transitions, the Till transitions and the Znll 
transitions. From the table of the frequency of occurrence of these 
transitions in Murph y et al.H2003[ l, aX z < 1.8 these transitions 
constitute about 3 percent of the total number of transitions used. 
On the other hand, in the VLT sample these transitions consti- 
tute about 13 percent of the sample used. The significantly greater 
prevalence of these weak transitions at low redshifts in the VLT 
sample may explain the lack of evidence for differential scatter be- 
tween high and low redshifts. Effectively, the greater prevalence of 
weak species in the low-z VLT sample may increase the scatter at 
low redshifts in that sample, making any low-z/high-z difference 
appear smaller. We retain the high/low contrast distinction when 
analysing the Keck sample. 
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5.1.1 LTS method applied to the \Murphy et al.\\2004\ results 

If we apply the LTS method to the high contrast sample to estimate 
the extra error needed about a dipole model, we find that an ex- 
tra error term of o"rand = 1.630 x 10~^ is needed. With this extra 
term, = 1.13, indicating that the distribution is mildly leptokur- 
tic (fat-tailed). The low-contrast sample data are already consistent 
under a dipole model with the LTS method (errand = 0). 

We then combine the high-contrast Aa/a values (with error 
bars increased) with the low-contrast Aa/a values to form a new 
sample under a dipole model (equation [Tsj. The LTS method ap- 
plied to this set reveals that the Aa/a values are consistent about 
dipole model. Additionally, xi = 1.04. Nevertheless, we identify 
one possible outlier from this set: the absorber with z ~ 2.84 to- 
wards Q1946+7658, with Aa/a = (-4.959 ± 1.334) x 10"^ 
and remove this absorber from the sample. This point has a resid- 
ual of — 3.6ct about the LTS fit. We refer to this sample as "Keck04- 
dipole". 

A dipole fitted to this sampl^ yields RA = (16.0 ± 2.7) hr, 
dec. = (-47 ± 29)°, and A = 0.41 x 10^^. la confidence limits 
on A are [0.29,0.78] x 10"^. The monopole is m = (-0.465 ± 
0.145) X 10"^ This fit has xl = 0.96. The dipole model is pre- 
ferred over the weighted mean model at the 36 percent confidence 
(0.47cr). 

The monopole offset appears to be significant at the 3.2cr con- 
fidence level, but this is related to the fact that the Keck results 
alone do not clearly support a dipole interpretation. 

For dipole model with no monopole (Aa/a — AcosS), 
the fitted parameters are A = 1.06 X lO"'^ (1(7 confidence lim- 
its [0.82, 1.34] X lO"'^), RA = (-16.4 ± 1.2) hr, dec. = (-56 ± 
12)°. This model is significant at the 72 percent confidence level 
(1.1a). 



5.2 Combined weighted mean 

We create a combined weighted mean fit by combining the VLT- 
dipole sample with the Keck04-dipole sample. The VLT sample 
has had errors increased in quadrature with (Trand = 0.905 x 10^^, 
whereas the Keck high-contrast sample has had errors increased in 
quadrature with errand = 1.743 x 10~'\ The same points identified 
as outliers have been removed. 

This leads to a weighted mean of {Aa/a)^ = (—0.216 ± 
0.086) X 10"^, with xl = 1-03. However, a weighted mean 
model does not appear to adequately capture all the information 
in the data (see figure [T4|. Comparing the weighted mean of the 
z > 1.6 points for both samples yields a simple demonstra- 
tion of north/south difference. For the VLT sample, Aa/a^,(z > 
1.6) = (0.533 ± 0.172) x 10"^, whereas for the Keck sample 
Aq/q„(2 > 1.6) = (-0.603 ± 0.224) x 10"^ The difference 
between these weighted means is significant at the 4cr level. 



5.3 Combined dipole fit 

To create our combined dipole fit, we combine the VLT-dipole sam- 
ple with the Keck04-dipole sample to create the "combined dipole" 
sample, which is our main sample. This sample consists of 293 MM 



We note that in [Murphy et al.|j2003} the reported position of a dipole 
model fitted to that data set was (5 hr, —48°). The RA value quoted there 
is eiToneous due to a typographical error. The value should have been stated 
as (17 hr, —48°). The value reported there in Galactic coordinates (the co- 
ordinate system in which the computations were done) was correct. 



absorbers. Importantly, both of these sets exhibit no jr^ j 3 residu- 
als, and thus a combined fit is unlikely to exhibit any large residuals 
provided that both data sets are well described by the same model. 

For an angular dipole fit to these Aa/a values (Aa/a — 
Acose + m), we find that m = (-0.178 ± 0.084) x 10"^ 
A = 0.97 X 10"^ (la confidence limits [0.77, 1.19] x 10"^), 
RA = (17.3 ± 1.0) hr, dec. = (-61 ± 10)°, with x^ = 280.6 
and xl = 0.97. A weighted mean (i.e. monopole only) fit to the 
same Aa/a values and uncertainties yields x^ = 303.8, and so a 
dipole model yields a reduction in x^ of 23.2 for an extra 3 free pa- 
rameters. With our bootstrap method, we find that the dipole model 
is preferred over the weighted mean fit at the 99.995 percent con- 
fidence level (4.06a), thus yielding significant evidence for the ex- 
istence of angular variations in a. Using the method of [Cooke &| 
[Lynden-Bell[ j2010^ , the significance of the dipole is found to be 
4.07a. 

Importantly, the combination of the Keck04-dipole Aa/a val- 
ues with the VLT-dipole Aa/a values yields xl ^ ^ about a dipole 
model. If inter-telescope systematics were present, we would ex- 
pect the combination of the Keck and VLT data to yield a xl that 
is significantly greater than unity under the dipole model, despite 
xl being ~ 1 when that model is fitted to the samples individu- 
ally. Thus, there is no significant evidence based on x^ th^t inter- 
telescope systematics are present. 

the values of Aa/a for both Keck and 



We show in figure 



11 



VLT against the best-fitting dipole model. We give binned values 
there, which yields a visual demonstration of the dipole effect. We 
also give there a plot of the standardised residuals about the fit, 
which demonstrates that the fit is statistically reasonable. We also 
show binned values of Aa/a for the Keck, VLT and combined 
samples in figure [14] We show an unbinned version of these data 
for \ Aa/a\ < 5 x 10"^ in figure 
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For a model with no monopole (Aa/a = ylcosO), the 
fitted parameters are A = 1.02 X lO"'^ (la confidence lim- 
its [0.83,1.24] X 10"^), RA = (17.4 ± 0.9) x 10"^ dec. = 
(—58 ± 9)°. This model is significant at the 99.996 percent level 
(4.14a). 

In figure [T2] we show the confidence limits on the dipole lo- 
cation for the model Aa/a = AcosO for the Keck, VLT and 
combined samples. The individual symbols illustrate the weighted 
mean of Aa/a along each sightline. 

There are several significant points to consider from these re- 
sults: 

/) The dipole is statistically significant. Even after accounting 
for random errors in a conservative fashion, the statistical signif- 
icance of the dipole is greater than 4a. This is strong statistical 
evidence for angular and therefore spatial variation in a. 

ii) Dipole models fitted to the Keck and VLT Aa/a values 
yield consistent estimates for the pole direction. This is important, 
and would be very surprising if one assumes that a dipole effect 
is not present. If two different systematic effects were operating in 
each telescope so as to produce a trend in Aa/a, then: a) it is un- 
likely that these effects would be correlated with sky position (the 
most likely systematics relate to problems with wavelength calibra- 
tion), and b) even if systematic effects existed in both telescopes 
which were correlated with sky position, it is very unlikely that 
such effects would occur in such a way as to yield very consistent 
estimates of the dipole position between the two telescopes, with 
a similar amplitude, particularly when the two telescopes are in- 
dependently constructed and separated by ~ 45° in latitude. Any 
attempt to ascribe the observed variation in a to systematics must 
account for the good alignment of the dipole vectors from dipole 
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models fitted independently to the Keck and VLT samples. Note 
that telescope or instrumental systematics which depend only on 
wavelength cannot produce observed angular variation in a for a 
sufficiently large sample of absorbers. 

Hi) The VLT and Keck Aa/a values appear consistent near 
the equatorial region of the dipole. From the middle panel of fig- 
ure [TT| both the VLT and Keck results show large variation from 
Aq/q = near the pole (6 = 0°) and anti-pole (9 = 180°) 
of the dipole, but show much less variation in the equatorial re- 
gion (0 = 90°). So, at least visually, the Keck and VLT points 
are not inconsistent in the region where they overlap. This issue is 
addressed quantitatively in the caption to figure [T3] 

iv) The dipole effect is not being caused by large residual 
points. The bottom panel of figure [TT| clearly shows that there are 
no |ri| > 3(T points present. 



5.4 Potential effect of differences in atomic data and q 
coefficients 

If the atomic data or q coefficients we used were significantly dif- 
ferent to those used by |Mutphy et al.| j2004| l, this could spuriously 
create differences in Aa/a between VLT and Keck. This has the 
potential to mimic spatial variation in a. To check the influence of 
this, we re-fit the VLT spectra using the same atomic data used by 



Murphy et al. and then combine the Aa/a values with the Keck 



values. Where we use transitions that were not available to lMur-l 
|phy et ar| (e.g. Mn II and Ti II) we make no modification to the 
atomic data or q coefficients. The frequency of occurrence of these 
transitions in the sample is small and therefore this is of little conse- 
quence. When we proceed in this way, the parameters for the model 
Aa/a = Acose + m are: A = 0.97 x 10"^ RA = (17.5 ± 
1.0)hr, dec. = (-60±10)° andm = (-0.168±0.084) x lO'^ 
The significance of the dipole-l-monopole model with respect to the 
monopole-only model is 4. ISa. We conclude that the impact of any 
variations between atomic data or the q coefficients used for our fits 
and those used by |Murphy et al.| ( (2004[ > is negligble. 



5.5 Alignment by cliance between Keck and VLT 

One can pose the question: "Given the distribution of sightlines 
and values of Aa/a in each sample, what is the probability of ob- 
serving alignment as good or better than that observed between the 
Keck and VLT samples by chance?" To assess this, we undertake a 
bootstrap analysis, where at each bootstrap iteration we randomly 
reassign the values of Aa/a in both the Keck and VLT samples 
to different sightlines within those samples, keeping the redshifts 
of the absorbers fixed. That is, we do not mix the two samples. We 
then calculate the best-fitting dipole vectors for each sample, and 
calculate the angle between them. We then assess over many itera- 
tions in what percentage of cases is the fitted angle smaller than the 
angle for our actual data. 

For our actual Keck and VLT samples, the angle between the 
fitted dipole vectors is 24 degrees, and the chance probability is 
~ 6 percent. We show the results of this bootstrap analysis in fig- 
ure [TS] Thus, it seems unlikely that inter-telescope systematics are 
responsible for the observed effect. The good consistency between 
the results also qualitatively supports the notion that the measured 
effect is real. 



5.6 Low-2 vs high-z 

We divide our sample into low-z and high-z absorbers to examine 
the contribution of the different redshifts to the dipole detection. 
Although there is no clear delineation between which transitions 
are fitted for a given redshift, we can generally say that the \ow-z 
sample is dominated by the Mg/Fe combination, that intermedi- 
ate redshifts display a wide range of transitions, and that high red- 
shift systems are dominated by the Si Il/Al Il/Fe II A1608 combina- 
tion. In particular, Mg II, Mg l A2852 and the Fe II transitions with 
Ao > 2200 are not generally used when fitting absorbers at high 
z because they are either beyond the red cut-off in the observed 
spectral range, or the transitions are affected by sky absorption or 
emission. 

If the observed dipole effect was due to chance or caused by 
a systematic effect which affects some combination of transitions, 
then we would not expect dipole fits to absorbers from high and low 
redshift to yield the same location on the sky. Conversely, if dipole 
models fitted to high and low redshift samples point in a similar 
direction, this lends support to the dipole interpretation of the data. 

To check the consistency of the estimates for low and high red- 
shift samples, we cut the data into a z < 1.6 sample (low-z) and 
a. z > 1.6 sample (high-z). This divides the data approximately 
in half, with 148 points in the low-z sample and 145 in the high- 
z sample. We show in figure [16] the results of separate fits to the 
low-z and high-z sample, and demonstrate that they yield consis- 
tent estimates of the dipole location. We give the parameters to the 
model Aa/a = A cos{&) + m in tablejs] In particular, the dipole 
vectors are separated by 13 degrees on the sky. 

To assess the probability of obtaining alignment as good or 
better than is seen by chance, we use a bootstrap method where we 
randomly re-assign values of Aa/a to different absorbers within 
the z < 1.6 and z > 1.6 samples. That is, we do not allow values 
of Aa/a to mix between the z < 1.6 and z > 1.6 samples. Given 
the distribution of Aa/a values, absorbers and sightlines in each 
sample, the chance probability is ~ 2 percent. Given that the tran- 
sitions used at low and high redshift are significantly different, this 
consistency further supports the dipole interpretation of the data. 
It is also clear that the dipole signal is significantly larger at high 
redshift, although the low redshift sample contributes. 

There is no significant evidence for a high-z monopole, but 
the low-2 monopole is significant at the 3.6a level. We discuss the 
significance of the low-z monopole in section[7!8l 



5.7 Joint probability 

The probability of obtaining alignment between the dipole vectors 
for dipole models fitted to the Keck and VLT samples separately as 
good or better than is seen by chance is about 6 percent. The chance 
probability of obtaining alignment between the dipole vectors for 
dipole models fitted to low- and high-redshift samples is about 2 
percent. Through a bootstrap method we have calculated the joint 
probability of obtaining alignment that is at least as good as seen 
for both of these conditions by chance, and it is « 0.1 percent. 



5.8 Significance of the low-z monopole 

In section [S^l we noted that the low-z sample shows evidence for 
a statistically significant monopole at the 3.6a level. In figure [14] 
the existence of the monopole in both the Keck and VLT samples 
can be seen at low z: the mean value of Aa/a at z < 1.6 is sig- 
nificantly negative in both samples. Note in particular the top panel 
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Figure 11. The top panel shows the combined results for the Keck and VLT samples, plotting Aa/a against angle from the fitted dipole location for the 
combination of the Keck and VLT Ao / a values, binned together. The middle panel shows the data for the top panel, but separating out the different telescopes. 
Points in the top panel contain approximately 25 absorbers per bin, whereas the points in the middle panel contain approximately 12. In the middle and bottom 
panels. Keck points are indicated by squares and VLT points are indicated by circles. The model shown (red, solid line) is Aa/a = j4cos(0) + m. The 
parameters for this model are: m = (-0.178 ± 0.084) X lO'^, A = 0.97 X 10"^ (lo- confidence limits [0.77, 1.19] X 10"^), RA = (17.3 ± 1.0) hr, 
dec. = (—61 it 10)°. The dashed, blue lines indicate the Icr uncertainty on the fit, including the uncertainty in determining the position of the dipole, the 
amplitude of the dipole and the monopole value. In both the top and middle panels, the dotted horizontal line indicates the monopole value. The bottom panel 
indicates the standardised residuals (r ^ = [data — model] /error) about fitted model. The vertical dotted line in all panels show 90° from the dipole pole. 
The presence of no points with \ri\ > 3 indicates that the fit is not being dominated by a small number of large residual points. Note that because of the 
direction of the dipole axis in relation to the positions of the telescopes, many of the Aa/a values fall near the equatorial region of the dipole (0 = 90°). 
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Figure 12. Sky map in equatorial (J2000) coordinates showing Aa/cj values aggregated for each sightline to a quasar, where the value of A«/ci along a 
sightline is taken as the weighted mean of the Aa/o values for that sightline. This necessarily obscures any z dependence that may be present. The symbols 
show the residual of the value of Aa/o with respect to the null model, that is r = {l\a/ a) / a ^. Larger symbols indicate greater deviation. Circles 
represent Keck quasars, squares represent VLT quasars and triangles represent quasars common to both samples. The colour coding indicates the difference 
between the sightline value of Ao / a and zero. The green shaded region indicates the la confidence limit for a dipole of the form Aa /a = A cos{©) fitted 
to values of Aa/a from all Keck absorbers. The blue and red regions show the same for dipole models fitted to the VLT and VLT+Keck sample respectively. 
The anti-pole is marked with an "A". The blue, dashed line indicates the equatorial region of the dipole. The grey shaded region indicates the galactic plane, 
with the bulge indicating the galactic centre. The dipole is visible as more large, blue points near the pole and more large, pink points near the anti-pole. 



Table 3. Pai'ameters for the model Aa/a = Acos(0) + m for z < 1.6 and z > 1.6 samples. The column "5A" gives la confidence limits on A. The 
column labelled "significance" gives the significance of the dipole model with respect to the monopole model. Although it is clear that most of the significance 
comes from the 2 > 1.6 sample, the z < 1.6 sample also contributes. Additionally, a dipole model for the z < 1.6 sample points in a similar direction to that 
of the z > 1.6 sample. 



Sample 


A (10-^) 


<5A(10--0 


RA (hr) 


dec. (°) 


»n(10"^) 


significance 


2 < 1.6 
2 > 1.6 


0.56 
1.38 


[0.38,0.85] 
[1.12, 1.74] 


(18.1 ± 1.8) 
(16.5 ± 1.4) 


(-57± 22) 
(-63 ± 11) 


(-0.390 ±0.108) 
(0.097 ± 0.138) 


lAa 
3.5a 



in figure [14] showing the combined results from both telescopes, 
where the mean value of Aa/a shifts from approximately zero at 
high z to negative values at low z. 

An obvious question is whether the monopole arises from 
one of the Keck or VLT samples. For a dipole-l-monopole model 
(Aa/a = Acos(0) + m), the Keck sample yields a 2 < 1.6 
monopole of m = (—0.404 ± 0.171) x 10"^, which differs from 
zero at the 2.4o- level. However, the same model fitted to the VLT 
Aa/a values at z < 1.6 yields m = (-0.373 ± 0.295) x 10"^ 
differing from zero at the 1.3cr level. There are three important 
considerations from these values: /) Both data sets yield consis- 
tent monopole values for Aa/a at low redshift, differing at only 
the 0.09a level. Therefore, whatever is generating the monopole 
appears to affect both the Keck and VLT samples, ii) Because there 
is no significant difference between the monopole values in the 
Keck and VLT samples, the monopole cannot be responsible for 



mimicing angular variation in a across the sky. Hi) Additionally, the 
dipole signal is dominated by Aa/a values at 2 > 1.6; there, the 
dipole-l-monopole model is significant at 3.5(7 over the monopole- 
only model. As such, the significance of the dipole signal is not 
strongly affected by the presence of the low-2 monopole. That is, 
to first order, the dipole and low-2: monopole signals appear decou- 
pled. 

There are several possible explanations for this, and we dis- 
cuss each of them in turn: 

i) Telescope systematics. Perhaps the most obvious and im- 
portant concern is that there is some simple, wavelength-dependent 
systematic effect between the Keck or VLT telescopes or spectro- 
graphs. However, such wavelength-dependent telescope systemat- 
ics seem difficult to support given the inter-telescope consistency 
of the best-fitting dipole directions on the sky. Nevertheless, such 
effects are difficult to rule out entirely. Therefore, we address this 
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Figure 13. Aa/a against the angle from the fitted dipole location under the model Aa/a = A cos(0) + m for the VLT and Keck Aa/a values. In contrast 
to figure [77] these data are not binned. Blue circles are VLT absorbers and pink squares are Keck absorbers. Error bars have been omitted; instead, larger 
points indicate points with greater statistical weight, according to the key provided. The precision includes the effect of o-j-and ■ The dipole trend is visible as 
the presence of more and larger points in the upper left and lower right quadrants. The visual cluster of points at © < 47° is due to 4 quasars which contribute 
14 values of Aa/a (2 points not shown because they lie outside the vertical range of the plot). One can investigate the consistency of the VLT and Keck 
Aa/a values in the region near the dipole equator (defined here as 80° < © < 100°) by comparing the weighted mean of the Aa/a values. In this case, 
Aq/«i„(VLT) — Aa/au,{Keck) = (0.32 it 0.19) x 10^^, giving no significant evidence for a difference between the two samples. In this region, the 
VLT sample contributes 39 points and the Keck sample contributes 43 points. The difference here is calculated so as to include the effect of (Tjand- 



potential concern in detail in section[6|by directly comparing spec- 
tra of quasars taken with both the VLT and Keck telescopes. 

ii) Errors in the laboratory wavelengths. EiTors in the labora- 
tory wavelengths of transitions which feature predominantly at low 
redshifts could cause a statistically significant monopole at low red- 
shifts. However, this seems particularly unlikely. The Mg I/II wave- 
lengths have been accurately measured on an absolute scale gener- 
ated using a frequency-comb calibration system. The Fe II wave- 
lengths used at z < 1.6 have also been precisely measured (the 
A1608, 1611 transitions are more difficult to measure accurately, 
but these transitions are used infrequently at low redshifts due to 
their short rest wavelengths). For instance, the absolute velocity 
uncertainty in the Fe II A2382 transition is « 14ms~^, which is 
significantly smaller than the ~ 82 ms~^ which would be needed 
to generate a monopole value of —0.39 x 10~*. This implies a sys- 
tematic error some six times larger than the existing error budget, 
which seems unlikely. Additionally, the relative wavelength scales 
of the different experiments which measured the transitions used at 
lower redshifts are likely to be significantly better than this. 

Hi) Time evolution of a. The functional form for variation of a 
(if a varies) is unknown. Recent observations confirm the apparent 
acceleration of the universe at late times ( [Astier et al.|2006l l, for 
which dark energy is posited as an explanation. For z < 0.5, dark 



energy dominates over matter and radiation jRiess et al.|2004^ . If 
a couples to dark energy, then late-time evolution of a might be 
possible. Monotonic evolution of a cannot by itself be an explana- 
tion for a \ow-z monopole, because this would imply that Aa/a 
should approach zero for 2 — >■ 0, with the greatest divergence of 
Aa/a from zero at high redshift. If a oscillates with time then a 
pattern such as that observed in figure[74]could arise. However, this 
would require the period of oscillations to be ~twice the age of the 
Universe, with the present day at a node of the oscillation, in order 
to obtain (Aa/a^>i,6) ~ 0, (Aa/Qo.2<2<i.6) ~ -0.4 x lO"'^ 
and Aa/az=o = 0. This may be possible, but this case seems 
rather contrived. 

/ V ) Dependence of a on the local environment. If the value of a 
depends on the local environment (e.g. matter density, gravitational 
potential, or gradient of the gravitational potential) then this could 
produce an offset between the value of a measured in the quasar 
absorbers and the value measured in the laboratory, even as z — )■ 0. 
If this was the case, we would expect a similar magnitude monopole 
to also be present at high redshift, which is not seen. It is possible 
that a combination of the time evolution of a and dependence of 
a on the local environment could explain the low-z monopole, but 
this requires two different mechanisms. Additionally, in this cir- 
cumstance the magnitude of the environmental dependence must 
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Figure 14. Binned values of Aa/« by redshift in the VLT-dipole sample (bottom panel, circles), the Keck04-dipole sample (middle panel, squares) and the 
combination of the two (top panel, triangles). The value of Aa/o for each bin is calculated as the weighted mean of the values of Aa / a from the contributing 
absorbers. The statistical errors for certain points have been increased prior to binning, as described in the text. Note that for z > 1.5, the Keck values generally 
indicate Aa/a < 0, whereas the VLT values indicate Aa/cx > 0. As Keck is located in the northern hemisphere, and VLT is in the south, this is a rough 
visual demonstration of the dipole effect. However, given the overlap between the samples, the proper procedure is to directly fit a dipole (see figure [Tl} . 
Interestingly, both the VLT and Keck values seem to to support Ack/q < for 2 < 1.5. This effect is considered in more detail in sectionp^ 
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Figure 15. Results of the bootstrap analysis described in section [53] to as- 
sess the probability of obtaining alignment of the dipole vectors between 
the VLT and Keck samples as good as we have observed by chance. The 
vertical red line shows the angle between the dipole vectors of our sample 
(24 degrees). The area to the left of the red line indicates the probability of 
interest, namely 6 percent. 



be very similar to the magnitude of the time evolution from z ~ 4 
to 2 = in order to obtain the observed distribution of Aa/a with 
z, requiring significant fine-tuning. 

v) Non-terrestrial isotopic abundances in the absorbers. The 



isotopic spfitting scales as AtJi oc ujo /rn^, where is the mass of 
the species under consideration. Mg is the lightest atom used in the 
MM method, and therefore the isotopic splitting for the Mg transi- 
tions is relatively large. If the abundance of the three Mg isotopes 
differs significantly in the quasar absorbers to terrestrial values, this 
would mimic a change in a. The low-z sample is dominated by the 
Mg Il/Fe II combination, which is particularly sensitive to the ef- 
fect of differences in the abundance of the Mg isotopes (Murphy| 
et al.|2001b) . At higher redshifts, the effect on Aa/a of isotopic 
abundance variations in other light species (e.g. Si) has been shown 
to be negligible ( |Murphy et al.|2001b[ l. The fact that we do not see 
a significant monopole signal at high redshift lends additional cre- 
dence to this interpretation. 

On balance, this last possibility - evolution in the abundance 
of the Mg isotopes - seems like the most likely of the explanations 
above. We therefore explore this isotopic abundance effect in more 
detail in sectionl6l2l 



5.9 Iterative clipping of potentially outlying Aa/a values 

We have attempted to be conservative in presenting our results 
when accounting for extra scatter in the Aa/a values about a 
model by adding a term, arand, in quadrature with the error bars. 
This effectively functions as an interpolation between a fit where 
the error bars are believed to be correct and an unweighted fit, 
where the error bars are unknown. 

However, another option is to assume that the statistical er- 
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Figure 16. Sky map in equatorial (J2000) coordinates showing the 68.3 
percent (Icr equivalent) confidence limits of the location of the pole of the 
dipole for the z < 1.6 combined sample (blue region) and 2 > 1.6 com- 
bined sample (red region) under the model Aa/a = A cos(0) + m. The 
location of the CMB dipole and anti-pole are marked as Pcmb ™d Acmb 
respectively for comparison jLineweaver| 1997) . This figure demonstrates 
that the low-2: and high-z absorbers produce consistent estimates of the 
dipole location, despite generally using significantly different combinations 
of transitions. The dipole vectors for the z < 1.6 and z > 1.6 sample are 
separated by 13 degrees. The probability of getting alignment this good or 
better by chance is 2 percent. 

ror bars on most Aa/a values are a good representation of the 
total uncertainty for those absorbers, and then remove points one- 
by-one ("clipping") until xi = 1. In our sample it is difficult to 
determine to what extent different random processes affect differ- 
ent absorbers, and therefore to determine to what extent clipping 
is justified. Adding some (Jrand in quadrature with all Aq/q val- 
ues, as we have done, is a conservative option. Nevertheless, we 
explore the effect of clipping here to investigate the robustness of 
our results to the removal of Aq/q values. 

Traditionally, data clipping involves iteratively removing the 
point with the largest residual and then re-fitting. However, for 
the reasons given in section [333] this has the potential to incor- 
rectly remove points. Therefore, we use a modified method. At 
each iteration, we calculate the LTS fit using the model Aq/q = 
A cos O + m to the Aq / q values with their raw statistical errors, 
and then remove the point with the largest residual about the fit. 
However, we choose = n — 1 in this case. Effectively, at each 
stage, we want to identify only one point to remove, and therefore 
it makes sense to calculate a fit to n — 1 points. We then calculate 
a weighted fit, using only the statistical error bars, and calculate 
the significance of the dipole model. For efficiency of calculation, 
we avoid bootstrapping, and so we use the method of |Cooke &| 
|Lynden-Be ll ( 2010) to calculate the significance of the dipole fit 
with respect to the monopole fit. We then repeat the process. At 
any iteration, if > 1, we multiply all entries of the covariance 
matrix by in order to account for excess scatter about the model. 
If < 1> we do not adjust the covariance matrix. 

Initially, one expects the significance of the fit to improve, as 
one discards a few points which are not consistent with the gen- 
eral trend of the fit. Eventually, one will remove enough Aq/q 
values that the significance must decline. If the significance de- 
clines rapidly, this implies that the dipole effect is dominated by 
a few points. Conversely, if the significance of the fit is sustained 
or improved for the removal of small fractions of data (e.g. ~ 10 
percent), this qualitatively implies robustness of the dipole effect. 

We show in figure[T7]the results of this process. We find that 
we must remove large numbers of absorbers to destroy the signifi- 



cance of the dipole. In particular, the significance does not decrease 
rapidly with the number of Aq/q values clipped initially, suggest- 
ing that the observed dipole effect is not being caused by a few 
outlying points. If we clip until = the significance of the 
dipole is almost 7cr. 

In figure [Is] we show the effect of clipping Aq/q values on 
the location of the dipole. One expects that if the dipole effect is 
real, then the position of the dipole should not change dramatically 
with the removal of small amounts of data (that is, A0 = 0i — Oo 
should be small). To assess how likely it is that this seemingly re- 
stricted path is typical for our distribution of data, we apply a boot- 
strap method to generate and iteratively trim 300 new samples, and 
examine the distribution of AG at each point. We cannot use a 
traditional bootstrap, which resamples the data with replacement, 
because how the data is trimmed depends crucially on the distribu- 
tion of residuals. Therefore, we resample the residuals of the fit to 
generate new samples. To do this, we use the following process to 
generate one sample: /) calculate the model prediction for each ab- 
sorber given the model, Pi = A cos(0) + ra; ii) calculate the resid- 
uals about the fit for each absorber, = (Aa/ui— Pi)/ ai; Hi) ran- 
domly reassign the calculated r,; to different absorbers, generating 
r'j; iv) generate a new set of Aq/q values as Aa/a'j = pj +rjaj. 
In this way, we generate new values of Aq/q which represent dif- 
ferent possible realisations of our sample where the actual distri- 
bution of residuals is preserved. This is demonstrated in figure [TS] 
We see that the bootstrapped samples do not wander very far even 
when much of the data is removed (A6 < 20°). 

To contrast this with the effect on a random sample, we also 
show in figure [TS] the effect of trimming random samples of data. 
To do this, we generate 300 new samples by randomly reassigning 
values of Aq/q to different sightlines, and iteratively trimming un- 
der the model Aq/q = A cos(0) +m. We see here that our actual 
sample is not typical of the random samples, therefore suggesting 
that the actual sample is significantly dissimilar to random samples. 

5.10 Removal of spectra 

A further question one might ask is how sensitive our results are to 
the inclusion of particular spectra. We would like to know whether 
the dipole result could be dominated by a small number of spectra 
which, if removed, would destroy the result. 

We therefore explore this question through a jack-knife 
method, where we remove one quasar at a time and recalculate 
the statistical significance of the fit. We show the results of this 
exploration in figure [T9] The figure clearly demonstrates that, un- 
surprisingly, our result is not due to a single quasar spectrum. We 
extend this in figure |20] to show the effect of removing 5 spec- 
tra at random. We chose 5 as a number somewhat larger than 1, 
and also in order to potentially include the cluster of 5 quasars at 
RA ~ 22hr, dec. ~ —45°, where all of these sightlines demon- 
strate Aq/q > 0. Under this circumstance,the probability of ob- 
taining a dipole result which is insignificant (< Scr) is small. This 
suggests that the dipole effect is not being created by a small num- 
ber of spectra. 

5.11 Translation from angular variation to a physical model 
including a distance measure 

We now explore simple phenomenological parameterisations of the 
dipole effect which attempt to account for distance dependence. In 
all of these models, the same Aq/q values identified as outliers 
have been removed from consideration. 
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Figure 17. Effect of iteratively clipping Aa/a values on the statistical sig- 
nificance of the dipole for the combined sample, as described in section 
|5.9| for the model Aa/a = j4cos(©) + m. The vertical axis shows 
the statistical significance of the dipole as determined by the method of 
|Cooke & Lynden-Bell|j20l"o| given in terms of a (solid line) and ^t that 
point (blue, dotted line). A dashed horizontal line is drawn at 3cr for refer- 
ence. The vertical red (dashed) line indicates the point at which our clipping 
method reduces to below unity. We note that we have to remove more 
than 40 percent of absorbers before the significance of the detection drops 
to about 3(T. As there is no good reason to remove so much data, this implies 
that the dipole result is not being drive by a few outlying Aa / a values. The 
actual significance given here is probably overstated compared to the "true" 
significance, given that no attempt has been made to account for systematic 
or random eiTors. 



5.11.1 z-dipole fit 

To model potential distance dependence directly with the observ- 
able quantity, z, we fit a power-law relationship of the form 



Aa/a = Cz^^ cos(O) + m, 



(18) 



for some /3 and amplitude C to the combined Keck -I- VLT data, 
where the Aa/a values have been increased in quadrature by the 
amounts set out for sample 11 in table [2] This sample is the "z^ 
dipole" sample. 

We use the Levenberg-Marquardt algorithm ( [Press et al.|1992^ 
to fit equation 



dipole sample. This fit yields RA — 
62 ± 10)°, C = 0.81 (Icr confidence 



to the 2 

(17.5 ± 1.0) hr, dec. = ( 

limits [0.55, 1.09] x 10"^), m = (-0.184 ± 0.085) x lO""^ and 
/3 = 0.46 ± 0.49. The fact that the amplitude grows as a low power 
of z, and the fact that it is statistically consistent with zero, is the 
reason that the approximation A Cz" yields reasonable results 
earlier We show the results of this fit in figure [21] This dipole -1- 
monopole model is statistically preferred over the monopole-only 
model at the 99.99 percent confidence level (3.9cr). The reduction 
in significance relative to the angular dipole model occurs as a re- 
sult of the uncertainty in determining /3, but is relatively small. 



5.11.2 r-dipolefit 

Another plausible alternative is to try to relate the amplitude of the 
dipole to some explicit distance metric. For simplicity, we use the 
lookback-time distance. This is defined by r = ct, where c is the 
speed of light and t is the lookback time to the absorber. Thus, we 
try a fit of the form 
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Figure 18. Effect of iteratively clipping Aa/a values on the location of the 
dipole, as described in section [s!9] The vertical axis shows the deviation of 
the fitted angle from the untrimmed model (A0) as a function of the per- 
centage of absorbers removed. This figure compares the results of trimming 
on our actual Aa / a values, bootstrapped samples designed to emulate our 
data, and random samples. Actual data: The solid black line shows the re- 
sults of trimming for our l\a/cx values. If the fit is stable and not due to 
the presence of a small number of highly significant points, we expect to 
see that A0 should not grow rapidly with the amount of data removed. 
This is what is seen. Bootstrapped samples: The dashed lines show the Itr 
range for 300 bootstrap samples (generated as described in the text). This 
shows the typical range of variation at fraction of absorbers removed given 
distribution of sightlines, values of Aa/a, statistical errors and distribu- 
tion of residuals in the sample. The region given reflects the Icr range for 
the bootstrapped samples at each point; each individual sample may wander 
substantially more than is indicated by this range, and so the deviation of 
the path for our actual sample outside the region is not indicative of any 
problem with our Aa/a values. Random samples: The blue, dotted lines 
show the 1(T range for 300 samples where we have randomised /\a/a over 
the sightlines. We see that A0 in this case grows rapidly with increased 
trimming for these samples. Our real sample does not do this, which sug- 
gests that our real sample is significantly dissimilar from a random sample. 



To calculate lookback times, we use the standard ACDM model (A 
Cold Dark Matter), with parameters given by the 5-year WMAP 
(Wilkinson Microwave Anisotropy Probe) results l |Hinshaw et al.| 
|2009). We note that this calculation is derived from the FLRW 
metric, which assumes isotropy of the universe. Our model im- 
plies anisotropy of the universe, and therefore use of the FLRW 
(Friedmann-Lemaitre-Robertson-Walker) metric is strictly incor- 
rect. Nevertheless, as Aa/a ^ 1 we assume that the FLRW metric 
is a good approximation to the actual metric, and therefore that our 
lookback times are approximately correct. The ACDM parameters 
used are: (Ho, Qm, ^^a) = (70.5 kms"^ Mpc"\ 0.2736, 0.726). 

the fit of Aa/a for a combined VLT 



We show in figure 
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Aq /a = Br cos(O) + m. 



(19) 



+ Keck sample ("combined r-dipole sample") against r cos(O) — 
ct cos(0). The parameters for this fit are: B = 1.1 x 10~® GLyr~^ 
(la confidence limits [0.9, 1.3] x 10"'^ GLyr"!), RA = (17.5 ± 
1.0) hr, dec. = (-62±10)° andm = (-0.187±0.084) x 10"^ 
Using the bootstrap method we assess the statistical significance of 
this fit with respect to the monopole-only fit as 4.15(t. In figure[23] 
we show the confidence regions on the dipole location for the VLT, 
Keck and combined samples on the sky. 

In galactic coordinates, the pole of this fit is at approximately 
{I, b) — (330, —15)°. The fact that the pole and anti-pole are close 
to the Galactic Plane explains the relative lack of absorbers near to 
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Figure 19. Effect of removing quasar spectra on the statistical significance 
of tlie dipole, as assessed tlirough a jaclc-knife method. Each spectrum is re- 
moved one at a time, and the value of the statistical significance of the model 
Aa/« = A cos(0) + m is calculated with respect to the monopole model, 
using the bootstrap method. We use a Gaussian kernel density estimator 
to construct the approximate probability density function of the effect of 
quasar spectrum removal, where the width of the Gaussian basis functions 
has been chosen to be the inverse of the number of spectra. The cumulative 
distribution function is plotted as a dashed, red line. This demonstrates that 
the angular dipole effect is not due to a single spectrum. 
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Figure 21. Binned values of Aa/o plotted against 2'' cos 0, for /3 = 0.41. 
The 2*^ dipole+monopole model is preferred over the monopole-only model 
at the 3.9a" level. Importantly, this graph only covers cos(0)| < 1. 
Given that it is possible to probe up to redshift z < 4 with the MM method, 
judicious choice of observational targets close to the dipole axis might be 
able to extend this horizontal range of this graph up to ~ ±2, thereby po- 
tentially increasing sensitivity to the effect substantially if it is real. 
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Figure 20. Effect of removing quasar spectra on the statistical significance 
of the dipole, as assessed through a sampling method. For 100,000 sam- 
ples, 5 quasar spectra are randomly removed from the combined Keck -I- 
VLT sample, and the statistical significance of the dipole model Aa/o = 
A cos( 0) + m is calculated through the method of |Cooke & Lynden-Bell| 
Gaussian basis functions with width 1 / 316 (=10'' ''■^) are used. 
This graph demonstrates that, in the absence of particular knowledge about 
problematic spectra, the chance of obtaining a dipole model where the sta- 
tistical significance of the dipole is less than 3(t is small as a result of ran- 
domly removing 5 spectra. 

the pole and anti-pole in both the Keck and VLT samples, a fact 
made obvious in figure [T2|earlier. 

If we adopt a dipole-only model, 

Aa/a = Br cos{e), (20) 

we derive parameter estimates B = 1.1 x lO"'^ GLyr~^, RA — 
(17.5±0.9) iir, dec. = (-58±9)°. Icr confidence limits on B are 
[0.9, 1.3] X 10~^. The statistical significance of the dipole model is 
99.998 percent (4.22o-). The confidence limits on the dipole loca- 
tion for this fit for the VLT, Keck and combined samples are shown 
in figure [23] 



6 SYSTEMATIC ERRORS 

6.1 General effect of systematic errors on monopole and 
dipole terms 

[Murphy et al.| ( (2001b||2003| considered a number of potential sys- 
tematic effects, including: quasar absorber isotopic abudances for 
Mg that are different from terrestrial values; differential isotopic 
saturation; spectrograph temperature variations; the possibility that 
a single transition could be spuriously generating a non-zero Aa / a 
(which could occur if the laboratory wavelength for this transition 
was significantly inaccurate); atmospheric dispersion; stretching or 
compression of the quasar wavelength scale relative to the ThAr 
calibration scale; absorption cloud kinematics; systematic blends; 
random blends; and others. 

Generating a significant angular (dipolar) variation in a from 
systematic effects is, however, rather harder than producing an off- 
set from Aa/a — 0. Whatever this effect is - if it exists - it must 
be well correlated with sky position or must be a combination of 
systematics that by coincidence mimics angular variation in a. 

The most significant potential concern relates to inter- 
telescope wavelength calibration systematics, and we treat this is- 
sue in detail in section [63] However, in section [STS] we raised the 
presence of a consistent offset from Aa/a = in both the \ow-z 
Keck and VLT samples (the low-2: monopole), and suggested that 
this may be due to evolution in the abundance of the Mg isotopes. 
We deal with this concern in the next section. 



6.2 Effect of isotopic abundances 

Most of the atomic species we use have a number of stable iso- 
topes, and each of these isotopes exhibits a slightly different rest 
wavelength for a given transition. The isotopic spacing depends on 
the transition and species under consideration, but scales accord- 
ing to the inverse square of the mass. That is, Aui oc a;o/mf. 
The isotopic shifts of Mg, as the lightest of the elements species 
under consideration, are relatively significant. The terrestrial abun- 
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Figure 22. Binned values of A«/o plotted against r cos(0) = ct cos 0, where t is the look-back time to a redshift 2. ACDM parameters are from |Hinshaw| 
|et al.|(2009^ . The top panel (triangles) shows the combined sample, binned with approximately 25 absorbers per bin. The bottom panel shows VLT (circles) 
and Keck (squares) points, binned with approximately 12 absorbers per bin. The red (solid) line in both cases shows the model, Aa/a = i3rcos(0) + m. 
The parameters for the fit are: B = 1.1 X IQ-"^ GLyr"! (Icr confidence limits [0.9,1.3] X 10"" GLyr"!), RA = (17.5 ± 1.0) hr, dec. = (-62± 10)° 
and m = (—0.187 it 0.084) X 10~^. It is interesting that this simple model is a reasonable representation of the data. 



dance of the Mg isotopes is ^"Mg^^Mgr^Mg = 79:10:11 (Ros 



[man & Taylor|1998] >. We define the heavy isotope fraction as F 
(^''Mg + '"^Mg)/Mg, which has a terrestrial value of Tt = 0.21. 

We have assumed for our final fits that the quasar absorber iso- 
topic abundances are the same as the terrestrial abundances. How- 
ever, if the abundances in the absorbers differ from the terrestrial 
abundances, this will introduce a small but potentially significant 
shift in the quasar absorption lines compared to laboratory mea- 
surements. Mg will be most affected by this, due to its low atomic 
mass compared to the other species. Previous work l |Murphy et al.| 
|2003 ) noted that the effect could be particularly significant for low- 
z absorbers, as these predominantly consist of the Fe/Mg combina- 
tion. High-2 absorbers are less likely to be affected due to the use 
of more massive anchors (Si and Al), for which this effect is less 
relevant. Additionally, the use of many transitions with differing q 
coefficients at high redshift will tend to reduce the importance of 
this effect jMurphy et al.|2004| l. 

Both observations 1 Gay & La mbert|2000| l and theoretical es- 
timates ( T immes et al.| |l995 1 of stellar abundances for Mg sug- 
gest that the heavy isotope abundance of Mg (i.e. the '^^Mg and 



Mg isotopes) decreases with decreasing metallicity. Murphy et al 



(j2003j noted that the low-z Mg/Fe systems considered in the Keck 
sample have relative metal abundances, [Fe/H], in the range —2.5 
to 0.0, whereas the high-z DLA systems have relative metal abun- 
dances of about —1.0. Therefore, the quasar absorbers we con- 
sider may also have sub-solar metallicities. However, observations 
of some low metallicity red giants show significant enrichment of 
the heavy Mg isotopes. [Ashenfelter et al.| ( |2004a|b| l considered a 



"modest" enhancement of the stellar initial mass function (IMF) for 
intermediate mass stars (M « SM©), and showed that this could 
produce V ~ 0.4 for [Fe/H] ~ -1.5. |Fenner et al.|p005l ) argued 
that such an IMF would substantially overproduce nitrogen rela- 
tive to observations, and therefore that this mechanism of creating 
r > Ft does not seem possible. 

However, it appears that the link between stellar evolution 
and the likely nitrogen abundance in quasar absorbers is not fully 
understood. [Centurion et al.| j2003[ > described observations of ex- 
tremely low relative abundances of nitrogen in DLAs, and thus 
argued that nitrogen production cannot be dominated by massive 
stars. In a detailed study, [Dessauges-Zavadsky et al.| ([2007 ) argued 
that "no single star formation history explains the diverse sets of 
abundance patterns in DLAs". [Melendez & Cohen] ( |2007[ l claimed 
(in contrast to previous analyses) that heavy Mg isotope enrich- 
ment due to AGB stars in the Galaxy halo does not occur until 
[Fe/H] > — 1.5. [Levshakov et al.| ( |2009| examined II metal-rich, 
high-redshift (1.5 < z < 2.9) quasar absorbers and argued that 
the nitrogen abundance is uncorrelated with the metallicity, which 
implies that nitrogen enrichment has several sources. 

From the arguments above, it appears that the observational 
situation concerning F at high redshift is uncertain. There are no 
stringent, independent observations which constrain F in our sam- 
ple. We therefore treat F as unknown and explore what happens if 
we vary it. 

We first consider F < Ft . To place an upper limit on the effect 
of F < Ft, we re-fit all the VLT absorbers with no ^'^Mg or ^'^Mg, 
and similarly re-fit the absorbers in [Murphy et al.| ( [2004^ using no 
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Figure 23. Sky map in equatorial coordinates showing tlie 68.3 percent (Icr 
equivalent) confidence limits of the location of the pole of the dipole for 
a fit to the Keck Aa/a values (green region), VLT l\a/a values (blue 
region) and combined h.a/a values (red region), for a fit of Aa/a = 
i?rcos(©) + m (top figure) and A«/o = Brcos(©) {bottom figure), 
where r = ct and t is the lookback time to the absorber. The pole and anti- 
pole of the CMB dipole are marked as PcMB ™d Acmb respectively. 
In the model which includes a monopole (top figure), the Keck confidence 
region is large due to a relative degeneracy with the monopole; the region 
is much smaller in the bottom figure on account of no monopole term being 
included. 



^^Mg or ^""Mg. We give the parameters for the fits to the Keck, VLT 
and combined samples in this situation in table H] The confidence 
regions on the dipole location are shown in figure|24| Importantly, 
the dipole model remains statistically significant at the 3.5(j level. 
The reduction in significance from 4. Icr is primarily due to extra 
scatter introduced into the Aa/a values about the model. The ex- 
tra scatter implies that the F = model is not a good model for the 
absorbers. Additionally, the monopole becomes statistically signif- 
icant at the 5.7a level. Thus, a lower heavy isotope abundance in 
the quasar absorbers is unable to explain the dipole effect, and ad- 
ditionally increases the significance of the monopole term. The in- 
crease in significance of the monopole term mirrors the result in 
[Murphy et all l llOOSl l. 

We now consider the impact of increasing the Mg heavy iso- 
tope fraction (F > Ft). In section [F!8| we discussed the presence of 
a low-z monopole in both samples, where the difference between 
the two samples is remarkably small. Explaining this result via al- 
terations to the Mg isotope abundance would require enrichment 
of the heavy isotope fraction relative to terrestrial values. If we as- 
sume that all of the z < 1.6 monopole is due to relative enrichment 
of the heavy Mg isotopes, we can extrapolate from the F = Ft and 
r = cases to estimate {Fz<i.6) using a simple linear model. A 
linear model may be used as the response of Aa /a to changes in F 
is linear ( [Murphy et al.|2004l >. This model assumes that the ratio of 
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Figure 24. la confidence regions for the Keck (green), VLT (blue) and 
combined (red) dipoles in the circumstance where absorbers containing Mg 
are fitted with no ^^Mg or ^''Mg, to mimic the maximum possible effect of 
a lower heavy isotope fraction in the quasar absorbers compared to ten'es- 
trial values. Although the confidence regions are enlarged as a result of extra 
scatter introduced into the data, the reasonable alignment between the sam- 
ples is still maintained. The separation between the Keck and VLT dipole 
vectors increases to 32°, which has a chance probability of 11 percent. 

^^Mg/^^Mg is fixed. For z < 1.6, m = (-0.390 ± 0.108) x lO"'' 
for F = Ft, and m = (-0.884 ± 0.115) x 10"^ for the case 
F = 0. Under our linear model, {rz<i.6) ~ 0.32 in order to make 
m = 0. If we take a^n ~ 0.108 x 10~^ as a representative error, 
this yields (T.^i.e) = 0.32 ± 0.03. 

In summary, variations in the magnesium heavy isotope frac- 
tion have the potential to significantly impact the monopole com- 
ponent of the angular dipole -I- monopole model, but cannot explain 
angular variations in a. 

6.3 A test for systematic kinematic segregation of species 

In the absorption system modelling described earlier in this paper 
we have tied corresponding velocity components to have the same 
redshift. We are thus assuming that the velocity structure for each 
complex is the same in those species which are tied. The justifi- 
cation, apart from empirical, is that we restrict the assumption to 
species with similar ionisation potentials. If there were gross de- 
partures from this assumption we would be unable to derive sta- 
tistically acceptable fits in all cases. However, this is not found. 
Thus, if there are kinematic segregations between species, these 
must therefore be small compared to the parameter error estimates 
for individual redshifts. 

Kinematic segregation between different ions might in princi- 
ple result from several causes, including ionisation potential differ- 
ences, mass segregation due to gravitational effects, and a possible 
dependence on the formation site. Kinematic segregation clearly 
applies only to different species or potentially to different ionisation 
stages of the same element. Pairs of transitions in the same species 
are immune to the effect. Thus we formed two groups of pairs of 
transitions, intra- and inter-ion using a previously studied subset of 
our entire sample, that present in [Chand et al.[p004| l. We chose 
that sample for this experiment because it was initially selected on 
the basis of simple velocity structure, which usually means that all 
the transitions studied are not saturated. However, two of the ab- 
sorption systems from the [Chand et al.[ sample towards the same 
quasar (J000344-232355, Zabs = 2.1854, 2.1872) were not used 
because this complex failed the initial selection tests described in 
Section lT2l 
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Table 4. Effect of removal of ^^Mg and ^^Mg isotopes on the model i^a/a = Acos(S) + m. Generally speaking the effect is to push Aa/a to more 
negative values. The column "5 A" gives Icr confidence limits on A. The column labelled "significance" gives the significance of the dipole+monopole model 
with respect to the monopole model. This also introduces extra scatter into the Aa / a values about the dipole model, which implies (unsurprisingly) that fits 
with no heavy Mg isotopes are not a good representation of the absorbers. Despite the extra scatter, the dipole model is still significant at the 3.5(t level. 
Additionally, the monopole term here becomes significant at the 5.7(t level. iTrand is given for the different samples; HC refers to the Keck high-contrast 
sample. 



Sample 




A (10-^) 


5A(10-5) 


RA (hr) 


dec. (°) 


m(10-5) 


significance 


VLT 


1.04 


1.20 


[0.78, 1.75] 


18.1 ± 1.4 


-65 ± 14 


-0.439 ±0.197 


1.9o- 


Keck 


1.63 for HC 


0.42 


[0.32,0.88] 


16.6 ± 2.2 


-35 ±35 


-0.835 ±0.156 


0.4(7 


VLT + Keck 


As above 


0.98 


[0.77, 1.23] 


17.3 ± 1.0 


-59 ± 10 


-0.528 ±0.092 


3.5(7 



There is an additional advantage of using this particular sam- 
ple for a kinematic segregation test. In this test we measure veloc- 
ity shifts, Av, between different absorption transitions to identify 
kinematic segregation. A non-zero Aa/a may manifest itself as 
such a velocity shift. The redshift distribution and location on the 
sky of these absorption systems happens to be such that the values 
of Aa/a predicted by our dipole model are small. The velocity 
splitting for each pair is then 

Aa 

Av = vi - V2 ~ -2cAq , (21) 

a 

where Aq — 52/1^2 — Qi/i^i- Here the subscripts 1 and 2 refer 
to the different members of the pair. Column 7 in table |5] shows 
that the dipole-predicted values are small for these absorption sys- 
tems. Hence, for this sample, whether or not there is a cosmological 
dipole in a, these particular absorption systems provide a clean test 
for kinematic segregation. 

In order to measure the velocity shifts between pairs of 
transitions, we fitted all species simultaneously. Redshifts and b- 
parameters of corresponding velocity components were tied in the 
usual way. Aa/a was fixed at zero. All column densities of all 
velocity components were free to vary. However, in order to al- 
low complete freedom for one transition to shift relative to another, 
i.e. to allow us to identify any kinematic segregation, if present, for 
a system having m different absorbing transitions, m — 1 new free 
parameters were introduced into each fit, where each of the m — 1 
parameters allowed for a potential relative velocity shift for that 
transition relative to the mean system redshift. (There were m — 1 
velocity shifts to avoid degeneracy with redshift). 

In doing it this way, we can allow for, and measure most sen- 
sitively, bulk shifts between one ion and another, which, if present, 
could emulate a non-zero Aa/a for any particular absorption sys- 
tem. However, we note we would require an extremely unusual uni- 
verse if kinematic segregation were to emulate a spatial dipole in 
Aa/a since each gas cloud would need to be aligned with some 
particular direction in the universe. 

The m — 1 measured velocity shifts were paired in order to 
obtain a roughly equal number of inter- and intra-ion pairs. The way 
in which this was done depended on transitions available. However, 
the systematic pairing procedure we followed was to match pairs of 
transitions of different strengths, such that at least one member of 
the pair provided a good estimate of the global velocity structure. In 
order to derive a consistent velocity sign convention over all pairs, 
the velocity difference for each pair is in all cases measured relative 
to the bluest line in the pair. 

The individual velocity measurements are presented in Table 
|5] giving 34 inter-ion and 28 intra-ion pair measurements. In Fig- 
ure[25] we show the distributions for intra- and inter-ion velocities. 
A random kinematic segregation between inter-ion pairs would be 
revealed as a broadening of the inter-ion distribution relative to the 



intra-ion distribution. The F test is sensitive to such a broadening 
and we use it to test the equality of the variances of the intra- and 
inter-ion distributions. We calculate an F statistic of 1.91, which 
implies that the probability that the variances of the two distribu- 
tions differ is only ~9%. In the unlikely case that kinematic segre- 
gation is non-random the effect would be revealed as a shifting of 
the mean of the inter-ion distribution with respect to the intra-ion 
distribution. The Mann- Whitney U test is sensitive to such a shift. 
Comparing these 2 distributions the Mann- Whitney U statistic is 
403, which for a two-sided test implies that the probability that the 
two distributions differ is ~30%. We thus detect no significant dif- 
ference between the inter- and intra-ion distributions. 

In summary, we have devised a method which checks for kine- 
matic segregation and which is similar to the way in which we have 
measured Aa/a. Although we have tied species together, we al- 
low for bulk velocity shifts between species by introducing the ad- 
ditional velocity shift parameter. It is exactly this kind of bulk kine- 
matic segregation which, if present, could emulate a variation in a. 
An a spatial dipole may only be produced by kinematic segrega- 
tion if and only if the kinematic segregation is non-random. Our 
Mann- Whitney U test results show no evidence for non-random 
segregation. If random kinematic segregation were present in the 
sample used for this test (our F test), that would be indicated by 
an increase in arand (see Section [3.5.3| earlier in this paper). At 
present there is no evidence for any kinematic segregation amongst 
the atomic species used in the current sample and thus we find no 
evidence that the dipole signal we have found is a result of kine- 
matic segregation. 



6.4 UVES, a dual-arm spectrograph 

UVES is a dual-arm spectrograph, where the incoming light is split 
into a red arm and a blue arm using a dichroic mirror. In princi- 
ple, misalignment of the slits in the blue arm relative to the red arm 
would produce a distortion of the wavelength scale between the 
two arms, which could mimic a change in a if transitions were fit- 
ted simultaneously from the blue and red arms. [Molaro et al.1p008l > 
investigated the possibility that such misalignment might cause ve- 
locity shifts between the blue and red arms, using measurements of 
asteroids, and argue that the two arms do not show separation by 
more than ~ 30 ms~^ in the situation where the science exposures 
are bracketed by the ThAr exposures. 

However, we note that |Molaro et al.| used a slit with of 0.5", 
which is rather different to the ~ 0.7" to --^ 1.0" typical of the 
quasar exposures. The UVES archive indicates that, for the obser- 
vations of [Molaro et al.[ the seeing was always poorer than the slit 
width. If the slits for the blue and red arms are misaligned, one 
would expect the induced effect on wavelength calibration to de- 
pend on slit size. In the seeing-limited regime, the slit is relatively 
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Table 5. Results for Av, the velocity separations between inter-ion and intra-ion pairs. The emission redshift of the quasar and the absorption redshift are 
given by Zem and Zabs respectively. Columns five and six show A-u and the corresponding Icr error for each pair. We can predict the expected for each 
pair given a model of a variation. Colvmm seven shows the predicted value of Ad corresponding to Aa/a = ylcos(©) + m. 



Quasar Name 




■^abs 


Species 


/\v (kms~^) 


<TAi;(kins ^) 


predicted (kms~^) 


J134427- 103541 


2.35 


0.873 


Fell2344/Mgll2796 


-0.446 


0.557 


0.031 








Mgll2803/Mgl2852 


-1.808 


0.744 


0.001 








Fell2374/Fell2383 


0.553 


0.867 


0.003 






1.278 


Fe II 2344/Mg II 2796 


-0.446 


0.557 


0.031 








Mg II 2803/Mg 1 2852 


1.171 


0.992 


0.001 








Fe II 2374/Fe II 2600 


0.034 


0.950 


0.003 








Mn II 2576/Mn II 2606 


-1.476 


3.759 


0.012 






1.915 


Sill 1526/Fe II 2600 


0.104 


0.268 


-0.041 








Fell2344/Fell2383 


0.135 


0.390 


-0.004 


J012417-374423 


1.91 


0.822 


Fell2600/Mgll2796 


0.237 


0.325 


0.017 








Mg II 2803/Mg 1 2852 


-0.792 


0.401 


0.000 








Fell2344/Fell2383 


-0.189 


0.883 


-0.002 






0.859 


Fe II 2600/Mg II 2796 


0.252 


0.473 


0.017 








Mg II 2803/Mg 1 2852 


0.561 


0.958 


0.000 








Fell2344/Fell2383 


0.191 


0.708 


-0.002 






1.243 


Fe II 2374/Mg II 2796 


-0.413 


0.409 


0.019 








Fe II 2260/Mg 1 2852 


0.545 


0.790 


0.017 








Fell2344/Fell2383 


-0.124 


0.542 


-0.002 








Aim 1854/Al III 1862 


-0.065 


0.532 


0.002 


J024008-230915 


2.23 


1.635 


Mg II 2803/Mg 1 2852 


0.584 


0.478 


-0.004 








Siiil526/Mgii2803 


0.328 


0.252 


0.003 








Fe II 2344/Fe II 2587 


-0.164 


0.471 


0.009 








Aim 1854/Al III 1862 


0.292 


0.480 


-0.005 






1.637 


Fe II 2344/Mg II 2796 


-0.126 


0.292 


-0.033 








Mg II 28()3/Mg 1 2852 


-1.347 


0.712 


-0.001 








Fell2587/Fell 2600 


0.662 


0.453 


-0.004 






1.657 


Allll670/Mgll2803 


0.131 


0.244 


-0.006 








Sill 1526/Mgl2852 


0.332 


0.780 


0.002 








Fell 2344/Fe II 2383 


-0.211 


0.451 


0.005 


J000344-232355 


2.28 


0.452 


Fell 2600/Mg II 2796 


0.001 


0.142 


0.005 








Mgll 2803/Mg 1 2852 


0.163 


0.659 


0.000 








Fell 2344/Fe II 2383 


0.436 


0.280 


-0.001 


JOl 1143-350300 


2.40 


1.182 


Fe II 2600/Mg II 2796 


-0.038 


0.297 


0.013 








Mg II 2803/Mg 1 2852 


-0.229 


0.788 


0.000 








Fe II 2344/Fe II 2383 


-0.424 


0.653 


-0.002 






1.350 


Fell 2600/Mg II 2796 


0.006 


0.090 


0.013 








Mg II 2803/Mg 1 2852 


0.252 


0.164 


0.000 








Fell 2344/Fe II 2383 


-0.134 


0.134 


-0.002 


J222006-280323 


2.41 


0.942 


Fell 2600/Mg II 2796 


-0.015 


0.225 


0.047 








Fell2587/Mgll2803 


0.357 


0.278 


0.056 








Fell 2344/Fe II 2383 


0.312 


0.529 


-0.006 






1.556 


Allll670/Mgll2796 


-0.285 


0.224 


0.004 








Fe II 2383/Fe II 2600 


0.135 


0.170 


0.000 








Fe II 2344/Fe II 2587 


0.344 


0.458 


-0.011 


J135038-251216 


2.61 


1.439 


Fell 1608/Mg II 2796 


0.622 


0.802 


-0.055 








Siiil808/Mgii2803 


0.919 


0.824 


0.014 








Znll2026/Mnll2594 


0.335 


1.272 


0.052 








Crll2056/Crll2066 


-0.802 


1.319 


0.012 








Fe II 2260/Fe II 2344 


0.163 


1.153 


0.000 


J045523-421617 


2.67 


0.908 


Mg II 2803/Mg 1 2852 


0.234 


0.293 


0.000 








Fe II 2344/Fe II 2600 


0.150 


0.175 


0.000 






1.858 


AllI1670/MgII 2796 


-0.507 


0.818 


0.000 








Fe II 2383/Fe II 2600 


-0.030 


1.220 


0.000 


J000448-415728 


2.76 


1.542 


Fe II 2383/Fe II 2600 


-0.102 


0.896 


0.000 








Fe II 2344/Fe II 2374 


-0.199 


1.125 


-0.010 






2.168 


Silll526/Allll670 


0.171 


0.199 


-0.012 








Mgll2796/Mgll2803 


0.002 


0.091 


0.004 








Fell 1608/Fe II 2344 


-0.753 


0.624 


-0.078 






2.302 


Si II 1526/Fe II 1608 


-0.205 


0.117 


0.030 








Nilll709/Nilll741 


0.142 


0.540 


0.037 


J212912-153841 


3.28 


2.022 


Sill 1808/Mg II 2803 


0.533 


0.401 


-0.001 








Fell 2344/Fe II 2587 


-0.395 


0.578 


-0.017 
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Figure 25. Histograms of velocity differences between intra-ion (solid line) 
and inter-ion pairs (dashed line). A Mann-Whitney U test and the F test find 
no significant differences between these 2 distributions (see text). Thus, we 
have not found any evidence for kinematic segregation between the species 
used. 



uniformly illuminated, and therefore the observed science wave- 
lengths should be well calibrated through the ThAr exposure. On 
the other hand, when the seeing is much better than the slit, one 
might expect to see larger differences, if such differences exist. 



6.5 Inter-telescope systematics & the A« test 

We have previously argued that the observed angular variation in a 
is unlikely to be caused by systematic effects. Any such systematic 
must be well correlated with sky position in order to explain the 
good alignment between dipole directions in dipole models fitted 
to the Keck and VLT samples. Nevertheless, accurate wavelength 
calibration for quasar spectra remains an issue of importance given 
the limitations of the existing ThAr calibration process. Therefore, 
in this section we directly explore wavelength-dependent system- 
atics using common observations of quasars from both Keck and 
VLT. 

One can ask whether there is some systematic effect which is 
constant in time between the two telescopes which might gener- 
ate the observed dipole effect. All of the Keck spectra used in the 
analysis in this paper were acquired whilst HIRES had only one 
CCD chip. In this configuration, multiple exposures are needed to 
yield full wavelength coverage. If the quasar image is not precisely 
centred in the spectrograph slit for every exposure, velocity offsets 
between spectral segments obtained at different times are possible. 
This issue should be substantially mitigated at VLT, as UVES can 
acquire almost the entire spectral range in a single observation. The 
effect could be exaccerbated in conditions of good seeing and could 
include an additional small effect due to the seeing profile decreas- 
ing slightly towards the red end of the spectrum. 

It turns out that the VLT and Keck samples have 7 quasars in 
common. We give a list of the quasars common to the VLT and 
Keck samples in table |6] The use of common sources allows one 
to search for problems with wavelength calibration; absorption fea- 
tures should be found at the same barycentric vacuum wavelength 
between different exposures. Note that the number of absorption 
lines which can be used for this purpose is much larger than is used 
for analysing Aa/a. Whilst for Aq/q many absorption lines are 
needed to yield a single measurement of Aa/a, in principle each 



Table 6. List of quasars common to the VLT and Keck samples. Keck names 
are given by their B 1950 designation as they appear in Murp hy et al.,^2003) , 
whereas VLT names are given by their J2000 designation. 



Keck sample name VLT sample name 



0216+0803 


J021857+081727 


0237-2321 


J024008-230915 


0940-1050 


J094253- 110426 


1202-0725 


J120523-074232 


0528-250 


J053007-250329 


1337+1121 


J 134002+ 110630 


2206-1958 


J220852- 194359 



absorption line in the spectrum yields one constraint on potential 
wavelength distortion. 

To explore potential wavelength scale distortions, we propose 
a method which we refer to as the An test. The method proceeds 
as follows: i) for each common quasar, visually identify regions of 
non-terrestrial absorption, typically having width of a few A; ii) for 
each of these regions, perform a Voigt profile fit to the VLT spectral 
data (identification of the transition responsible is unimportant); Hi) 
fit corresponding spectral regions of the Keck and VLT simultane- 
ously, but with an extra free parameter. An, which allows for a ve- 
locity shift between the two spectral regions. The VLT spectral data 
for these regions were fitted using an automated Voigt profile fitting 
routine designed to fit regions of the forest automatically. An is de- 
fined hereafter as the velocity difference An = n (VLT) —n (Keck) 
which must be applied to minimise between two comparable 
spectral regions. Each value of An provides an estimate of the 
velocity offset between the two telescopes at that observed wave- 
length, giving An(A). One can therefore examine the functional 
form of An(A)i, where i refers to the ith quasar pair under consid- 
eration. For each set of An values from a particular spectral pair, 
we use the LTS method to calculate the weighted mean of that set 
of An values, which we then subtract from the An values for that 
spectral pair. This is to remove any constant offset resulting from 
mis-centering of the quasar within the slit. We use k = 0.95n for 
the LTS fit (see section|333]l. 

Any relative wavelength scale distortion can in principle be 
removed by applying an inverse function based on the observed 
An data. To see this, consider the form of the distortion. For an 
absorption line with rest wavelength Ao, observed wavelength Ai, 
and velocity distortion An then 

A, = Ao(l + 2)(l + An/c), (22) 

where we have assumed that An is constant over the absorption 
profile under consideration. The effect of Aa/a can be ignored - 
whatever transition is being examined is the same in both spec- 
tra, and so any effect due to a change in constants will be ab- 
sorbed into the determination of z. There are two options to at- 
tempt to remove the wavelength scale distortion given some func- 
tion An(Aobs)- One could modify the spectral data, changing the 
observed wavelengths as Aobs — Aoba/[l + An(Aobs)/c)]. When 
one fits a particular transition, the other possibility is to perturb the 
rest wavelength of the transition fitted, as 

Ao^Ao[l + An(Aobs)/c]. (23) 

We use the second option for ease of implementation within VPFIT. 
Doing this means that the value of Aa / a derived from the fit will 
be the same as if the wavelength scale from the other telescope 
in the spectral pair had been used, thereby removing any inter- 
telescope differences (provided that An is correctly specified). 
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Figure 26. Binned values of Av points for six quasar spectral pairs ("core pairs"), with vertical range restricted for clarity. Some data points lie outside the 
visible range due to the contribution of gross outliers; these points are automatically ignored in our robust fit. 5 values of Av contribute to each bin. Av is 
defined as the velocity difference which must be applied to the Keck spectrum to minimise for ^ joint fit. Visually, it is not clear that a common trend 
exists for the six quasars shown. For 02 16+0803/J02 1857+08 1727, At; appears to increase with increasing wavelength. For 0237— 2321/J024008— 230915 and 
1202-0725/J120523-074232, it is not clear whether a trend exists. 0940- 1050/J094253- 110426 shows a non-linear effect. 0528-250/J053007-250329 
shows a clear trend for a decrease in Av with wavelength. The effect in 1337+1 121/J134002— 1 10630 is unclear given the limited wavelength coverage and 
relative lack of data. 



We show the Av data for 6 of the quasar spectral pairs ("core 
pairs"), which appear similar to each other, in figure|26] We analyse 
the Av data from these quasars in the following section. We noticed 
a problem with the 7th pair, 2206-1958/J220852-194359, which 
displays variations of Av with wavelength which are grossly dif- 
ferent from the other six pairs. A systematic trend in Av is seen, 
with a maximum difference in A v of ~ 2.5 km over the range 
4000 < A < 6OOOA. In section 6.5.3 we apply an inverse func- 



tion derived from the Av data seen in this spectral pair to the all the 
VLT spectra, and show that a distortion of this type cannot affect all 
the data. We consider the joint impact of the Av functions from the 
6 core quasars and from 2206- 1958/J220852- 194359 in section 
EM 

In all our analysis in this section we have removed those ab- 
sorbers which were previously flagged as outliers from consider- 
ation in the statistical analysis (the z = 1.542 absorber towards 
J000448-415728 in the VLT sample and the 2 ^ 2.84 absorber 
towards Q1946+7658 in the Keck sample). 



6.5.1 Av data & 6 core quasar pairs 

We note the presence of significant outliers within the set of Av 
values from these 6 spectral pairs. Therefore we rely wholly on 



robust statistical methods to estimate parameters for phenomeno- 
logical models of Aw (A). 

We discuss the particular trends for each spectral pair in the 
caption for figure |26] but note that no common trend is seen be- 
tween these spectra, as would be expected if no systematic wave- 
length distortions exist. Importantly, the functional form of Av ap- 
pears to differ in both magnitude and sign between quasars. This 
suggests that any relative wavelength distortion is likely to average 
out over a large number of absorbers. Additionally, the wavelength 
coverage of the Av data for most spectra is significantly smaller 
than the wavelength range within which MM absorbers are fitted. 
This means that from each spectral pair it is impossible to tell what 
the wavelength distortion might be over large amounts of the spec- 
tral range. 

Due to the fact that the Aw values from each spectral pair do 
not densely span the whole spectroscopic wavelength range, we 
combine the Av values from each of the six core pairs together in 
order to estimate a common function which spans the full wave- 
length range. The functional form of this is unknown, however a 
high-order polynomial cannot be statistically supported. We use a 
linear function as a first approximation. We fit the linear function 
with the LTS method, using k = 0.95n. We show this linear fit in 
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Figure 27. LTS linear fit to the Ati data for the 6 core quasar pairs. The 
dashed blue lines show the la confidence limit on the linear fit (red line). 
We have utilised k = 0.95n for this fit to obtain robustness against some 
outliers. The data shown are binned with 5 points per bin, but only those 
points contributing to the fit have been included. For a fit of Au = a\ + b, 
a = (-7±14) X lO-^kms-i k'^ andb = 0.38±0.71 kms^i. Most 
transitions used in the MM analysis fall in the range 4000 < A < 7000A; 
from this graph there is no significant evidence for a significant wavelength 
distortion in this region. 



figure [27] For the form 

/\v = a\ + b, (24) 

a = (-7±14) X lO^'^'kms-i and 6 = 0.38 ± 0.71 kms'^ 
Note firstly that a is statistically consistent with zero. Therefore, it 
is difficult to conclude that a common linear systematic exists in 
the A?; data. Nevertheless, in the next section we apply an inverse 
function of this type to the VLT spectral data to determine the ef- 
fect that a distortion of this type and magnitude would have on the 
detected dipole effect. 

A legitimate question to ask is whether the choice of fc = 
0.95n is reasonable. We show the effect of difference choices of k 
in figure |28] Our estimate of the slope is not overly sensitive to a 
choice of k. With the exception of a small region around k — 0.7n, 
the general trend is for the slope to decrease with decreasing k. The 
fact that the slope decreases with increased trimming implies that 
the underlying trend may be less than what we have estimated. 



6.5.2 Application to VLT sample 

To investigate the effect of the potential wavelength distortions 
from the 6 core pairs, we apply an inverse Av function (equation 
\24\ to all the VLT absorbers by perturbing the rest wavelengths of 
the transitions fitted in each absorber, as described in section [631 
We apply the same linear function in every VLT spectrum fitted. 
This therefore puts the VLT and Keck spectral data on a common 
wavelength scale. Any observed angular variation in a which sur- 
vives the inverse function can not be due to stable inter-telescope 
wavelength calibration differences. 

Because we apply the inverse function by perturbing the rest 
wavelengths of transitions in absorbers fitted, we can only do this 
where each fitted transition occurs in only one spectral region in a 
particular fit. There are two pairs of absorbers where we have fit- 
ted both absorbers in the pair simultaneously, because a transition 
from one absorber in the pair overlaps with a transition from the 
other absorber (at a different redshift) in the pair. In this case, a 
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Figure 28. Slope of At; vs A under the LTS method with different amounts 
of data excluded from the fit. The red (solid) line shows a weighted LTS 
fit, and the blue (dashed) line shows the unweighted fit. Importantly, the 
magnitude of the slope decreases as fc — > 0.5n, implying that the slope we 
utilise may be an over-estimate of any underlying effect. 

particular transition can be fitted twice (in two absorbers, in two 
widely separated spectral regions). However, we cannot apply two 
different perturbations to a single rest wavelength. Therefore, we 
remove these two pairs of absorbers to form a "VLT reference 
set". Thus, we compare a VLT set of Aa/a values where the An 
inverse function has been applied with the Aa/a values from a 
VLT reference set. The two pairs of absorbers which are removed 
are the z ~ 2.253 and z ~ 2.380 absorbers associated with 
I2I4225-4420I8 and the z ~ 1.154 and z ~ 0.987 absorbers in 
the same spectrum (i.e. 4 absorbers are removed to form the VLT 
reference set). This means that the VLT reference set contains 149 
absorbers. 

In table]?] we give the results of applying the Av inverse func- 
tion above to those absorbers in the VLT reference set. The effect 
is generally to push Aa/a to more negative values. We show an 
updated plot of the confidence regions of the Keck, VLT and com- 
bined dipole locations in figure ]30] Although the statistical signif- 
icance of the dipole-l-monopole model (Aa/a = y4cos(0) + m) 
over the monopole-only model decreases from 3.9o- (reference set) 
to 3.1(7, the position of the VLT (and therefore combined) dipole is 
effectively unchanged. This accords well with our earlier argument 
that because the detection of a dipole is a differential effect, it is dif- 
ficult to emulate through any simple systematic. The Keck and VLT 
dipoles in this case are separated by 25°, which has a chance prob- 
ability of 7 percent (see section ]53] ). Also note that introducing this 
modification to the wavelength scale of the VLT spectra does not 
significantly change the good alignment between the z < 1.6 and 
z > 1.6 samples. The dipole directions in this case are separated 
by 13°, which has a chance probability of 2 percent. 

We would expect that any good model for a wavelength- 
dependent systematic should quantitatively improve the fit of the 
dipole model to the Aa/a values. To see if the Av model sig- 
nificantly improves the fit, we compare the AICC of the model 
Aa/a — ylcosS + m fitted to different sets of Aa/a val- 
ues: i) AICCvlt.Ad, the AICC of the angular dipole model fit- 
ted to the VLT absorbers in the VLT reference set with the lin- 
ear Av function described above applied; ii) AICCvLT.rot, the 
AICC of the angular dipole model fitted to the VLT absorbers 
in the VLT reference set; Hi) AICCvLT,Au+Kcck, the AICC for 
the angular dipole model fitted to the absorbers in set (i) com- 
bined with the Keck absorbers in the Keck04-dipole set, and; iv) 
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AICCvLT,rof+Keck, the AICC for the angular dipole model fitted 
to the absorbers in set (ii) combined with the Keck absorbers in 
the Keck04-dipole set. In all cases we apply the same values of 
Crand to the VLT Acv/a values in order to compare points on a 
like-with-like basis. We find that AICCvlt.Au — AICCvLT,rcf ~ 
—0.8, indicating that the set of VLT absorbers with the in- 
verse function applied is preferred, but not significantly. Compar- 
ing the VLT-FKeck set to the equivalent reference set, we find that 
AICCvLT,Au+Kcck — AICCvLT.rcf+Kcck ~ —3.2, indicating that 
the VLT-l-Keck set where the Aw inverse function has been applied 
to the VLT absorbers is weakly preferred. However, when compar- 
ing the reference sets and the Au sets, the AICC does not account 
for the extra two parameters for the linear model of An vs A. Thus, 
with a two-parameter model for the Au function there is no sig- 
nificant preference for the At; results, and thus there is no strong 
evidence in the Aa/a values themselves for a wavelength distor- 
tion of this type. 

In deriving the results above, we have assumed that Au values 
from different spectral pairs may be legitimately combined in order 
to estimate a common systematic. This may not be a good assump- 
tion, given the differences in the signal-to-noise of the spectral data, 
the spectral range which the Au values cover, the potential func- 
tional form of Ai;(A)i and number of exposures. We then proceed 
as follows: i) fit a linear model to the Aw values from each spectral 
pair using the LTS method; ii) from each model, estimate Aw (A) 
along with an uncertainty on the estimate; Hi) for the six estimates 
of At; at each A, form a weighted mean of the estimates, At;,i,(A), 
and calculate the associated uncertainty, and; iv) plot At;u,(A) as a 
function of wavelength. We show the result of this in figure[29] Im- 
portantly, under this model we can find no wavelength where At; is 
statistically different from zero. 

In summary: we are unable to detect a statistically significant 
linear wavelength distortion common to the 6 core spectral pairs. 
Applying to the entire VLT spectral sample a simple linear model 
for At; (A) from the six core pairs reduces the statistical signifi- 
cance of the dipole, but the statistical significance still remains high 
enough to be of interest. The systematic applied here does not de- 
stroy the good alignment between the fitted Keck and VLT dipole 
vectors. We are therefore unable to remove the dipole effect from 
the combined Keck and VLT sample. 



6.5.3 2206- 1958/J220852~ 194359 

We show the At; values from the pair in figure[3T| The Au test only 
examines calibration differences between Keck and VLT, and so we 
cannot tell whether Keck or VLT is responsible for significant trend 
in Au for this spectral pair. 

The limited spectral range (~ 4000 < A < 6000 A) of the At; 
data means that we simply do not know what the functional form of 
At; is for this spectral pair at A > 6OOOA. In order to estimate the 
potential impact of the distortion present on Aa/a values in the 
whole sample, we need knowledge of At;(A) at all observed wave- 
lengths. One possibility is to assume that the relationship is linear, 
but then the extrapolation over the whole spectral range results in 
a total change in At; of ~ 5kms~'^, which is comparable to the 
velocity width of the spectrograph slit; this seems too extreme. Ad- 
ditionally, the At; values do not appear by-eye to be linearly related 
to A; the gradient of At; with A visually appears steepest in the mid- 
dle of the At; values, and flatter near the edges of the At; data. We 
therefore try an arctangent model. 
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Figure 29. At;™ (A), a joint estimate of At; (A) from the 6 core quasar pairs 
made without combining the data into a single linear fit. The red (solid) line 
shows the estimate for Ad, and the dashed (blue) lines show the 95 per- 
cent confidence interval on the estimate. One can see that over the range 
where most of our At; data are obtained (4000A < A < 7000A) that 
At) is relatively flat. Ati diverges from zero for A < 35OOA (not shown) 
and for A > 7000A, however there are few Ad measurements to con- 
strain Ad in these regions. At no wavelength is Ad significantly different 
from zero. Because the linear model used for each quasar fit is unlikely to 
be a true description of any underlying wavelength distortion, there is also 
uncertainty due to model specification, which is naturally not included in 
the confidence region shown above. As such, the confidence region shown 
is under-estimated. These considerations show that there is no statistically 
significant evidence for a common systematic from consideration of the 6 
core quasar pairs. 
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Figure 30. Icr confidence regions for the Keck (green), VLT (blue) and 
combined (red) dipoles, where the linear inverse function derived from the 
Ad data from the 6 core spectral pairs has been applied to the VLT set. 
Although the statistical significance of the dipole decreases from 3.9a" to 
3. Icr, the positions of the VLT and combined dipole are effectively un- 
changed. 



Applying the LTS method to this fit, with k = 0.95n, yields: A — 
(-0.98 ± 1.03) kms"^A-\ k = (1.9 ±2.6) x 10"^ A"\ A^ = 
4547±526 A and & = 0.48±0.74 kms-\ where errors are derived 
from the diagonal terms of the covariance matrix at the best fit. The 
uncertainty on each At; values has been increased in quadrature 
with 0.26 km s^^ to account for overdispersion about the LTS fit. 

The impact of this wavelength perturbation on the values of 
Aa/a when this function is applied to the VLT absorbers is se- 
vere. In particular, large numbers of points are scattered away from 
Aa/a = 0, inducing a highly significant detection of Aa/a at 
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Table 7. Results of applying the inverse /S.v function from figure|27]to the VLT absorbers under the model A« /a. = A cos © + m. The column "5 A" gives 
1(T confidence limits on A. The column labelled "sig" gives the significance of the dipole+monopole model with respect to the monopole model. The origin 
of the reference VLT set is described in section [6.5.1| The impact of the (non-significant) linear distortion modelled from the 6 core quasar pairs overall is a 
reduction in the statistical significance of the dipole from 3.9(t to S.lcr. 



Sample 


Crand (10 




A (10-5) 


SAilQ-^) 


RA(hr) 


dec 




) 


m (10" 




sig 


VLT reference 


0.88 




1.21 


[0.80, 1.72] 


18.3 ± 1.1 


-61 


± 


13 


-0.110 ± 


0.179 


2.2ct 


VLT with Ad function applied 


0.95 




1.02 


[0.65, 1.52] 


18.6 ± 1.3 


-61 


± 


16 


-0.262 ± 


0.183 


l.Scr 


VLT reference + Keck 


Keck HC= 


1.63 


0.97 


[0.57, 1.39] 


17.4 ± 1.0 


-61 


± 


10 


-0.177 ± 


0.085 


3.9o- 


VLT with Av function + Keck 


Keck HC= 


1.63 


0.79 


[0.59, 1.07] 


17.4 ± 1.2 


-60 


± 


12 


-0.273 ± 


0.085 


3.1o- 




4000 



6000 



8000 



Wavelength (A) 



Figure 31. Au values for the 2206- 1958/J220852- 194359 spectral pair. 
The At) values demonstrate a substantial distortion, with |(5(Ad)| 
2.5kms~^ over ~ 2000A. The trend is substantially larger, and different 
to the trends seen in the six core pairs (see figure[26). Note that although dif- 
ferent trends are seen individually in each of the six core pairs, no common 
trend is seen. 



z ~ 0.8 and z ~ 1.5. For instance, a formal weighted mean of all 
the 1.3 < 2 < 1.8 points yields Aa/a = (-3.96 ± 0.12) x lO"'^ 
- a 33(T "detection". If one multiplies the error by \/xI^ to ac- 
count for = 14.7 about the weighted mean, then one obtains 
Aq/q = (-3.96 ± 0.46) x lQ-^ still an S.Sa "detection". Such 
a signal is seen in neither the Keck or VLT samples, which im- 
mediately implies that this particular relative wavelength distortion 
can not possibly apply to the all of either the VLT or Keck spectra. 
That is, the distortions seen in this particular KeckA'LT spectral 
pair appear not to be representative of a significant fraction of the 
entire sample. However, the fact that we have identified this distor- 
tion demonstrates the power and utility of the quasar pair analysis 
in identifying systematic errors, even when their actual origin re- 
mains unknown. 



6.5.4 Overall effect of wavelength systematics using the Au test 

We now investigate whether a diluted form of the above effect (i.e. 
the effect from the 2206- 1958/J220852- 194359 pair) could ex- 
ist in the spectral data in combination with the (non-significant and 
much smaller) effect observed from the 6 core quasars. To do this, 
we use a Monte Carlo approach where we apply the inverse func- 
tion derived in section [6!5. 1 | (equation|24[l to 6/7 of the quasar spec- 
tra selected at random in the VLT sample, and the arctan function of 
equation 25 to the remaining 1/7 of the spectra. We then apply the 
LTS method to estimate a new (Jrand- We then add the Keck sample 



tistical significance of the dipole using the bootstrap method. The 
mode of the distribution obtained is ~ 2.2a, with quite substantial 
variation between iterations. 

To determine whether the Au function significantly improves 
the goodness-of-fit in the VLT sample, we compare the AICC at 
each iteration in the Monte Carlo simulation for a dipole model 
(Aa/a = AcosQ + m) fitted to the VLT Act/oi values in that 
iteration with the AICC from a dipole model fitted to the Aa/a 
values in the VLT reference set, where in each iteration we use 
frand ~ 0.88 X 10^^ iu Order to compare Aq/q values on a like- 
with-like basis. We show this distribution in figure[32] In only 3.5 
percent of iterations is the AICC lower than in the reference set. 
This implies that it is unlikely that a wavelength distortion of this 
type is present in our data set. However, in almost all of the itera- 
tions the AICC is much larger than the AICC from the VLT refer- 
ence set; the median AAICC = 43.5. Importantly, in no case is the 
AAICC > 10. Thus in no case can we say that there is very strong 
evidence in favour of the model with the Aw function applied. Ad- 
ditionally, the AICC does not account for the 6 parameters used in 
deriving the Ai; model - we should expect a significant reduction 
in the AICC if the A?; function is a good model. From this argu- 
ment, we thus conclude that a wavelength distortion of this type is 
unlikely to be present in the VLT spectral data. 



to this new VLT sample. At each iteration, we calculate the sta- 



6.6 Comment on'Griest et al. 72010^ I2 and ThAr 
measurements on Keck/HIRES 

The path that the quasar light takes through the telescope is simi- 
lar but not identical to that of the thorium-argon (ThAr) calibration 
lamp. In particular, the ThAr lamp nearly uniformly illuminates the 
slit, whilst the quasar image will be centrally concentrated, particu- 
larly if the seeing is smaller than the slit width used. On account of 
the path differences, it is possible for distortions in the wavelength 
scale to emerge. Such distortions may be long range or short range. 

|Griest et al. ' compared the wavelength calibration scale of 
exposures derived from I2 absorption spectra and the standard 
thorium-argon (ThAr) calibration exposures. The distortions are 
identified by comparing the calibration of the spectrum derived 
from a ThAr lamp to the calibration obtained when an I2 absorption 
cell is used. The iodine cell is placed in the quasar light path, and 
therefore the wavelength scale obtained for the quasar light using 
the iodine cell suffers from no path differences. They reported drifts 
between the I2 and ThAr calibration scales of up to 2000 ms~^ 
over several nights, and claimed that that "this level of systematic 
uncertainty may make it difficult to use Keck HIRES data to con- 
strain the change in the fine-structure constant". 

The Au test above explicitly includes the effect of any drifts in 
the wavelength calibration both within a single night and between 
observation nights. From figure vrn the RMS of the 66 binned 
points about the fit is ~ 250 ms" . The mean wavelength sep- 
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Figure 32. AICC for 1000 iterations of a Monte Carlo simulation 
where the Ad function from the 6 core pairs (equation \1A) is applied 
to 6/7 of the VLT quasars at random, and the Av function from the 
2206- 1958/J220852- 194359 pair (equation|25) is applied to the remain- 
ing 1/7 of the quasars. The AICC here is calculated for the Aa/a data with 
respect to the model /\a/a = y4cos(0) + m for the VLT points only. 
The dashed red line shows the actual AICC for the VLT reference set (see 
section [6.5.2) . Only in 3.5 percent of the iterations is the AICC lower than 
for the reference set, which suggests that it is unlikely that a distortion of 
this type is present in the VLT data. The median value of the AICC here is 
191.4, corresponding to a median AAICC = 43.5. For the median case, 
the odds against a wavelength distortion of this type being present in the 
data are ^ 3 X lO" : 1. 



aration between these points is comparable to the echelle order 
width. This RMS can therefore be compared directly to the spread 
in Wsiiift seen in figure 5 of |Griest et al.| In contrast to their spread of 
~ 2000 ms^^, we see see typical wavelength distortions between 
VLT and Keck which are some 8 times smaller. We have directly 
quantified the impact of this on our measurements of Aa/a in sec- 
tion [63]2] Our results demonstrate that it is possible to reliably 
use Keck/HIRES data to constrain the fine-structure constant from 
quasar observations. In the next section, we deal with the effect of 
'intra-order distortions' reported by |Griest et al.| 



6.7 Intra-order wavelength distortions 

A pattern of wavelength distortions within each echelle order has 
been identified both within Keck/HIRES spectra l |Griestetal.|2010^ 
and VLT/UVES spectra ( jWhitmore et al.|2010| l. In this section we 
attempt to estimate the impact of the extra scatter that has already 
been introduced into the VLT Aa/a values as a result of the intra- 
order distortions reported by |Whitmore et al!] 

Griest et al. (2010) identified an apparently repeating pattern 
of distortions that occur within each echelle order ('intra-order dis- 
tortions'), where the wavelength scale of pixels near the centre of 
the order is displaced with respect to the calibration scale for pix- 
els at the echelle order edges. The peak-to-peak amplitude of the 
distortion is ~ 500 ms~^ at ~ 5600A. As the distribution of MM 
transitions is random with respect to the location of the echelle or- 
ders, the effect of these distortions will be random from absorber 
to absorber. A distortion of this type, with no monotonic long- 
range component, constitutes a random effect (section |33l.|MuF 



phy et al.jj2009[ l applied a model of the distortion found by [Griest 
et al.| ( |2010| l to the 2004 Keck results, and found that the impact on 
the weighted mean of the Aa/a values was effectively negligible. 
It is also worth noting that systems which utilise a large number of 
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Figure 33. An function used for the investigating the wavelength distor- 
tions found by Whitmo re et ah] based on a Fourier analysis of their data. 
This function is repeated to longer and shorter wavelengths. 



MM transitions will be less sensitive to an effect of this type. This 
is because with many transitions, the distortion is sampled in many 
locations; if the distortion does not have a long-range component, 
the average distortion must tend to zero. 

[Whitmore et al.| ( [2010| l have identified similar intra-order dis- 
tortions in VLT/UVES spectra. Whilst the distortions found by |Gri- 
lest et al. 'display a repeating pattern, the distortions found by |Whit- 
more et al. are much more irregular, displaying little consistent pat- 
tern. It is important to note that because the VLT spectra used here 
are the result of the co-addition of many exposures, taken with dif- 
ferent echelle grating settings and over many nights, it is expected 
that any distortions of the wavelength scale due to light path differ- 
ences should reduce in magnitude. Thus, we consider the possible 
estimate of the impact of the effect we present to be an upper limit. 

To construct a model for the [Whitmore et aT] distortions, we 
carried out a Fourier analysis of the velocity shift data presented 
in that paper to produce a Au model which we can apply to the 
VLT absorbers. The iodine cell absorption lines used to establish 
the intra-order distortion results only occur over the wavelength 
range ~ 5000-6200A. We are therefore forced to assume that our 
model of these distortions applies to much bluer and redder wave- 
lengths as well. Clearly, the important part of the model is the am- 
plitude rather than the period of the distortions; our model has a 
maximum peak-to-peak distortion of ~ 300 ms^^. We show this 
model in figure [33] 

In table J8] we show the result of applying the function shown 
in figure]33]to the VLT absorbers using equation|23] The impact on 
the location of the dipole and the value of the monopole is minimal, 
as expected. However, we note that the Urand required is somewhat 
larger, which means that this model of the wavelength distortion has 
introduced extra scatter into the Aa/a values. Any good model of 
the systematic should reduce, not increase, the scatter. The extra 
scatter reduces the significance of the dipole, but does not destroy 
the good alignment between Keck and VLT, nor between low and 
high redshift samples. In particular, the chance probability of align- 
ment for the Keck and VLT samples (where the VLT sample has 
been altered with this A?; model is 6 percent, the chance probabil- 
ity of alignment between low and high redshift samples is 4 percent, 
and the joint chance probability for these two factors is 0.3 percent. 

The presence of intra-order wavelength distortions would 
serve to increase the scatter of the Aa/a values about the true 
values. These distortions can only randomise but not bias Aa/a 
values. They can not manufacture a dipole or monopole. Were we 
able to make the same quasar observations without the presence 
of any wavelength scale distortions, the scatter in the Aa/a values 



Spatial variation in the fine -structure constant 35 



Table 8. Results of applying the inverse Ati function from figure[33]to the VLT absorbers. The column "5 A" gives Icr confidence limits on A. The column 
labelled "sig" gives the significance of the dipole model with respect to the monopole model. The origin of the reference VLT set is described in section [6.5.1| 
For 7 = 2, we remove the absorber at z = 1.6574 toward J024008— 230915, as it is identified as an outlier with the LTS method. Note that the estimates of 
the dipole location and monopole value do not differ greatly between samples 2 and 3. For errand. HC refers to the Keck high-contrast sample. Application of 
the Whitmore et. al. Av function causes a reduction in significance of the dipole model of Ki O.Scr. This reflects the maximal amount by which the significance 
has already been reduced; the 3.3a- figure is not the significance after applying a correction for a systematic. 
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Sample 


O-rand (10 ^) 


A(10-^) 


(5A(10-^) 


RA (hr) 


dec. (°) 


m (10-^) 


sig 


1 


VLT with Whitmore Ati function z 


1.090 


0.87 


[0.51, 1.43] 


19.0 ± 1.6 


-54 ± 23 


-0.025 ±0.215 


1.2o- 


2 


#1 with z = 1.6574 absorber removed 


1.067 


1.03 


[0.63, 1.56] 


18.4 ± 1.1 


-51 ± 19 


-0.090 ± 0.204 


1.6o- 


3 


VLT reference 


0.882 


1.21 


[0.80, 1.72] 


18.3 ± 1.1 


-61 ± 13 


-0.110 ± 0.179 


2.2cr 


4 


#2 + Keck sample 


Keck HC=1.630 


0.88 


[0.68, 1.12] 


17.4 ± 1.0 


-56 ± 11 


-0.217 ± 0.090 


3.3o- 


5 


VLT reference + Keck 


Keck HC=1.630 


0.97 


[0.57,1.39] 


17.4 ± 1.0 


-61 ± 10 


-0.177 ± 0.085 


3.9o- 



about the model should be smaller (and so (Trand would be smaller). 
We would therefore expect that this would increase the significance 
of the dipole model. Our analysis in this section suggests that the 
maximal reduction in statistical significance of the dipole which 
may have occurred as a result of intra-order wavelength distortions 
present is ~ 0.6a. 



7 DISCUSSION 

7.1 Summary of results 

We have presented the results of a new analysis of 154 quasar ab- 
sorbers over the range 0.2 < z < 3.7 using the many-multiplet 
method, derived from VLT/UVES spectra of 60 quasars. These ab- 
sorbers allow us to constrain changes in a at often better than the 
10"^ level over much of the observable universe. 

Assuming that the VLT Aa/a values are described by a sim- 
ple weighted mean, we found that {Aa/a)^ = (0.21 ± 0.12) x 
10~^. This appears at first to be inconsistent with the results of 
[Murphy et aL] l |2004| , who reported ( Aa/a)^ = (-0.57±0.11) x 
lO"'' from 143 absorbers observed with Keck/HIRES. 

We fitted a simple angular dipole-l-monopole model, Aa/a — 
AcosQ + m, to the VLT Aa/a values and found some (2.2ct) 
evidence for angular (and therefore spatial) variations in a, with 
the dipole direction (i.e. the direction of maximal deviation in a 
from the current laboratory value) towards RA — (18.3 ± 1.2) hr, 
dec. = (-62 ± 13)°, with an amplitude A = l.lStoll x 10"^ 
and a monopole term m = (-0.109 ± 0.180) x 10"^ A dipole- 
only model (Aa/a — AcosO) points in the direction RA — 
(18.4 ± 1.3) hr, dec. = (-58 ± 15)°, with an amplitude A = 
0.99lo.2g X 10"^ 

A similar angular dipole-l-monopole model applied to the Keck 
Aa/a values gives model parameters RA = (16.0 ± 2.7) hr, 
dec. = (-47 ± 29)°, with an amplitude A = l.lSto'tl x 10"^ 
and a monopole term m — (—0.465 ± 0.145) x 10"''. The signif- 
icance of the dipole-l-monopole model over the monopole model is 
low in this case, just O.Scr, and is due to the relatively low sky cover- 
age of the Keck/HIRES quasar lines of sight, leading to a stronger 
degeneracy of the dipole and monopole terms. For a dipole-only 
model we found model parameters RA = (16.4 ± 1.2) hr, dec. — 
(-56 ± 12)°, with an amplitude A = l.OGtoH x 10"^ This 
dipole-only model has a slightly higher significance (I.Ict, or 72 
percent), and is consistent with the direction and amplitude found 
from a dipole-only model fitted to the VLT Aa/a values. 

For the combined Keck-l-VLT sample, the dipole-l-monopole 
model (Aa/a — ^cos0 + m) is preferred over the monopole- 
only model at the 4.06a level, with the dipole pointing in the di- 
rection RA = (17.3 ± 1.0) hr, dec. = (-61 ± 10)°, with an 



angular amplitude A — 0.97^0 20 x 10 ^ and monopole term of 
m = (-0.178 ± 0.084) x 10"^ Fitting a dipole-only model in- 
creases the significance to 4.15a with the dipole pointing in a very 
similar one to the best-fit dipole-l-monopole model. Thus the com- 
bined VLT-l-Keck sample shows significant evidence for angular 
(and therefore spatial) variations in a. 

It is important to realise that, while the weighted mean Aa/a 
values from the VLT and Keck samples disagree, the two sam- 
ples are consistent when potential spatial variations in a across 
the sky are considered. The numerical results given above show 
that dipole models fitted to the VLT and Keck Aa/a values point 
in very similar directions. This good alignment would be surpris- 
ing if the dipole effect were not real. Under the model Aa/a — 
A cos O + m, the chance probability of obtaining as close an align- 
ment as observed (or closer) between the dipole vectors in both 
samples is 6 percent. Additionally, dipole vectors fitted to low- 
redshift (z < 1.6) and high-redshift (z > 1.6) subsamples point 
in a very similar direction, with a chance probability of 2 percent 
under the model Aa /a = A cos 6 + m. A joint analysis for both 
of these conditions finds a chance probability of 0.1 percent. 

In section [F!8] we examined the presence of a statistically sig- 
nificant monopole (3.6a) in the z < 1.6 subsample of absorbers 
in the combined sample. There is consistent evidence that this off- 
set of Aa/a from zero exists in both the VLT and Keck samples 
independently. This implies not only that the low-z monopole may 
be a real effect (i.e. not generated by intstrumental systematic er- 
rors), but that it appears not to be responsible for the observed 
angular variations in a. We suggested that evolution in the rela- 
tive abundances of Mg isotopes may be responsible for the low-« 
monopole. After examining this possibility in more detail (section 
\6.2\ , we suggested that an average heavy Mg isotope fraction of 
(TzKLe) ~ 0.32 ± 0.03 could account for the observed \ow-z 
monopole. 

To further check the robustness of the apparent spatial vari- 
ation in a, in section |5.9| we investigated the effect of itera- 
tively clipping potentially outlying Aa/a values from a combined 
Keck-l-VLT sample in which the raw statistical errors on Aa/a 
were retained. When outliers are clipped iteratively from the best- 
fit dipole-l-monopole model at each iteration until xi ~ the sig- 
nificance of the model (over a monopole-only model) increases to 
almost 7a. Importantly, its significance does not decline rapidly, 
even after removal of more than ~20 percent of the most outlying 
Aa/a values, suggesting that the dipole effect is robust. Similarly, 
the dipole direction is stable under removal of potentially outlying 
Aa/a values. 

In section [5.1 1| we applied simple phenomenological models 
to attempt to account for the fact that the amplitude of the dipole 
appears larger at higher redshifts. Under a model where the ampli- 
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tude is proportional to we find that /? = 0.46 ± 0.49, and that 
the dipole+monopole model is preferred over the monopole-only 
model at the 3.9(t confidence level. We also investigated a model 
where the amplitude grows in proportion to the lookback-time dis- 
tance, r = ct, and find that the dipole+monopole model for the 
combined Keck+VLT sample is significant over the monopole-only 
model at the 4.15(t level. Removing the monopole from consider- 
ation increases the significance to 4.22cr. The fitted dipole direc- 
tions are very similar to that found in the strictly angular models 
(i.e. with no redshift dependence). Under the r-dipole+monopole 
model, the amplitude of the dipole is (1.1 ± 0.2) x 10"*' GLyr"^ 
with a monopole term of m = (-0.187 ± 0.084) x 10~^. 

In section |6] we investigated various potential systematic er- 
rors. We showed that temporal evolution in the abundance of Mg 
isotopes cannot be responsible for the observed angular variation 
in a. We demonstrated a powerful quasar pair test, where by com- 
paring spectra of common quasars taken with both the VLT and 
Keck telescopes we can constrain relative wavelength scale dis- 
tortions between the two telescopes. Considering 6 quasar pairs, 
we could not detect a statistically significant common linear wave- 
length distortion. Applying an estimate of the (non-significant) dis- 
tortion observed, the significance of the dipole-l-monopole model, 
Aq/q = ^cos(O) + m, over the monopole-only model was re- 
duced to 3.1ct. The sign and functional form of the wavelength 
scale distortion differs in each of the quasar pairs, so any resid- 
ual systematic in the full sample may be smaller than what was 
found. We also considered a seventh quasar pair which shows a 
very large relative wavelength scale distortion (which is quite un- 
like that seen in the other 6 pairs). We argued that it is unlikely that 
a distortion of this magnitude is found in a significant fraction of 
either the Keck or the VLT absorbers because of the large increase 
in the scatter of the Aa/a values it would imply, a scatter grossly 
inconsistent with that actually observed. In section [677 ] we showed 
that the intra-order wavelength distortions reported by fCriest et al.| 
p0T0| l for Keck/HIRES and |WhitmoreetaL]pOTOt for VLT/UVES 
do not plausibly explain the observed dipole effect. 

As we discussed in section [53| any systematic which mimics 
angular variation in a must not only be well correlated with sky 
position, but must do so in a way which is consistent between the 
two telescopes. A systematic which is correlated with zenith angle 
is not sufficient to produce the observed dipole effect; such an ef- 
fect should produce a variation in a that is approximately symmet- 
ric about the latitudes of the telescopes projected onto the sky (i.e. 
dec. ~ 20° for Keck and dec. ~ —25° for VLT), which is not what 
is seen. Importantly, such an effect is unable to produce the consis- 
tency observed between the dipole locations. The fact that dipoles 
fitted to low-2; and high-2; subsamples of the Aa/a values (com- 
bined from the two telescopes) also point in the same direction on 
the sky makes a simple systematic which explains the results even 
more difficult to identify. The low-z and high-2: absorbers utilise 
different transitions, which have a different pattern of line shifts as a 
function of Aa/a. Thus, on average, the low-2: and high-z subsam- 
ples will respond differently to simple systematic errors, making it 
difficult to explain consistent low-z and high-z dipole directions. 
We know of no systematic at present which is able to explain the 
observed dipole effect. 

In summary, the results presented here provide significant and 
seemingly robust evidence for angular variations in a across the 
sky. We find evidence that the dipole amplitude grows with increas- 
ing distance. Taken together, these findings suggest a simple picture 
where there is a gradient in the value of a through the Universe, al- 



though this is strictly an observationally motivated interpretation 
and may be just the (seemingly) simplest many possible others. 

7.2 Implications 

The dipole-like variation in a presented here, if confirmed, would 
be a detection of new physics at the most fundamental level. It 
would directly demonstrate the existence of a preferred frame in 
the universe; it may be that this anisotropy could be detected in 
other cosmological measurements. Additionally, it would demon- 
strate that the Einstein Equivalence Principle is only an approxima- 
tion. 

Some consider that the universe is fine-tuned for life (see for 
instance [Davies|2003^ , in that the values of the fundamental con- 
stants appear to be set in such a way that small variations (of order a 
few percent) in some fundamental constants would seemingly make 
some ingredients for life as we understand it (e.g. requiring wa- 
ter) unstable. If the fundamental constants vary throughout space, 
this has the potential to resolve this possible fine-tuning problem: 
the universe need not be globally fine-tuned for life. Instead, Earth 
may simply be located in a region of space where the constants are 
amenable to life. This may imply the existence of regions where 
the constants take on values where life as we understand it is not 
possible. 

7.3 Relationship with other experiments 

[Berengut & Flambaum|pblO^ explicitly demonstrate that the re- 
sults presented here are consistent with atomic clock measure- 
ments, the natural nuclear reactor at Oklo, and measurements of 
/?-decay in meteorites. [Berengut et al.| ( |2010) considered other po- 
tential cosmological dipoles (for example, in the primordial deu- 
terium abundance and other dimensionless constants). To the best 
of our knowledge, the results presented here are not in conflict with 
the results of any experiment. 



8 CONCLUSIONS 

Our work significantly extends the work of [Murphy et al.H2003| 
|2004[ l from Keck/HIRES by adding a data set of comparable size 
using VLT/UVES. 

We find that the combination of both data sets yields a de- 
tection of spatial variation of ct. That is, we observe variation in 
a across the sky and with an increasing amplitude at greater dis- 
tances from Earth. There are several lines of argument which sug- 
gest that the dipolar effect observed may be real: ;) The mild ev- 
idence (~ 2cr) for angular variations in a within the VLT sample 
(which cannot be the result of inter-telescope systematics). ii) The 
close alignment between dipole vectors in both dipole+monopole 
and dipole-only models fitted to the VLT and Keck samples sepa- 
rately. In the dipole+monopole case (Aa/a = A cos(0) +m), this 
has a chance probability of 6 percent. Hi) The close alignment be- 
tween dipole+monopole models fitted to the 2 < 1.6 and 2 > 1.6 
cuts separately, which has a chance probability of 2 percent, iv) The 
robustness of the result, both in terms of statistical significance and 
stability of dipole direction, to the removal of significant numbers 
of absorbers, v) The general difficulty of generating angular varia- 
tions in a through simple systematic effects. 

If the observed result is due to systematics, it must be a com- 
bination of finely-tuned systematics, the nature of which is un- 
known. We have used spectra of a small number of quasars, all ob- 
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served with both the VLT and Keck, to investigate inter-telescope 
wavelength-dependent systematics. We conclude that the (gener- 
ally insignificant) relative wavelength distortions observed are un- 
able to explain the dipole effect observed. The origin of the wave- 
length calibration issues detected at both Keck and VLT is not well 
understood, and as such the analysis presented here reflects our cur- 
rent state of knowledge. Clearly, the power of comparing spectra of 
the same quasars from different telescopes for ruling out - or de- 
tecting - relevant wavelength distortions, from any source, even 
unknown ones, should be exploited in a much larger sample. 

The apparent dipole-like variation in a presented here is man- 
ifestly verifiable and falsifiable. Due to the orientation of the dipole 
relative to the observations, there is a relative lack of absorbers near 
its pole and anti-pole. Future observations toward the pole and anti- 
pole will have increased sensitivity to detect a significantly positive 
or negative Aa/a, respectively, if the effect is real. Furthermore, 
we have demonstrated the power of detecting systematics, even of 
unknown origin, using the quasar pair analysis; a campaign to ob- 
serve quasars near the dipole pole and/or anti-pole would best be 
made with two telescopes targeting the same quasars so that any 
confirmation of any dipole-like variation in a is maximally com- 
pelling and reliable. Because of the orientation of the dipole, these 
two telescopes would have to be in the same hemisphere if one 
wanted to apply a similar quasar pair analysis to that demonstrated 
here. 

We noted the presence of a statistically significant monopole 
in the z < 1.6 sample. Evidence for this effect is present in both the 
VLT and Keck samples. Because the effect is common to both sam- 
ples, it cannot be the cause of the angular variations in a observed. 
Evolution in the abundance of the heavy Mg isotope fraction may 
be the most fikely cause for this, but further work is needed to re- 
solve the issue. 
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APPENDIX A: VLT Aa/a VALUES 

We give the Aa/o values from our MM analysis of the 154 VLT ab- 
sorbers in table [aT] An ASCII version of this table is available at http : , 
I //astronomy . swin . edu . au7 ~mmurphy /pub . html| 



version of this table, including newly measured transitions, which should 
be used for future analyses is available from the authors on request. 



fits for the VLT many-multiplet systems analysed. Each absorber is plotted 
on a velocity scale, such that con'esponding components align vertically. 
Velocities are given as differences from an arbitrary redshift, which is cho- 
sen to be close to the maximum optical depth of the absorber. The positions 
of fitted components are indicated by blue tick marks. Plotted above each fit 
are the residuals of the fit, that is [fit-data]/en'or, where the error is the Icr 
uncertainty associated with each flux pixel. The two red lines indicate it lo", 
within which the residuals are expected to occur about 68% of the time if 
the errors are Gaussian, the error array is correct and the fitted model is a 
good representation of the data. 

Each plot contains a maximum of 16 regions. In the event that there 
are more fitting regions than this, the fit is split into several parts. Each part 
may contain common transitions so as to provide a common reference, and 
to illustrate the velocity structure more clearly. 

The full set of Voigt profile fits for all absorbers can be found at 
[http : / /astronomy . swin . edu . au/ -mmurphy/pub . htmlj 

This paper has been typeset from a T^f/ ETbX file prepared by the 
author. 



APPENDIX B: ATOMIC DATA 




APPENDIX C: VOIGT PROFILE FITS 
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(N 


C^ 




(N 


O) 


(N 




ID 


ID 


ID 


ID 


ID 


ID 


ID 








5^ 


5^ 


00 


00 


OO 


00 


C") 


m 


m 


ID 


ID 


ID 


ID 


ID 


00 






cS 


cS 


cn 


cs 


cS 




cS 


cS 


cS 


cS 


oi 


c4 


(N 


c^i 


(N 


CM 


oi 








(N 














cn 




(N 


(N 


CNl 


r4 


(N 


cn! 



o 


ID 


ID 


ID 


00 


00 


00 


00 


00 


00 


00 


ID 


ON 


On 


On 


ON 


On 


On 


ON 


m 


CNl 


Ol 


CM 
















ID 


ID 


ID 


ID 


ID 


ID 


ID 


ID 


o 


m 


CS 


m 


o 


o 


o 


o 


o 


o 


o 


NO 


cn 


C^ 


CS 


m 


C-) 


CI 


m 


-t 








CNl 


04 


04 


04 


04 


CM 


04 


CI 


-t 


-t 




-t 


-t 


-t 


-t 


NO 
-t 


1 1 


i 




5 5 5 5 5 
1 1 1 1 1 


5 


o 


ON 


On 


On 


ON 


On 


On 


ON 


1 


1 

On 


1 

ON 


1 

On 


>D 


1 

>D 


1 

ID 


1 

>D 


1 

>D 


1 

ID 


1 

>D 


1 

-t 


1 

<N1 


1 

04 


1 

Ol 


1 

CN 


1 

04 


1 

CN 


1 




>D 


ID 


ID 


04 


C^l 


04 


04 


O^l 


04 


04 


m 


ID 


ID 


ID 


ID 


ID 


ID 


ID 










04 


Ol 


04 


04 


Ol 


04 


04 




OO 


00 


00 


OO 


00 


00 


OO 


cn 








^ 


^ 


^ 


^ 


^ 


^ 


^ 


o 


O 


o 


o 


o 


o 


o 


o 
























Ol 


04 


04 


OnI 


04 


04 


Ol 


04 




04 


CnI 


CM 


04 


OnI 


Ol 


04 


OnI 


04 


04 


OnI 


04 


04 


OnI 


04 


04 


OnI 


04 



C-J c^ 

04 04 

CI m 

o o 

00 00 

OnI 04 

I I 

NO NO 

o o 

o 

04 04 

01 04 
OnI 04 



C-) c-1 

04 CM 

C-) C-J 

O O 

00 00 

04 OnI 

I I 

NO NO 

o o 

o o 

04 CM 

04 Ol 

04 OnI 



Ol 04 

o 
C-) 

-t -t 

04 04 

01 04 
I I 

NO NO 

ID ID 

r- 

04 04 

04 04 

04 04 



04 Ol 

o 

C-J 

-t -t 

01 04 

OnI 04 

I I 

NO NO 

ID ID 

04 04 

Ol 04 

OnI 04 



Ol 04 04 

00 00 00 

o o o 

On On ON 

O O O 

I I I 

NO NO NO 

5 5 5 

CO CO CO 

CO CO CO 

<N (N 



c-j 


C-J 


00 


00 


OO 


-t 


-t 


ID 


ID 


ID 


r~ 


1^ 


00 


00 


00 


-t 


-t 








04 


Ol 


04 
















+ 


+ 


+ 


+ 


+ 


ID 


ID 


00 


00 


00 


04 


OnI 


04 


04 


O^l 


NO 


NO 


NO 


NO 


NO 


^ 


^ 


^ 




^ 


m 


c<-j 








04 


OnI 


04 


04 


OnI 



Spatial variation in the fine -structure constant 43 



Table Bl. Atomic data for transitions usable in many-multiplet or alkali-doublet analyses, i.e. transitions with precise laboratory wavelengths. An updated 
version of this table, including newly measured transitions, which should be used for future analyses, is available from the authors on request. Information for 
isotopic and hyperfine components is given in italics. Columns 1 and 2 show the common names used for the transitions. Column 3 shows the mass number 
for each ionic species. The derivation of the laboratory wavenumbers, u)q, is summarized by the value of X as follows: - Measured wavenumber; 1 - 
Infen'ed from measured component wavenumbers; 2 - Inferred from measured composite wavenumber and measured component splitting; 3 - Infen'ed from 
measured composite wavenumber and calculated component splitting. Column 6 gives the reference(s) for the wavenumber measurement and/or calculations 
(specified below the table). Vacuum laboratory wavelengths, Aq, are derived from the wavenumbers. Columns 8 and 9 show the lower and upper/excited state 
electronic configurations. The ID letters in column 10 offer a simple shorthand for labelling transitions used to fit absorption systems. Column 11 s hows the 
ionization potential for the relevant ion, TP"*", and for the ion with a unit lower charge, IP~ . Column 12 shows the oscillator strengths, /, taken from |Morton| 
j2003| l or the relative strengths of the hyperfine or isotopic components. The latter are taken from |Rosman & Taylorlirgg SV The q coefficients and their 
uncertainties are from Dzuba et al. 1 1999b, a 2001 2002) and Berengut et al. 1 2004a b). Note that uncertainties in the q coefficients are representative, not 
statistical. Wavenumbers are on the Wh aling et al.|^l995^ Aril calibration scale; the Fell A1608/1611 and Nill wavenumbers have been scaled from their 
original values to account for the calibration difference between the Aril scales of Norlen^t,1973j and [Whaling et al^jl995). The exceptio ns to this are the 
Mg I/II wavenumbers which are on a highly accurate absolute scale generated using a frequency-comb calibration system. The |Whaling et al.H1995^ scale best 
agrees with this absolute scale. 



Ion 


Tran. 


A 


uJa rem "'"1 


X 


Ref. 


Ao [A] 


Lower state 


Upper state 


ID 


IP^ , 1P+ [eV] 


/ or % 


q [cm ^] 


Mg I 


2026 


24.31 


49346.772611(36) 


1 




2026.4749792(15) 


3s2 iSo 


3s4p ^PJ 




— ,7.65 


0.113 


87(7) 






26 


49346.854173(40) 





a 


2026.4716298(16) 










11.0 








25 


49346.807724(40) 





a 


2026.4735372(16) 










10.0 








24 


49346.756809(35) 





a 


2026.4756281(14) 










79.0 






2852 


24.31 


35051.28076(19) 


1 




2852.962797(15) 




3s3p ipj 


CL2 




1.83 


86(10) 






26 


35051.32015(25) 





b 


2852.959591(20) 










11.0 








25 


35051.29784(25) 





b 


2852.961407(20) 










10.0 








24 


35051.27311(17) 





b 


2852.963420(14) 










79.0 




Mgll 


2796 


24.31 


35760.85409(20) 


1 




2796.353794(16) 


3s 


3p ^P3/2 


bl 


7.65, 15.04 


0.6155 


211(10) 






26 


35760.940387(5) 





c 


2796.3470457(4) 






11.0 








25 


35760.92474(64) 


3 


c 


2796.348269(50) 


F = 2 


F = 1, 2, 3 






4.2 








25 


35760.86700(64) 


3 


c 


2796.352784(50) 


F = 3 


F = 2,3,4 






5.8 








24 


35760.837397(5) 





c 


2796.3550990(4) 










79.0 






2803 


24.31 


35669.30439(20) 


1 




2803.530983(16) 




3p ^Pi/2 


b2 




0.3058 


120(2) 






26 


35669.390571(5) 





c 


2803.5242094(4) 










11.0 








25 


35669.37651(64) 


3 


c 


2803.525314(50) 


F = 2 


F = 1,2,3 






4.2 








25 


35669.31684(64) 


3 


c 


2803.530004(50) 


F = 3 


F = 2,3,4 






5.8 








24 


35669.287670(5) 





c 


2803.5322972(4) 










79.0 




AllI 


1670 


26.98 


59851.976(4) 





d 


1670.78861(11) 


3s2 iSo 


3s3p iPi 


Cl 


5.99, 18.83 


1.74 


270(30) 


Aim 


1854 


26.98 


53916.554(1) 


1 


d 


1854.717941(34) 


3s2Si/2 


3p2P3/2 


di 


18.83,28.45 


0.559 


464(30) 






27 


53916.8149(8) 





d 


1854.708966(28) 


F = 2 






41.7 








27 


53916.3574(6) 





d 


1854.724704(21) 


F = 3 








58.3 






1862 


26.98 


53682.884(2) 


1 


d 


1862.791127(69) 


3s2Si/2 


3p2Pi/2 


d2 




0.278 


216(30) 






27 


53683.1953(15) 





d 


1862.780325(52) 


F = 2 






41.7 








27 


53682.6692(12) 





d 


1862.798581(42) 


F = 3 








58.3 




Sin 


1526 


28.09 


65500.4538(7) 





d 


1526.706980(16) 


3s23p2p°^2 


3s24s2Si/2 


ei 


8.15, 16.35 


0.133 


50(30) 






iO 


65500.441994 


3 


e 


1526.7072550 








3.1 








29 


65500.448002 


3 


e 


1526.7071150 










4.7 








28 


65500.454492 


3 


e 


1526.7069637 










92.2 






1808 


28.09 


55309.3404(4) 





d 


1808.012883(13) 




3s3p2 2D3/2 


62 




0.00208 


520(30) 






30 


55309.435938 


3 


f 


1808.0097601 








3.1 








29 


55309.387116 


3 


f 


1808.0113560 










4.7 








28 


55309.334806 


3 


f 


1808.0130660 










92.2 




Si IV 


1393 


28.09 


71748.355(2) 





d 


1393.760177(39) 


2p63s2Si/2 


2p63p 2P3/2 


/i 


33.49, 45.14 


0.513 


862(20) 






30 


71748.551629 


3 


e 


1393.7563579 






3.1 








29 


71748.451219 


3 


e 


1393.7583084 










4.7 








28 


71748.343484 


3 


e 


1393.7604012 










92.2 






1402 


28.09 


71287.376(2) 





d 


1402.772912(39) 




2p63p2Pi/2 


/2 




0.254 


346(20) 






30 


71287.574290 


3 


e 


1402.7690098 










3.1 








29 


71287.473031 


3 


e 


1402.7710024 










4.7 








28 


71287.364387 


3 


e 


1402.7731402 










92.2 




Till 


3067 


47.87 


32602.627(2) 





9 


3067.23750(19) 




3d24p z4D°^2 


91 


6.82, 13.58 


0.0489 


791(50) 






50 


32602.651577 


3 


h 


3067.2351837 








5.2 








49 


32602.640059 


3 


h 


3067.2362673 










5.4 








48 


32602.628061 


3 


h 


3067.2373961 










73.7 








47 


32602.603236 


3 


h 


3067.2397316 










7.4 








46 


32602.615933 


3 


h 


3067.2385371 










8.3 






3073 


47.87 


32532.355(1) 





9 


3073.86293(9) 




3d24p z-^DJ^j 


92 




0.121 


677(50) 






50 


32532.379612 


3 


h 


3073.8606027 








5.2 








49 


32532.368077 


3 


h 


3073.8616926 










5.4 








48 


32532.356062 


3 


h 


3073.8628278 










73.7 
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Table Bl - continued. Atomic data for transitions usable in many-multiplet or alkali-doublet analyses 



Ion 


Tran. 


A 


ujQ [cm ^] 


X 


Ref. 


Ao[A] 


Lower state 


Upper state 


ID 


IP-,IP+ [eV] 


/or% 


q [cm ^] 


Till 


3230 


Al.%1 


30958.586(1) 





9 


3230.12169(10) 




3d^4pz^Fg/2 


93 




0.0687 


673(50) 






50 
49 


30958.610542 
30958.599041 


3 
3 


h 
h 


3230.1191252 
3230.1203251 








5.2 
5.4 








48 


30958.587059 


3 


h 


3230.1215753 










73.7 








47 


30958.562268 


3 


h 


3230.1241619 










7.4 








46 


30958.574948 


3 


h 


3230.1228389 










8.3 






3242 


41.S1 


30836.426(1) 





9 


3242.91797(11) 




3d24p z4F°/2 


94 




0.232 


541(50) 






50 


30836.450997 


3 


h 


3242.9153410 








5.2 








49 

48 


30836.439283 
30836.427080 


3 
3 


h 
h 


3242.9165729 
3242.9178562 










5.4 
73.7 












3 
J 


I. 
ii 




















46 


30836.414/40 


3 


h 


3242. 91 91540 










8.3 






33o4 


47.87 


Zyj44.4j4(l ) 





9 


I^QA TIAOI CI 1 ^ 

J J54. /.?UUI(1 1 ) 




3d 4pz G°^2 


95 




U.3jo 


090(00 j 






50 


29544.480532 


3 


h 


3384.7269676 








5.2 








49 
48 


29544.468409 
29544.455781 


3 
3 


h 
h 


3384.7283564 
3384.7298032 










5.4 
73.7 








47 


29544.429586 


3 


h 


3384.7328042 










7.4 








46 


29544.442984 


3 


h 


3384.7312692 










8.3 




Cm 


2056 


52.00 


48632.058(2) 





9 


2056.256801(85) 




3d*4p 6p°^2 


hi 


6.77, 16.50 


0.103 


-1110(150) 




2062 


52.00 


48491.057(2) 





9 


2062.235929(85) 




3d*4p6pg/2 






0.0759 


-1280(150) 




2066 


52.00 


48398.871(2) 





9 


2066.163899(85) 






hs 




0.0512 


-1360(150) 


Mnll 


2576 


54.94 


38806.689(3) 





9 


2576.87534(20) 


3d54s a^Ss 


3d54p z^'P" 


h 


7.44, 15.64 


0.361 


1420(150) 






55 


38806.974333 


3 


i 


2576.8563955 


F = 0.5, 1.5 


F = 1.5,2.5 






14.3 








55 


38806.879265 


3 


i 


2576.8627082 


F = 2.5 


F = 1.5,2.5,3.5 






14.3 








55 


38806.768508 


3 


i 


2576.8700627 


F = 3.5 


F = 2.5,3.5,4.5 






19.0 








55 


38806.625155 


3 


i 


2576.8795818 


F = 4.5 


F = 3.5,4.5,5.5 






23.8 








55 


38806.451511 


3 


i 


2576.8911123 


F = 5.5 


F = 4.5,5.5,6.5 






28.6 






2594 


54.94 


38543.121(3) 





9 


2594.49669(20) 




3d54p z''P| 


12 




0.280 


1148(150) 






55 


38543.399993 


3 




2594.4778464 


F = 0.5, 1.5 


F = 0.5, 1.5,2.5 






14.2 








55 


38543.306507 


3 




2594.4841392 


F = 2.5 


F = 1.5,2.5,3.5 






14.3 








55 


38543.198206 


3 




2594.4914294 


F = 3.5 


F = 2.5,3.5,4.5 






19.1 








55 


38543.058612 


3 




2594.5008260 


F = 4.5 


F = 3.5,4.5,5.5 






23.8 








55 


38542.888064 


3 


i 


2594.5123064 


F = 5.5 


F = 4.5,5.5 






28.6 






2606 


54.94 


38366.230(3) 





9 


2606.45886(20) 




3dS4p z'^Pg 


is 




0.198 


986(150) 






55 


38366579688 


3 




2606.4351009 


F = 0.5, 1.5 


F = 0.5,1.5,2.5 






14.3 








55 


38366.467831 


3 




2606.4426999 


F = 2.5 


F = 1.5,2.5,3.5 






14.3 








55 


38366.330202 


3 




2606.4520498 


F = 3.5 


F = 2.5,3.5,4.5 






19.1 








55 


38366154185 


3 




2606.4640078 


F = 4.5 


F = 3.5,4.5 






23.8 








55 


38365.943000 


3 


i 


2606.4783550 


F = 5.5 


F = 4.5 






28.6 




Fe II 


1 AftQ 

16Uo 


JJ.OJ 


621 / L.tZy{3) 





7 

t 


160o.4DU6y/(/o) 


ou 4s aa JJ9y'2 


3d54s4py6p°/2 


h 


/.o/, 16.16 


0.0577 


— iaUU(^oUUj 






58 


62171.673196 


3 


k 


1608.4495536 








0.3 








57 


62171.652492 


3 


k 


1608.4500892 










2.1 








56 


62171.631049 


3 


k 


1608.4506440 










91.8 








54 


62171.585779 


3 


k 


1608.4518152 










5.8 






1 1 
101 1 


JJ.OJ 


62U6j.j32(j) 





i 


1 6 i 1 .zUOz^y ( /is) 




3d64py4F°/2 


h 




O.UUljo 


ilUU(oUUj 






58 


62065.503440 


3 


k 


1611.2009805 








0.3 








57 


62065.516819 


3 


k 


1611.2006332 










2.1 








56 


62065.530676 


3 


k 


1611.2002735 










91.8 








54 


62065.559929 


3 


k 


1611.1995141 










5.8 






2260 


55.85 


44232.534(6) 





9 


2260.77936(31) 




3d64p z4f°^2 


33 




0.00244 


1435(150) 




2344 


55.85 


42658.243(2) 





9 


2344.21282(11) 




3d64p z6p°/2 


ji 




0.114 


1210(150) 






58 


42658.217800 


3 


k 


2344.2142020 








0.3 








57 


42658.229605 


3 


k 


2344.2135533 










2.1 








56 


42658.241832 


3 


k 


2344.2128814 










91.8 








54 


42658.267643 


3 


k 


2344.2114630 










5.8 






2374 


55.85 


42114.836(2) 





9 


2374.46015(11) 




3d64p z'3Fg/2 


35 




0.0313 


1590(150) 






58 


42114.804727 


3 


k 


2374.4619178 








0.3 








57 


42114.819377 


3 


k 


2374.4610918 










2.1 








56 


42114.834550 


3 


k 


2374.4602364 










91.8 








54 


42114.866583 


3 


k 


2374.4584303 










5.8 
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Table Bl - continued. Atomic data for transitions usable in many-multiplet or alkali-doublet analyses 



Ion 


Tran. 


A 


uio [cm 1] 


X 


Ref. 


Ao [A] 


Lower state 


Upper state 


ID 


IP-,IP+ [eV] 


/or% 


q [cm ^] 


Fell 


2382 


55.85 


41968.065(2) 





9 


2382.76413(11) 




3d64p z«F°^,2 


ie 




0.320 


1460(150) 






5S 


41968.040382 


3 


k 


2382.7655304 








0.3 








57 


41968.051914 


3 


k 


2382.7648756 










2.1 








56 


41968.063859 


3 


k 


2382.7641975 










91.8 








54 


41968.089075 


3 


k 


2382.7627658 










5.8 






2586 


55.85 


38660.052(2) 





9 


2586.64939(13) 




3d''4p z'^D°^2 






0.0691 


1490(150) 






58 


38660.025896 


3 


k 


2586.6511386 








0.3 








57 


38660.038124 


3 


k 


2586.6503204 










2.1 








56 


38660.050790 


3 


k 


2586.6494730 










91.8 








54 


38660.077528 


3 


k 


2586.6476840 










5.8 






2600 


55.85 


38458.991(2) 





9 


2600.17222(14) 




3d''4p z6d°^2 


js 




0.239 


1330(150) 






58 


38458.965068 


3 


k 


2600.1739730 








0.3 








57 


38458.977216 


3 


k 


2600.1731517 










2.1 








56 


38458.989798 


3 


k 


2600.1723011 










91.8 








54 


38459.016359 


3 


k 


2600.1705053 










5.8 




Nill 


1709 


58.69 


58493.075(4) 





m 


1709.60409(12) 


3d9 2D5/2 


3d^4p z2f°^2 


ki 


7.64, 18.17 


0.0324 


-20(250) 




1741 


58.69 


57420.017(4) 





m 


1741.55295(12) 




3d84p z2Dg^2 


k2 




0.0427 


-1400(250) 




1751 


58.69 


57080.377(4) 
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Figure CI. Many-muldplet fit for the z = 0.452 absorber toward J000344-232355. 
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Figure C2. Many-muldplet fit for the z = 0.949 absorber toward J000344-232355. 



